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ABSTRACT 


The  objective  of  this  study  was  to  describe  the  mechanics  of  wind  wave 
generation  and  spectral  development.  Intermittency ,  high  frequency  micro- 
structure  in  wind  and  wave  fields,  and  strong  nonlinear  coupling  involving 
a  wide  range  of  scales  are  shown  to  be  crucial  elements  in  the  transfer  of 
momentum  to,  from,  and  within  the  wave  field.  None  of  these  elements  are 
included  in  available  theories. 

Measurements  of  wave  height  and  of  the  turbulent  atmospheric  and  sub- 
surface boundary  layers  were  made  from  a  small  surface  following  platform 
and  from  a  stable  38.5m  spar  buoy.  The  structure  of  moving  gust  patterns 
(cat's  paws)  is  described  and  related  to  the  generation  of  surface  waves. 
Results  from  this  and  other  background  studies  are  then  applied  to  a  dis- 
cussion of  spectral  growth  during  a  two  day  period  of  active  wave 
generation. 

Cat's  paws  contain  'bursts'  of  intense  turbulent  stress  and  buoyancy 
fluctuations  separated  by  quiescent  'intervals'.   There  is  a  difference  of 
over  three  orders  of  magnitude  in  fluctuation  strength  between  these 
features.   Rapid  growth  rate  generation  of  high  frequency  surface  waves 
and  atmospheric  turbulence  occurs  during  the  bursts.   The  resultant  micro- 
scale  components  aid  the  growth  of  lower  frequency  instabilities  by  strong 
nonlinear  coupling  between  scales  of  motion  and  by  acting  as  drag  or 
roughness  elements.   Evidence  of  strong  coupling  between  frequency  bands 
and  of  weakly  resonant  capillary-gravity  wave  interactions  is  presented. 
Thermal  stratification  has  a  strong  influence  on  fluctuation  magnitude 
and  can  delay  the  onset  of  surface  wave  generation. 

Major  spectral  growth  is  highly  unsteady.   Much  of  the  momentum  flux 
from  air  to  sea  occurs  during  intermittent  events  that  are  similar  in 
nature  to  cat's  paws,  and  goes  directly  into  high  frequency  waves.   The 
bursts  occur  predominantly  over  large  groups  of  surface  waves  and  involve 
strong  nonlinear  interactions  between  media  and  frequency  bands.   The 
long-term  equilibrium  balance  between  wind  and  water  is  disrupted  by 
variations  in  surface  currents.  There  are  'critical'  wind  speeds  charac- 
terized by  anomalous  relationships  between  parameters  of  the  predominantly 
logarithmic  velocity  profile. 
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I.   INTRODUCTION 
Description  of  the  Study 

This  paper  describes  an  experimental  study  of  the  mechanics  of  wind- 
wave  generation  and  the  spectral  development  of  sea  surface  waves.  Em- 
phasis is  placed  on  the  time  history  of  generation  and  development,  and  on 
the  nature  and  time  scales  of  air-sea  coupling.   The  basic  objective  was 
to  describe  and  illustrate  the  physics  involved  in  the  transfer  of 
momentum  to,  from,  and  within  the  wave  field.   The  underlying  hypothesis 
is  that  this  process  is  highly  nonlinear  and  depends  strongly-  on  the 
microstructure  of  the  wave  and  wind  fields.   The  experiments  were  con- 
sequently focused  upon  investigation  of  intermittent  interactions  between 
variables  and  scales  of  motion. 

Two  separate  experiments  are  described.   The  first  is  a  detailed 
investigation  of  the  structure  of  moving  gusts.  The  second  is  an  anal- 
ysis of  spectral  growth  during  a  two  day  period  of  increasing  wind. 
Both  separately  and  together  their  results  indicate  that  wind-wave  gen- 
eration and  spectral  development  involve  strong  nonlinear  interactions 
with  the  wind  and  current  fields;   representation  of  these  processes  in 
terms  of  linear  or  weakly  nonlinear  coupling  is  inadequate. 

Background 

The  available  linear  and  weakly  nonlinear  wave  generation  models 
fail  to  yield  satisfactory  estimates  of  wave  growth.  In  addition  there 
exist  unexplained  anomalies  such  as  spectral  overshoot  and  the  degener- 
ative effect  of  surface  slicks  on  long  waves.  This  combination  of  qual- 
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itative  and  quantitative  unpredictability  indicates  that  the  physics  of 
wave  generation  and  spectral  development  is  still  inadequately  under- 
stood.  Until  the  correct  physics  is  described,  there  can  be  little  sig- 
nificant progress  in  either  theoretical  modeling  or  forecasting. 

Since  1957,  most  wind-wave  generation  experiments  have  been  designed 
around  the  linear  Miles-Phillips  models  (Miles  [1957],  Phillips  [1957, 
1966]).   Early  studies  compared  observed  spectral  growth  with  theoretical 
values;  they  soon  showed  that  spectral  shape  and  initial  growth  rates 
were  reasonably  well  predicted,  but  that  growth  rates  of  the  major  por- 
tion of  the  spectrum  were  severely  underestimated.   Consequently,  more 
recent  work  in  the  low  frequency  range  has  concentrated  on  the  direct 
relationship  between  wind  and  wave  fields  in  attempts  to  determine  the 
adequacy  of  the  theory's  representation  of  momentum  transfer  mechan- 
isms.  The  high  frequency  range  (f  >lhz)  has  been  treated  as  a  separate 
entity,  primarily  relevant  to  initial  generation. 

Work  based  on  the  simplified  linear  model  has  essentially  dealt 
with  the  general  or  time-averaged  characteristics  of  wind-wave  coupling 
and  spectral  development.   It  has  neglected  the  relationship  between 
portions  of  the  wave  spectrum,  and  in  particular  has  deemphasized  the  role 
of  high  frequency  wavelets  and  the  coupling  between  waves  and  turbulence. 
There  has  thus  been  very  little  work  encompassing  the  full  range  of  spec- 
tral growth  and  relating  it  to  interactions  within  and  among  the  wind, 
wave,  and  surface  current  fields.   Hasselmann's  Hamiltonian  type  model 
(1967)  was  so  general  (and  unfamiliar  to  field  workers)  that  it  also 
failed  to  inspire  experimental  investigation  along  these  lines  (with 
the  exception  of  some  ongoing  parts  of  JONSWAP) . 

A  second  experimental  approach  to  the  problem  of  spectral  growth 
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has  been  through  parameterization  of  wind  and  wave  fields  and-  the  devel- 
opment of  spectral  transfer  functions  for  use  in  forecasting.   Such  studies 
have  described  the  generally  logarithmic  nature  of  the  atmospheric 
boundary  layer  over  water,  specified  the  characteristics  of  wave  spectra 
in  some  detail,  indicated  the  significance  of  energy  transfer  between 
wave  components,  hinted  at  the  importance  of  high  frequency  elements, 
and  yielded  partially  valid  prediction  schemes.   They  provide  little 
direct  information  about  the  physics  of  momentum  transfer  processes,  but 
rely  upon  theory  and  experience   for  guidance  in  the  choice  of  empirical- 
ly applied  parameters. 

It  is  highly  unlikely  that  either  forseeable  modifications  to  the 
available  models  or  additional  parametric  type  studies  will  correct  the 
known  theoretical  deficiences.   Further  advances  will  require  a  more 
complete  understanding  of  the  physics  of  generation.   Important  elements 
in  a  reevaluation  of  the  wave  generation  process  must  include  the  detailed 
structure  of  the  turbulent  boundary  layers  and  nonlinear  couplings,  both 
weak  and  strong,  within  and  among  the  three  fields  of  motion.   Theoreti- 
cal and  experimental  studies  have  shown  that  such  couplings  do  exist. 
Turbulence  is  known  to  be  intermittent  and  strongly  dependent  on  boun- 
dary conditions ;  a  large  number  of  scales  are  involved  in  the  dynamics  of 
flow  at  typical  geophysical  Reynolds   numbers;  and  the  generation  of  waves 
is  inseparable  from  the  generation  of  turbulence.   These  considerations 
imply  that  wave  and  turbulence  generation  may  be  a  strongly  nonlinear  pro- 
cess involving  a  large  variety  of  forces  and  scales,  and  that  the  evalu- 
ation of  experimental  results  in  terms  of  mean  values  and  long  averaging 
times  can  easily  obscure  the  physics.   Viewed  in  these  terms  the  time  his- 
tory of  events  during  wave  generation  becomes  an  important  factor;  the 
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micros tructure  of  the  air  and  water  boundary  layers  and  the  wave  field 
reassumes  great  significance;  and  the  ability  to  distinguish  the  differ- 
ences between  generation,  transfer,  and  detritus  becomes  a  major  experi- 
mental goal. 

The  pair  of  experiments  described  herein  was  designed  to  investigate 
the  waves  and  near-naviface  boundary  layers  and  to  determine  the  nature 
and  time  scales  of  the  important  processes  in  wave  generation  and  spectral 
development.   Events  ranging  in  duration  from  less  than  one  second  to 
two  days  were  examined  to  obtain  answers  to  the  following  questions: 

1.  What  are  the  important  time  scales  in  the  transfer  of 
momentum  to,  from  and  among  surface  waves? 

2.  What  are  the  nature  and  magnitude  of  the  fluctuations 
involved  in  these  transfers? 

3.  What  is  the  nature  of  the  coupling  between  media,  and 
between  scales  of  motion  both  within  and  between  media? 

4.  What  is  the  'imprint'  of  wave  generation  in  the  air  and 
in  the  water;  and  can  this  imprint  be  used  to  distinguish 
events  of  importance  in  the  generation  process  (i.e.,  both 
to  predict  and  to  permit  conditional  sampling)? 

5.  Is  there  a  similarity  between  scales,  and  if  so,  for  what 
variables  and  over  what  range  of  scales  does  the  similarity 
exist? 

The  first  experiment  was  designed  to  isolate  and  intensively  analyze 
an  event  of  known  significance  -  the  'cat's  paw'.   The  cat's  paw  is  an 
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isolated  parcel  of  atmospheric  turbulence.   It  contains  anomalously  large 
velocity  and  temperature  fluctuations,  and  can  be  visually  delineated  by 
the  'dark',  moving  patch  of  high  frequency  waves  which  it  creates.   The 
second  experiment  yielded  data  from  a  larger  system.   It  was  designed  to 
extend  the  time  scales  and  number  of  forces  involved  to  values  typical 
of  open  ocean  spectral  development.   The  experiments  are  first  individually 
discussed  in  separate  papers;  then  the  above  questions  are  answered  and 
the  development  of  a  wave  field  is  described  in  terms  of  the  results.   An 
outline  of  experimental  details  is  contained  in  appendices;  complete  des- 
criptions and  all  significant  data  are  available  in  additional  reports  of 
the  M.I.T.  Buoy  Project. 
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II.   THE  STRUCTURE  OF  MOVING  GUST  PATTERNS  (CAT'S  PAWS) 
AND  THEIR  ROLE  IN  WAVE  GENERATION 


Abstract 


Measurements  of  wave  height,  turbulent  velocity,  and  temperature 
fluctuations  within  1  meter  of  the  sea  surface  were  taken  from  a  small 
surface  following  platform.   The  resultant  data  show  that  moving  gusts 
are  highly  structured,  consisting  of  bursts  of  turbulent  stress  separ- 
ated by  relatively  quiescent  intervals.   Magnitudes  of  stress  and  buoy- 
ancy fluctuations  vary  between  bursts  and  intervals  by  over  3  orders  of 
magnitude.   The  bursts  are  associated  with  the  generation  of  very  high 
frequency  wavelets  which  in  turn  serve  to  couple  the  wave  field  to  the 
mean  air  flow.   A  large  range  of  scales  of  motion  interact,  leading  to 
a  lognormal  distribution  of  turbulent  velocity  fluctuations  and  to  ener- 
gy transfer  among  wave  components. 


Introduction 

Gust  patterns  moving  over  the  sea  surface  (cat's  paws)  are  well- 
known  to  sailors,  both  those  seeking  a  strong  wind  and  those  trying  not 
to  capsize.   Cat's  paws  can  be  located  by  their  imprint  on  the  water 
surface,  namely  the  change  in  surf ace  reflectivity  caused  by  the  high  fre- 
quency wavelets  generated.   Associated  with  these  gusts  are  large  magni- 
tude, high  frequency,  turbulent  stress  and  buoyancy  fluctuations.    Cat's 
paws  are  patches  of  intense  air-sea  interaction  with  correspondingly  large 
transfers  of  momentum,  heat,  and  moisture.   This  study  reports  an  attempt 
to  measure,  describe,  and  find  relationships  among  the  fluctuations  with- 
in a  cat's  paw,  with  particular  reference  to  the  mechanisms  involved  in 
momentum  transfer  and  wave  growth. 

Incorporated  in  a  cat's  paw  are  three  elements  of  wave  generation 
and  spectral  development  that  are  neglected  in  the  available  modern  models 
(Jeffreys  [1925,  26],  Phillips  [1957,  1966],  Miles  [1957],  Hasselmann 
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[1967]).   These  are  intermittency  of  generation,'  the  role  of  wavelets  in 
spectral  development,  and  direct,  strong  nonlinear  coupling  between  scales 
of  motion  and  between  the  sea  surface  and  atmospheric  turbulence.   The  cat's 
paw  is  itself  an  intermittent  or  burst  type  phenomenon;  the  activity  with- 
in it  is  qualitatively  and  quantitatively  different  than  that  observed 
on  the  average.   As  an  experimental  objective,  a  gust  serves  as  a  realiza- 
tion one  type  of  wave  generation  event,  with  constraints  determined  by  the 
environment,  at  Reynolds  numbers  unattainable  in  the  laboratory. 

Although  the  above  three  elements  have  been  largely  neglected,  sub- 
stantial eveidence  of  their  significance  is  available  in  the  literature. 
Munk  noted  in  1955  that  while  high  frequency  waves  account  for  an  insig- 
nificant portion  of  total  wave  energy,  they  play  a  disporportionate  role 
in  both  wind  stress  and  the  development  of  the  spectrum.   This  had  been 
exemplified  by  Van  Dorn's  (1953)  experiment,  where  the  presence  of  a  sur- 
face slick  significantly  reduced  setup  in  an  artificial  pond.   This 
effect  of  slicks  was  recently  reemphasized  (Ruggles  [1969],  Barger,  et  al, 
[1970])  by  a  demonstration  that  removal  of  wavelets  resulted  in  a  re- 
duction of  wave  energy  accompanied  by  an  increase  in  mean  wind,  decrease 
of  roughness  length,  and  significant  reduction  of  coherence  between 
wind  and  waves . 

Stewart  (1961,  1967)  has  long  maintained  that  much  if  not  most  of 
the  momentum  transfer  from  air  to  water  is  effected  by  wave  drag.   Van 
Dorn  (1953)  and  Wu  (1969  a,  1969  b,  1970)  attribute  this  wave  or  form 
drag  to  the  high  frequency  wavelets  superimposed  on  larger  waves.   Wu 
further  relates  the  wavelets  to  roughness  height  and  states  that  the  mo- 
mentum flux  to  wavelets  accounts  for  a  substantial  fraction  of  the  wind 
stress;  the  stress  appears  as  a  normal  component  to  the  wavelets,  but 
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feels  like  a  tangential  stress  to  the  underlying  waves.   Longuet-Higgins 
(1969a)  has  shown  that  such  a  fluctuating  tangential  stress,  occurring 
selectively  near  the  crests  and  leading  faces  of  longer  waves  where  the 
shorter  ones  are  'trapped'  (Munk  [1955],  Cox  [1958]),  is  dynamically  equi- 
valent to  a  normal  pressure  fluctuation  in  quadrature  with  them. 

The  role  of  turbulence  and  turbulent  Reynolds  stress  in  wave  gen- 
eration is  commonly  viewed  as  minor  or  negligible.   Miles-Phillips  theory 
holds  that  the  waves  do  not  affect  turbulence;  Hasselmann  includes  a  weak 
interaction  through  a  disturbance  of  mean  flow  by  one  wave  component, 
interaction  of  the  disturbance  with  a  turbulence  component,  and  a  feed- 
back of  the  resulting  pressure  perturbation  into  another  wave  element. 
Stewart  (1967)  has  shown  that  wave  energy  is  typically  comparable  with 
the  energy  of  turbulence  in  the  entire  troposphere,  and  consequently  has 
suggested  that  while  turbulence  can  serve  as  a  catalyst  for  transfer  of 
momentum  from  the  mean  flow,  it  can  by  itself  add  very  little  to  the  wave 
energy . 

Two  of  the  more  important  properties  of  atmospheric  turbulence  with 
respect  to  wave  generation  are  intermittency  and  the  associated  large 
growth  rates  of  localized  high  frequency  instabilities  (Greenspan  and 
Benney  [1963]).   Laboratory  experiments  (e.g.,  Klebanoff,  et  alj  [1962], 
Kline,  et  al,  [1967])  have  shown  that  the  progression  of  events  leading 
to  instability  and  transition  of  laminar  boundary  layers  involves  suc- 
cessively the  formation  of  travelling  Tollmien-Schlichting  waves ,  develop- 
ment of  spanwise  disturbances  leading  to  formation  of  local  shear  layers, 
and  explosive  development  of  the  flow  into  short  turbulent  length  scales 
in  bursts  of  high  frequency  instability  fluctuations.   The  turbulent 
boundary  layer  is  itself  intermittent,  with  periods  of  intense  activity 
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separated  by  relatively  quiescent  intervals  (e.g.,  Frenkiel  &  Klebanoff 
[1971],  Rao  and  Narayanan  [1971]).   Reverse  transition  is  similarly 
abrupt  (Steiner  [1971])  and  is  due  to  the  restoring  action  of  buoyancy 
forces.   Furthermore,  Klebanoff,  et  al,  (1962)  demonstrated  that  the 
spatial  structure  of  turbulent  breakdown  is  strongjy  dependent  on  boundary 
irregularities,  and  Gupta,  et  al,  (1968)  showed  that  unstable  oscilla- 
tions iirhced  by  the  simultaneous  introduction  of  Tollmien-Schlichting 
waves  and  free  surface  waves  resulted  in  a  large,  unpredicted  amplification 
rate  of  small  wavelengths. 

The  combined  effects  of  intermittency ,  high  growth  rates,  and  boun- 
dary influence  were  illustrated  (Mollo-Christensen  [1970])  by  filtering, 
squaring,  and  plotting  selected  frequency  components  from  wave  gauge  and 
anemometer  signals  recorded  during  a  period  of  wave  generation.   The  plots 
showed  that  ripples  and  high  turbulent  wind  energy  were  associated  with 
groups  of  large  amplitude,  low  frequency  waves.   This  suggests  that 
generation  of  waves  and  turbulence  are  coupled,  intermittent,  and  asso- 
ciated with  the  presence  of  wavelets  on  longer  wave  crests.   The  basic 
roles  of  the  high  frequency  waves  are  the  coupling  of  wind  and  wave 
fields  and  the  enhancement  of  low  frequency  instabilities;  the  wavelets 
are  generated  by  bursts  of  high  frequency  turbulence,  and  growth  rates 
in  these  local  instabilities  are  extremely  rapid  and  due  to  strong  non- 
linear coupling.   In  turn  the  turbulence  patterns  are  controlled  by  the 
major  wave  field.   The  combined  result  is  a  strongly  coupled  system  with 
growth  rates  much  larger  than  those  predicted  by  linear  or  weak  coupling 
theory. 

Implicit  in  this  outline  of  spectral  development  is  the  participa- 
tion of  a  large  range  of  scales,  coupled  by  nonlinear  interactions.   The 
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high  frequency  waves  are  in  the  capillary-gravity  and  capillary  wave 
ranges.   As  pointed  out  by  Stewart  and  Manton  (1971),  the  important 
time  scale  for  generation  by  turbulent  pressure  fluctuations  (Phillips ' 
resonance  mechanism)  is  that  of  the  wave  group  generated,  i.e.,  the  wave 
frequency.   Similarly,  the  large  amplification  rates  in  shear  layer  in- 
stability (Greenspan  &  Benney, [1963] )  have  a  typical  time  scale  of  a 
half  period  of  the  primary  wave.   Water  surface  waves  of  the  capillary- 
gravity  type  have  second  order  resonant  interactions  (McGoldrick  [1965], 
Simmons  [1969])  with  effective  time  scales  of  two  to  three  waves  periods. 
Turbulent  bursts  or  spots  are  similarly  of  short  duration.   It  can 
therefore  be  anticipated  that  the  generation  of  wavelets  is  a  short- 
lived but  intensive  phenomenon  which  is  capable  of  rapidly  depleting  the 
locally  available  energy  of  the  mean  flow.   At  the  same  time,  as  pointed 
out  by  Mollo-Chriotensen  (1971),  "work  done  on  the  water  to  make  waves  must 
equal  a  stress  times  a  strain,  and  the  reactions  of  the  stresses  on  the 
water  act  on  the  air  to  produce  turbulence." 

The  intermittency  and  selectivity  of  the  location  of  high  frequency 
waves  on  the  underlying  long  ones  is  only  partially  attributable  to  the 
coupling  with  atmospheric  turbulence.'  Variations  in  surface  tension  near 
the  crests  of  steep  waves  produce  a  local  increase  in  normal  stress  which 
in  turn  generates  capillary  waves  (Longuet-Higgins , [1963] ,  Crapper  [1970]). 
The  short  waves  steepen  on  the  forward  faces  of  longer  waves  due  to  hor- 
izontal contraction;  this  leads  to  breaking  and  the  transfer  of  mass  and 
momentum  from  the  short  to  the  longer  waves  (Longuet-Higgins  [1969b]). 
Hasselmann  (1971)  has  shown  that  the  work  done  on  the  long  waves  by  the 
"interaction  stresses"  is  balanced  to  first  order  by  the  loss  of  poten- 
tial energy  involved  in  mass  transfer,  and  to  second  order  results  in 
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a  slight  attenuation  of  the  longer  waves.   The  basic  physical  idea  sug- 
gested by  Longuet-Higgins  has  not  been  negated  by  Hasselmann's  rescaling 
of  the  terms;  there  does  exist  a  considerable  energy  interchange  between 
scales  of  motion  within  the  water,  and  this  interchange  is  strongly  depen- 
dent on  the  presence  and  location  of  high  frequency  waves. 

A  similar  interplay  occurs  in  the  atmosphere.   Mollo-Christensen 
(1970,  1971)  has  shown  that  when  the  velocity  field  is  described  in  terms 
of  three  components  -  mean  flow,  long  wavelength  and  short  wavelength 
fluctuations  -  the  resultant  energy  equations  involve  interaction  or  pro- 
duction terms,  in  the  form  of  stress  times  strain,  among  all  three  com- 
ponents.  If  the  interactions  occur  with  the  proper  phase  relationship, 
a  primary  instability  can  be  enhanced  by  the  work  done  by  Reynolds' 
stresses  of  the  secondary  instability,  and  both  instabilities  can  for 
a  while  gain  energy  from  the  mean  flow.   This  type  of  interaction  can 
occur  similataneously  on  many  scales,  in  both  the  wind  and  the  wave 
fields.   The  essential  points  are  that  Reynolds  stresses  are  inseparable 
from  the  generation  of  fluctuations,  and  that  by  nonlinear  couplings  between 
scales,  the  growth  of  short  waves  can  enhance  the  growth  of  longer  ones. 

On  the  basis  of  these  considerations  from  the  literature,  it  is 
suggested  that  the  major  reason  for  the  inability  of  available  theories 
to  satisfactorily  predict  spectral  growth  is  that  they  fail  to  include 
the  correct  physics  of  the  momentum  transfer  process.   The  inclusion  of 
strong  nonlinear  couplings  tremendously  increases  the  complexity  of  the 
problem  (  and  makes  it  as  yet  mathematically  intractable) ,  but  appears 
to  be  necessary  for  a  proper  understanding  of  wave  generation. 

The  cat's  paw  is  a  readily  accessible  phenomenon  that  provides  a 
realistic  geophysical  framework  within  which  to  initially  investigate  these 


16. 


ideas.   The  wind  tunnel  is  too  small  for  studies  at  the  Reynolds  num- 
bers and  hence  the  range  of  scales  involved  in  the  proposed  mechanisms, 
and  the  dynamic  range  of  wave  height  on  the  open  ocean  makes  detailed 
near-naviface  measurements  exceedingly  difficult.   This  experiment  was 
designed  to  investigate  the  time  history  of  cat's  paws  and  to  determine 
the  scales  involved  in  them.   The  primary  features  sought  were  the  rel- 
ative magnitudes  of  the  fluctuations  in  air  and  water  and  an  indication 
of  coupling  between  media.   Thus  the  variables  of  major  concern  were 
turbulent  Reynolds  stress  and  buoyancy  fluctuations,  and  high  frequency 
wave  activity.   Interactions  between  scales  of  motion  were  examined  in 
terms  of  the  basic  cat's  paw  structure,  turbulent  velocity  fluctuations, 
and  nonlinear  wave  interactions.   Specifically,  answers  were  sought  to  the 
following  questions: 

Is  there  a  demonstrable  relationship  between  the  generation 
of  turbulence  and  the  generation  of  high  frequency  waves?   Is  the 
process  intermittent  or  is  it  typically  near  the  mean  or  average 
flux? 

What  is  the  evidence  for  coupling  between  scales  of  motion? 

On  what  time  scales  does  instability-associated  high  growth 
rate  generation  occur? 

What  is  the  imprint  of  a  generation  'event'  in  the  air  and 
water? 
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The  Experiment 

The  measurements  were  taken  from  a  small  (1  m  long)  surface  following 
catamaran  (cat;  Figure  1).   The  cat  was  free  in  pitch,  roll,  and  heave 
(design  natural  frequency  1.5  hz)  and  restrained  in  yaw  and  horizontal 
position.   It  could  be  oriented  into  the  wind  through  the  full  range  of 
acceptable  fetch  directions.   The  experiment  was  conducted  in  March  1971 
from  the  dock  of  the  M.I.T.  Yacht  Club  on  the  Charles  River,  Cambridge, 
Massachusetts. 

Instruments  included  6,  DISA  type  55A81  hot  film  anemometers,  a  capac- 
itance wave  gauge  with  .076  cm  diameter  sensor,  and  VECO  42A70  (air) 
and  Fenwall  K496  (water)  thermistors  balanced  against  precision  resis- 
tors in  a  differential  bridge  circuit.   Subsidiary  equipment  included 
an  accelerometer  for  sensing  catamaran  motion  and  cup  and  propeller  anem- 
ometeB  for  determining  mean  wind.   Data  were  recorded  on  magnetic  tape 
in  analog  mode,  accompanied  by  a  voice  narration  that  was  later  used 
for  subjective  conditional  digitization. 

Two  of  the  hot  films  were  orthogonally  mounted  to  permit  determina- 
tion of  horizontal  (  \J   )  and  vertical  ( W  )  velocities  and  hence  Reynolds 
stress  (-/»o'w  ).   The  sensors  were  attached  to  the  forward  vertical  cat 
support  to  minimize  tilt  (-^l0  mean,  <  3°  fluctuating).   Table  1  lists  the 
instrumental  arrangement  and  mean  conditions  for  the  seven  2-hour  runs 
that  were  analyzed.   Data  from  runs  2,  4,  and  6b,  which  represent  the 
full  range  of  conditions  encountered,  were  selected  for  presentation 
herein. 

The  data  were  analyzed  using  a  combination  of  digital  and  analog 
methods.  Pre-analysis,  demodulation,  and  filtering  (100  hz  low  pass) 
were  analog.   The  data  were  then  digitized  at  250  hz  and  recorded  on  7 
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track  digital  magnetic  tape.   Final  analysis  was  accomplished  using  lo- 
cally generated  software  on  the  M.I.T.  IBM  360.   All  available  channels 
were  analyzed  simultaneously  for  each  run  segment  of  interest.   Exten- 
sive use  was  made  of  graphical  output  from  a  Stromberg-Carlson  4020 
plotter. 

The  major  contraints  on  the  analysis  were  the  cat  motion  and  the 
nature  of  the  cat's  paw  itself.   Cat  motion  could  be  successfully  removed 
from  the  data,  but  the  very  fact  that  it  did  follow  the  low  frequency 
waves  prevented  meaningful  quantitative  estimates  of  spectral  growth. 
The  cat's  paw  will  be  shown  to  consist  of  patches  of  strong  fluctuations 
separated  by  relatively  quiescent  intervals.   Thus  time  averaging  and 
the  selection  of  data  segments  for  analysis  had  to  be  very  carefully 
done  to  avoid  distortion  of  fluctuation  magnitudes  and  to  prevent  ob- 
scuring events  of  importance.   The  time  history  of  events  proved  to  be 
of  crucial  importance.   Further,  the  momentum  transfer  processes  involve 
both  internal  and  external  interactions;  they  are  highly  non-stationary, 
non-homogeneous,  and  anisotropic,  so  that  the  standard  statistical  theo- 
ries of  turbulence  are  inapplicable  and  time  averaging  could  not  be  sub- 
stituted for  ensemble  averaging.   Therefore  data  from  any  one  gust  pos- 
sesses little  if  any  statistical  significance,  and  conclusions  could  be 
drawn  only  on  the  basis  of  the  change  with  time  in  the  individual  prop- 
erties of  a  variable  and  in  the  relationship  between  variables,  and  then 
tested  on  the  basis  of  similarity  among  a  number  of  identically  treated 
events.   Each  cat's  paw,  and  each  segment  within  a  cat's  paw,  must  be 
viewed  as  a  single  realization.   No  meaningful  error  bars  or  confidence 
limits  can  yet  be  placed  on  the  data  since  a  statistically  insufficient 
number  of  events  has  been  evaluated.   The  data  were  subjectively  viewed 
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in  terms  of  the  relative  changes  with  time,  the  sequence  or  pattern  of 
events,  and  the  relationship  between  variables  in  both  time  and  frequency 
domains.  Reported  magnitudes  can  be  considered  only  as  order  of  magnitude 
estimates.   This  is  particularly  true  of  wave  autospectra,  where  filtering 
and  the  necessarily  very  short  averaging  times  preclude  even  rough  estimates 
of  the  magnitude  of  spectral  growth  and  decay. 

Emphasis  was  placed  on  understanding  the  intermittent  nature  of  tur- 
bulent fluctuations,  their  interaction  with  the  free  surface,  and  changes 
in  the  presence  and  structure  of  capillary-gravity  and  pure  capillary 
waves.  These  features  were  then  evaluated  in  terms  of  their  significance 
to  the  generation  of  surface  waves. 

Results  and  Discussion  .... 

A.   Cat's  Paw  Structure  in  Time  and  Frequency.    Separating  velocity 
and  temperature  fields  into  steady  (U,T)  and  fluctuating  (u',T')  components, 
the  equation  for  the  total  kinetic  energy  of  turbulence  can  be  written 
(Lumley  and  Panofsky  [1964]) 

®        ®  © 

where  Einstein's  summation  convention  holds,  a     is  dynamic  viscosity,  and 

T0  the  temperature  of  the  undisturbed  atmosphere.  (R\    is  commonly  called 
the  production  term;  it  represents  the  turbulent  energy  produced  by  fluc- 
tuations detracting  energy  from  the  mean  motion  by  instability  mechanisms 
as  work=stress  x  strain.   Cb)   represents  work  done  against  buoyancy  in 

terms  of  heat  flux,  H.  =  p  C  T'u!  :  an  upward  heat  flux  acts  as  a  source, 

x       p     1 

a  downward  flux  a  sink  of  kinetic  energy.   (c)   is  the  rate  of  working 
against  viscous  forces  and  represents  dissipation  of  energy  to  heat  ( €    ). 
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The  term  -/?u'w'  from  A  is  the  Reynolds  stress  (  c  );  it  commonly 
represents  momentum  flux  from  the  air  to  the  sea  surface  (downward  flux 
for  "*  >   0).   This  transfer  across  the  boundary  acts  as  an  additional  sink 
term  as  seen  from  the  viewpoint  of  atmospheric  turbulence.   Similarly, 
it  is  apparent  that  the  reaction  stresses  caused  by  the  generation  of 
waves  acts  as  a  source  term  for  the  generation  of  turbulence.   Thus  gen- 
eration of  waves  and  generation  of  turbulence  are  inseparable. 

Horizontal  (U)  and  vertical  (W)  velocities,  termperature  (T) ,  wave 
height  ( V    ),  Reynolds  stress  (*"  =  -/^u'w'),  and  buoyancy  (/^  =  w'T')  were 
measured  or  computed  and  related  in  time  and  frequency  domains  to  depict 
cat's  paw  structure.   All  data  were  taken  at  a  single  horizontal  location, 
and  w'  and  T  were  measured  at  only  one  height  in  air  and  water;  it  was 
therefore  impossible  to  determine  a  meaningful  quantitative  energy  bal- 
ance, although  the  qualitative  interplay  of  production  and  stability  terms 
was  apparent . 

A.l.    Individual  Variables.   Figures  2  through  7  (Figs.  2  - 
7)  depict  the  imprint  of  three  gusts.   Fig.  2,  from  run  2,  shows  the 
fluctuations  in  tn    ,  U,  T  ,r  ,  and  jS    for  40  seconds  (4,  10-second  time  per- 
iods [T.P.])  starting  just  in  advance  of  the  leading  edge  of  a  gust  and 
ending  aft  of  the  trailing  adge.   Fig.  3  shows  samples  of  the  computer 
generated  (SC-4020)  plots  from  which  Fig.  2  was  drawn  (W  is  also  shown). 
(The  arrows  and  slant  numerals  on  the  r  plot  of  Fig.  2  indicate  the  2- 
second  and  1-second  segments  chosen  for  detailed  analysis.)   Fig.  4  and 

Fig.  5  are  similar  representations  of  a  cat's  paw  from  run  4.  «7  ,U,  W,  x   , 

2 
a  switch  or  intermit tency  function  of  t     ^  100  dynes /cm  ,  T,  T     ,  and 

water 

-r  are  shown.   The  three  non-consecutive  T.P.  are  from  before  the 

*-  water 

gust,  near  the  leading  edge,  and  about  5  seconds  later,  still  within  the 
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gust.   Fig.  6  and  Fig.  7  are  similar  diagrams  from  run  6B.   In  run  6B 
there  are  three  consecutive  10-second  T.P.  starting  just  at  the  leading 
edge,  and  one  T.P.  from  10  seconds  after  the  gust.   These  gusts  were 
selected  to  illustrate  the  full  range  of  velocity  and  bulk  stability 
values  encountered  (Table  1) . 

The  basic  features  in  all  three  runs  areas  follows.   The  leading 
edge  of  a  gust  is  marked  by  rapid  increases  in  both  mean  and  turbulent 
wind  velocities.   Within  a  very  short  time  (<  1  second),  the  air  temper- 
ature drops  as  the  boundary  layer  is  mixed  and  destabilized  by  the  tur- 
bulent stresses.   T  continues  to  drop  as  the  winds  increase;  it  starts 
to  increase  after  U  has  reached  its  maximum  value.   There  is  commonly 
a  secondary  increase  in  U  near  the  trailing  edge  that  causes  an  additional 
drop  of  T.   In  run  6B,  this  second  T  drop  is  more  rapid  than  the  first, 
since  the  boundary  layer  has  not  yet  recovered  its  initial  static  sta- 
bility and  the  stress  fluctuations  are  consequently  more  effective  in 
mixing.   After  gust  passage,  T  rapidly  rises  to  its  initial  value  as  the 
boundary  layer  restabilizes.   There  is  no  clear-cut  relationship  be- 
tween atmospheric  variables  and  the  low  frequency  wave  field  or  water 
temperature. 

Both  ~C   andyS  traces  are  highly  intermittent.   There  are  bursts  of 
high  frequency,  large  magnitude  fluctuations,  separated  by  quieter 
intervals.   Further,  there  are  well-defined  sub-bursts  and  intervals 
within  the  larger  bursts.   The  ^  fluctuations  are  a  maximum  during  in- 
creasing wind,  but  are  not  restricted  to  these  periods.   They  are  strong- 
est and  contain  most  high  frequency  components  not  at  the  gust  leading 
edge,  but  rather  in  the  central  bursts,  after  the  boundary  layer  has 

been  mixed  and  near  the  point  of  U    .   Bursts  in  t  are  generally  asso- 
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ciated  with  similar  but  shorter  and  weaker  events  in  /&     ;  such 
coupled  events  are  followed  by  T  decreases.   Bursts  in  ft   also  occur 
independently  or  with  relatively  stronger  magnitude  than  the  t:   fluctua- 
tions.  These  events  are  followed  by  a  levelling  off  or  rise  in  T.   The 
stress  dominated  fluctuations  destabilize  and  mix  the  boundary  layer,  and 
buoyancy  forces  act  rapidly  to  restratify  it.   The  buoyancy  forces  are 
relatively  stronger  during  run  6B,  with  high  bulk  stability.   Here  the 
sub-bursts  are  clearly  evident.   Comparison  of  the  three  runs,  and  of 
the  two  T.P.  from  run  6B,  indicates  that  there  are  both  stronger  and 
higher  frequency  *  and  /3   fluctuations  during  higher  winds.   In  the  stable 
run  6B,  /&     fluctuations  die  out  before  those  in  "*  ;  in  nearly  neutral 
run  2,  even  small  c  anomalies  cause  variations  in  /S>    .    ^  water  fluctu- 
ations are  also  seen  to  be  intermittent.   They  appear  to  occur  during 

atmospheric  bursts  and  under  breaking  waves.   T      fluctuations  are 

°  water 

very  small,  and  strongly  coupled  to  the  wave  field  (i.e.,  they  are  due 

to  motion  of  the  thermistor  through  the  vertical  stratification);  a 

small  mean  T      decrease  results  from  mixing  during  the  passage  of 

the  gust.   (The  spikes  on  the  T      traces  of  Fig.  5  are  due  to 
°  r  water  & 

least  bit    noise  in  the  digitizer.) 

The  overall  pattern  indicates  that  fluctuations  in  the  leading  edge 
of  the  gust  cause  initial  destabilization  of  the  stratified  boundary 
layer.   This  is  followed  by  intermittent  large  magnitude  instability 
fluctuations.  '  Buoyancy  forces  start  to  restabilize  the  boundary  layer 
during  quiescent  intervals,  and  fully  restore  stability  within  10-15 

seconds  after  the  gust  has  passed. 

2 
The  stress  mean  values, <r>foW  ,  are  on  the  order  of  101  dynes/cm 

for  run  2,  5  x  10°  for  run  4,  and  10°  for  run  6B.  <y?>lOSec  is  approximately 
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_IQ        cm-  °C/sec  for  run  2,  and  +10   for  runs  4  and  6B.   These  values 
are  consistent  with  the  differences  in  mean  conditions  and  comparable 
with  the  commonly  assumed  values  used  in  dynamic  modeling  of  the  exchange 

processes  (z-  =  1  dyne/cm2,  H  =  3x10   gm-cal/ cm-sec) .   The  'instantaneous' 

3        2 
(250  hz)  values  commonly  approach  or  exceed  f  -~  10  dynes/cm  and  /S^ 

2 
3x10  cm-°C/sec  in  the  bursts.   The  strong  high  frequency  fluctuations 

thus  have  amplitudes  typically  three  orders  of  magnitude  greater  than 

the  mean,  or  departures  from  the  mean  of  10  to  50  standard  deviations. 

Similar  types  of  measurements  taken  5.6  to  22.6  m  over  a  wheat  field 

during  unstable  stratification  (20  hz  sampling  rate)  (Haugen,  et  al,  [1971]) 

2         2         2 
yielded  burst  magnitudes  of  c  ~  10  dynes/cm  , fi  ~    3x10   cm  °C/sec.   Burst 

duration  was  roughly  10  seconds,  and  <   U>  ~  600  cm/sec.   This  data  also 
showed  the  presence  of  bursts  within  bursts.   McBean  (1970)  also  found 
similar  values  over  land,  with  a  strong  effect  of  stability  (momentum 
transfer  efficiency  decreased  and  heat  transfer  efficiency  increased  with 
decreasing  stability).   The  major  qualitative  patterns  of  intermittent 
gusts  are  thus  similar  over  land  and  water.   Quantitatively,  the  time 
interval  selected  for  averaging  is  deterministic  in  estimates  of  instan- 
taneous fluctuation  magnitude  (Section  C) .   A  15-minute  averaging  period 
for  the  cat's  paw  data  would  have  increased  the  above  estimates  by  at 
least  an  order  of  magnitude.   The  difference  in  magnitude  between  cat 
and  land  results  is  an  indication  of  the  significance  of  stability  and 
the  enhancement  of  frequency  and  magnitude  of  stress  fluctuations  caused 
by  the  participation  of  a  fluctuating  boundary  capable  of  supporting 
well-defined  high  frequency  waves. 

Fig.  8  shows  some  typical  histograms  of  the  individual  variables 
from  the  10  second  T.P.  of  all  three  runs.   Fluctuation  magnitudes  and 
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variance  are  much,  larger  during  than  before  or  after  the  gusts;  the  peak- 
edness  of  the  r  ,  yS  ,  and  w  histograms  is  due  to  the  high  frequency 
fluctuations  and  consequent  large  number  of  axis  crossings.   U  and  T' 
distributions  are  more  unimodal  after  than  during  the  gust.   There  are 
few  other  basic  differences  between  the  T.P.'s  for  any  one  variable  when 
the  relatively  long  10-second  intervals  are  used. 

Fig.  9  shows  typical  autospectra  based  on  the  10  second  T.P.  (6  degrees 
of  freedom;  90%  confidence  limits  are  .63  x  10    to  1.74  x  10  for  a  spec- 
tral density  estimate  of  10  ).   With  this  comparatively  long  averaging 
time  there  are  only  generalized  differences  in  frequency  content  between 
gusts  and  mean  conditions.   With  the  exception  of  the  wave  gauge  data, 
the  spectra  are  fairly  white,  with  (note  run  6B  data)  a  slightly  more 
rapid  drop-off  with  increasing  frequency  in  the  T.P.  after  the  gust.   As 
also  shown  in  the  time  plots ,  the  high  frequencies  are  augmented  during 
the  bursts.   The  differences  in  magnitude  between  gust  and  no  gust  con- 
ditions are  apparent;  the  peaks  of  y   (f)  and  §>**  (f)  are  3  decades 
lower  after  the  gust.   This  basic  similarity  of  pattern,  but  difference 
in  high  frequency  content  and  major  difference  in  intensity  reflects  the 
patterns  displayed  in  Fig.  7  where  the  t    and  /3    scales  have  been  deter- 
mined by  the  computer  on  the  basis  of  the  10-second  extrema.   Wave  propa- 
gation and  the  presence  of  long  period  wave  groups  prohibit  any  deductions 
from  magnitudes  of  the  wave  spectra. 

Fig.  10  displays  the  changes  in  the  0.5  second  vertical  profiles  of 
U  from  run  6B.   It  also  shows  <C  U  >  0.5  sec  and  a  measure  of  the  turbu- 
lent kinetic  energy  production  term.   Profile  perturbations,  the  devel- 
opment of  unstable  shear  layers,  and  turbulent  production  occur  predom- 
inantly during  bursts.   As  with  the  T    and^  fluctuations,  production  has 
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a  range  of  i.  3  to  4  orders  of  magnitude  and  is  highly  intermittent.   The 
comparatively  long  averaging  times  (0.5  sec)  are  incapable  of  showing 
the  extremely  rapid  fluctuations  within  the  bursts;  instead,  the  profile 
changes  are  similar  to  those  that  would  be  produced  by  a  set  of  eddies 
propagating  past  the  sensor  array.   Again,  the  maximum  turbulence  pro- 
duction occurs  after  the  leading  edge  in  T.P.  2,  i.e.,  after  the  atmos- 
pheric boundary  layer  has  been  destabilized.   The  profile  remains  nearly 
steady,  and  production  nearly  zero,  after  the  gust  has  passed  and  thermal 
stratification  has  been  reestablished.   Both  stability  and  intermittency 
are  of  major  importance  in  gust  time  history. 

A  sample  time  history  of  three  cat's  paws  thus  indicates  that  gusts 
are  highly  structured.   Averaging  techniques  such  as  spectra  and  histo- 
grams tend  to  obscure  the  significance  of  the  microstructure.  f  and  yS 
fluctuations  in  bursts  are  at  least  3  orders  of  magnitude  above  the 
mean.   The  localized  instabilities  also  contain  a  larger  range  of  fre- 
quency components  than  is  the  case  in  quieter  periods  and  intervals. 
Production  of  turbulence  and  mixing  in  the  boundary  layer  are  strongly 
associated  with  the  intermittent  bursts;  r  fluctuations  are  largest 
after  the  mixing  and  destabilization  associated  with  the  leading  edge. 
Participation  of  a  large  range  of  scales  of  motion  and  of  the  fluctuating 
water  surface  enhances  the  burst  magnitudes.   The  use  of  long  averaging 
periods  obscures  the  details  of  the  transfer  processes  that  are  dependent 
upon  these  strong  localized  instabilities. 

A. 2.   Relationships  between  Variables.    Relationships  were 
depicted  by  joint  histograms  (JE)  and  time  plots  in  the  time  domain,  and 
by  coherence  (Coh)  plots  in  the  frequency  domain.   This  duality  of  re- 
presentation helped  to  lend  confidence  to  conclusions  based  on  a  limited 
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amount  of  data.   Coherence  may  be  defined  aa 


It  is  represented  in  terms  of  its  magnitude  and  phase, 

as,(rt  •  V„<rt  ei6"«(0    . 

It  was  computed  using  fast  fourier  transform  techniques.   In  the  conven- 
tions adopted,  Q.J>0   means  <$    leads  "$   .  Averaging  was  over  25  adjacent 
frequency  bands  for  10  second  T.P.,  9  for  the  shorter  intervals."  As  shown 
by  McBean  (1970),  the 'spectral  correlation  coefficient'  (or  Coh)  is  an 
actual  measure  of  transfer  efficiency  as  a  function  of  frequency  (or  scale). 

a.   10  Second  Intervals.    Fig.  11  shows  JH(U,T') 
(with  10  second  mean  and  trend  both  removed  from  T)  from  run  6B.   Asso- 
ciated with  the  decrease  and  increase  of  <  T  >  10  sec  is  the  boundary  layer 
destabilization  and  mixing  by  the  leading  edge,  violent  mixing  during 
the  major  bursts,  and  subsequent  restabilization  of  a  nearly  constant  temp- 
erature after  gust  passage.   The  top  of  Fig.  12  shows  a  similar  pattern 
in  the  frequency  domain  using  CohT'„  from  run  4.   Before  the  gust  the 
waves  and  temperature  fluctuations  (above  5  hz)  are  quite  coherent  and 
nearly  out  of  phse.   This  is  the  pattern  to  be  expected  if  the  catamaran 
rides  up  and  down  through  a  thermally  stratified  layer  that  follows  the 
pattern  of  the  waves.  (Ramzy  and  Young  [1968]).   Near  the  leading  edge 

">r'r7    is  slightly  lower  and  extends  to  lower  frequencies  ;®t'o  begins  to 
shift,  particularly  in  higher  frequencies,  as  the  atmosphere  is  mixed. 
Later  in  the  gust  the  temperature  field  has  been  effectively  inverted,  and 
T  and  rj     are  more  nearly  in  phase.   The  lower  half  of  Fig.  12  shows  Cohra 
from  the  same  T.P.   Here  ^^B     ^s  mucb  lower  but  ^ip      shows  effectively 
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the  same  temperature  inversion. 

Fig.  13  and  Fig.  14  show  JK(r ,  6    )  and  Cohr^    from  run  6B;  Fig.  15 
shows  JK(  r  ,  A    )  from  run  2 .   These  diagrams  show  that  the  mixing  and 
tumbling  is  effected  by  T   fluctuations  destabilizing  the  boundary 
layer,  and  that  the  following  restabilization  is  due  to  the  predominance 
of  buoyancy  forces.   The  change  of  the  patterns  with  time  is  most  clear 
in  the  more  stable  run  6B   (note  the  change  in  scales);  in  the  nearly 
neutral  run  2  the  fluctuations  show  less  correlation,  are  more  randomly 
scattered.   Here  the  stress  fluctuations  predominate  and  remain  intermit- 
tently large  even  after  the  gust  has  passed,  while  in  run  6B  the  buoyancy 
forces  rapidly  debilitate  the  momentum  flux  as  the  gust  departs.   The 
coherence  diagrams  also  indicate  the  phase  reversals  of  stress  and  buoyan- 
cy fluctuations  as  the  boundary  layer  is  turned  and  mixed.   The  rapid,  nearly 
linear  shift  of  &x&   with  f  in  the  9.5-17  hz  range  of  T.P.  4  implies  that 
fluctuations  of  similar  frequency  are  'slipping'  past  each  other,  or 
effectively  that  in  the  presence  of  strong  stratification,  weak  ~c    fluc- 
tuations are  unable  to  overturn  and  mix  the  stable  boundary  layer. 

Boundary  layer  stability  is  effective  in  retarding  momentum  flux 
by  turbulent  velocity  fluctuations.   This  effect  is  particularly  large 
near  the  leading  edge  of  cat's  paws.   Wave  generation  during  high  fre- 
quency instabilities  will  thus  be  less  efficient  during  stable  than 
neutral  stratification,  and  delayed  or  reduced  by  the  conversion  from 
kinetic  to  potential  energy  involved  in  disruption  of  the  stability 
structure. 

Fig.  16  and  Fig.  17  show  the  JH  and  Coh  represetations  of  the  re- 
lationship between  U  at  two  heights  during  run  6B.   As  with  the  other 
variables,  fluctuations  are  largest  at  and  just  after  the  leading  edge, 
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with  a  collapse  of  the  histogram  after  the  gust  has  passed.  tfU(J  is  sig- 
nificant only  at  the  DC  level  (T.P.  1-3)  and  the  reversals  of    £>uo(D.C) 
demonstrate  the  large  variations  in  mean  wind  velocity  shear  that  occur 
during  the  gust.   The  generally  low  level  of  0U(J  and  lack  of  JH  patterns 
over  even  this  short  distance  are  indicative  of  strong  vertical  struc- 
ture and  the  anisotropic  character  of  the  high  frequency  intermittent 
fluctuations. 

Fig.  18  and  Fig.  19  are  JH(  T    tr7    )  from  run  6B  and  Coh  t„    from 
during  and  after  the  gusts  of  runs  2  and  6B  and  from  all  three  T.P.  of 
run  A.   It  is  apparent  that  there  is  no  effective  evidence  of  strong 
wind-wave  coupling  during  the  gust  when  viewed  from  this  averaging  period. 
The  effective  generation  mechanisms  are  obscured  by  the  'detritus'  or 
by  averaging  over  events  of  radically  differnt  character.   If  anything, 

~c    and  y  are  more  coherent  in  the  7-12  hz  range  (of  runs  2  and  6B) 
after  than  during  the  gust;  and  in  this  frequency  band  the  continuous 
linear  Or*       drift  is  evidence  that  the  fluctuations  of  similar  frequen- 
cy slip  by  one  another  rather  than  strongly  coupling.   The  turbulent 
stress-wave  relationship  that  will  be  demonstrated  has  again  been  obscur- 
ed by  considering  averaging  intervals  much  longer  than  the  effective  inter- 
mittent events. 

b.  2  and  1  Second  Intervals.  Two  second  segments 
containing  primarily  'burst'  or  primarily  'interval'  type  fluctuations 
were  selected  from  runs  2  and  6B.  The  location  of  these  2-second  T.P. 
is  shown  in  Fig.  2  and  Fig.  6  (with  the  exception  of  T.P.  1  from  run  6B 
which  was  taken  about  3  seconds  prior  to  arrival  of  the  leading  edge.) 
A  portion  of  the  results  is  shown  in  Fig.  20,  JH(  Z     ,  /3     )  from  run  2, 
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and  in  Fig.  21,  Cohry   from  run  6B.   These  figures  show  some  overall 
similarity  to  those  of  the  10  second  T.P.;  some  of  the  segments  still 
contain  both  bursts  and  intervals,  and  this  tends  to  partially  blur  the 
effects  of  the  very  high  frequency  fluctuations. 

Fig.  20  shows  that  the  mixing  and  destabilization  process  is  indeed 
associated  with  the  bursts.   Relatively  weak,  positively  correlated  fluc- 
tuations occur  near  the  leading  edge;  mixing  occurs  during  the  next 
strong  burst  and  is  followed  directly  by  negative  correlation  of  -c 
and  i&      and  consequent  further  destabilization  of  the  boundary.   With  the 
boundary  layer  inverted  the  buoyancy  forces  work  toward  restabilization 
and  rapidly  degrade  the  strength  of  the  momentum  fluxes.   Both  the  mixing 
and  the  restabilization  occur  rapidly;   the  intial  boundary  layer  struc- 
ture is  wiped  out  by  the  leading  edge  but  is  never  firmly  reestablished 

6. 

during  gust  passage. 

The  wind-wave  coupling  begins  to  be  apparent  on  this  time  scale. 
Fig.  21  shows  that  stress  fluctuations  are  intially  (before  the  gust)  out 
of  phase  with  the  waves  (7-15  hz  range).   They  are  almost  totally  inco- 
herent during  the  leading  edge,  i.e.,  during  destabilization.  Tri7   then 
increases  slightly  in  T.P.  3.   This  T.P.  contains  several  well-def ined  sub- 
bursts  and  intervals  so  that  the  picture  is  confused,  but  &z<?    is  relative- 
ly constant  at  the  YT^  peaks  at  5.9  and  12  hz.   During  the  major  burst, 
the  high  frequency  instabilities  and  high  frequency  waves  serve  to  strong- 
ly and  immediately  couple  the  lower  frequency  fluctuations  to  the  under- 
lying waves.   The  high  frequency  coupling  is  still  not  evident,  but  \zn     is 
high  at  and  below  the  main  wave  period,  with  ®tij  -  TT/'Z. ,  or  the  maximum 
stress  occurring  aft  of  the  crests.   In  the  immediately  following  inter- 
val *$xf?       at  tne  D»C.  end  drops  off  slightly  and  the  phase  shifts;  in 
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the  higher  frequencies  o^-  is  slightly  higher  but  the  propogating  cap- 
illary waves  now  'slip'  with  respect  to  stress  fluctuations  of  the  same 
frequency. 

In  order  to  further  decrease  the  blurring  of  burst  and  interval, 
segments  approximately  1  second  long  were  selected  from  within  the  2- 
second  T.P.  (see  Fig.  6).   The  results  from  run  6B  are  shown  in  Fig.  22- 
Fig.  24.   T.P.  1  is  prior  to  gust  arrival;  T.P.  2  just  behind  the  leading 
edge;  T.P.  3  during  a  short  but  intense  fluctuation  2  seconds  later;  T.P. 
4  during  the  major  burst  and  T.P.  5  from  the  quiescent  interval  3-4 
seconds  after  T.P.  4.   The  pattern  of  events  during  this  sequence  may 
be  described  as  follows. 

Prior  to  the  arrival  of  the  gust,  t  and  A    are  nearly  totally  co- 
herent and  in  phase.  X     fluctuations  are  small,  loosely  coupled  to  the 
main  wave  field,  and  out  of  phase  with  the  high  frequency  waves.   With 
the  arrival  of  the  leading  edge  of  the  cat's  paw  the  stress  fluctuations 
increase  in  frequency  and  magnitude  and  the  boundary  layer  is  tumbled; 
the  ~C  i  fi     fluctuations  decrease  in  coherence  and'  become  nearly  out  of 
phase;  little  effect  is  noted  immediately  on  the  water  surface,  as  the 
turbulent  kinetic  energy  is  utilized  to  stir  the  boundary;  waves  and  ~c 
are  incoherent.   After  the  intial  burst  the  air  is  mixed  and  high  fre- 
quency waves  are  generated,  associated  with  the  short,  strong  bursts  of 
momentum  flux.  T     and  A     are  again  coherent,  but  out  of  phase,  and  the 
wind-wave  coupling  is  apparently  random.   With  the  boundary  layer  de- 
stabilized, the  magnitude  of  the  stress  fluctuations  increases  (T.P.  4) 
and  the  high  frequency  instabilities  become  strongly  associated  with  the 
generation  of  surface  waves.  trg      is  quite  low  at  frequencies  above  15  hz 
with  Bis      more  variable  than  before  but  centered  around  O  .     *Zn 
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increases  considerably  in  the  <«>  -18  hz  range,  but  remains  fairly 
low  for  lower  frequencies.   ®rfc  slips  across  part  of  the  coherent  ranges 
but  shows  a  major  reversal  of  trend  in  the  center  of  the  coherent  band. 
During  the  quiescent  period  immediately  after  this  event  r  fluctuations 
are  much  smaller  and  less  strongly  coupled  to  the  short  wave  field,  with 
slightly  increased  Xz„        in  lower  frequencies,  z   f  /6>  fluctuations  are 
again  strongly  coupled,  and  out  of  phase  due  to  the  mixing  of  the  preceed- 
ing  burst.   This  process  is  repeated  during  subsequent  bursts,  until 
the  cat's  paw  has  passed  the  site  and  stability  is  reestablished. 

The  change  in  character  of  the  high  frequency  stress  and  wave  pat- 
terns (Fig.  23)  is  strong  evidence  of  the  large  growth  rates  associated 
with  intermittent,  large  magnitude,  high  frequency  instabilities.   The 
generation  of  high  frequency  waves,  as  in  T.P.  3  and  especially  in  T.P.  4, 
is  indeed  strongly  coupled  to  the  generation  of  turbulence.   The  large 
value  of  YTy      for  lower  frequencies  in  2-second  T.P.  4  indicates  that 
the  high  frequency  waves  can  rapidly  and  effectively  couple  the  lower 
frequency  waves  and  instabilities,  and  aid  their  growth.   It  is  obvious 
that  little  statistical  confidence  can  be  placed  in  low  frequency  measures 
based  on  such  short  time  periods;  nevertheless  the  change  between  time 
periods  and  averaging  lengths,  the  similarity  of  this  pattern  in  other 
bursts  and  gusts,  and  the  enhanced  high  frequency   Q™   in  1-second  T.P.  4 
lend  strong  subjective  credence  to  the  hypothesis. 

Figs.  23  and  24  also  clearly  show  how  wave  generation  is  delayed  by 
the  stability  structure.   Generation  processes  act  only  after  the  boundary 
layer  has  been  initially  destabilized.   The  continuous  Oxn     slipping  in 
T.P.  5  (across  the  full  width  of  the  coherent  band)  appears  to  be  charac- 
teristic of  'detritus',  in  this  case,  of  the  propagation  and  advection  of 
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the  waves  and  turbulence  generated  during  the  high  frequency,  localized 
instabilities.   The  slipping  pattern  is  similar  in  T.P.  4  except  for  the 
10  -  18  hz  range,  where  &  z  n       levels  off  and  reverses  direction  near  the 
zero  level.  Most  of  the  generation  and  momentum  transfer  occurs  in  this 
frequency  range  (see  Fig.  27). 

Generation  appears  to  be  a  very  strong  but  selective  and  short- 
lived phenomenon.   The  time  scales  associated  with  the  sub-bursts  and 
instability  associated  transfers  are  very  short,  typically  on  the  order 
of  1  to  1.5  seconds,  which  is  also  the  typical  period  of  the  immediately 
underlying  major  wave  field.   The  strongest  coupling  and  high  growth  rates 
associated  with  the  localized  instabilities  occur  in  a  narrow  high  fre- 
quency band.   Averaging  times  significantly  longer  than  these  scales  can 
depict  the  general  changes  in  atmospheric  structure  but  fail  to  show 
the  detailed  nature  of  the  coupling  and  transfer  processes.   The  impor- 
tant scales  are  thus  those  of  the  waves,  the  bursts,  and  the  internal 
fluctuations  within  the  bursts.   Momentum  is  transferred  to  the  water 
at  scales  associated  with  these  features. 

B.    Interaction  Between  Scales.      The  microscales  in  both  wind 
and  wave  fields  account  for  only  an  insignificant  portion  of  the  total 
energy,  yet  as  shown  above  they  are  the  active  participants  in  the  local- 
ized high  growth  rate  instabilities.   They  are  effective  in  aiding  the 
growth  rate  of  lower  frequency  components  and  in  coupling  wind  and  wave 
fields,  as  was  heuristically  explained  by  Mollo-Christensen  (1970,  1971) 
and  as  demonstrated  in  the  time  plots  and  coherence  diagrams  of  the  pre- 
ceeding  section.   Other  evidence  of  nonlinear  coupling  between  scales 
may  be  shown  by  an  examination  of  the  probability  distribution  of  turbulent 
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velocity  fluctuations  and  of  the  shape  of  the  high  frequency  wave  spectrum. 

B.l.   Scale  Interaction  in  Velocity  Fluctuations.   A  gust 
has  been  shown  to  consist  of  a  wide  range  of  sizes  of  bursts  or  eddies 
of  strong,  high  frequency  fluctuations,  separated  by  relatively  quiet  in- 
tervals.  The  cat's  paw  itself  is  merely  a  burst  on  a  larger  time  scale. 
The  absolute  magnitude  of  a  single  fluctuation  is  thus  dependent  not  only 
on  the  instantaneous  departure  from  adjacent  values,  but  also  on  the  de- 
parture of  the  fluctuation  and  its  neighbors  from  similarly  sized  parcels, 
and  on  the  departure  of  the  mean  magnitude  of  those  patches  from  segments 
on  yet  larger  scales;  and  so  on  up  the  range  of  scales  to  an  averaging 
period  covering  the  largest  perturbations  involved.   Thus  each  scale  of 
motion  involves  fluctuations  whose  magnitude  depends  both  upon  the 
strength  of  internal  instabilities,  and  upon  the  forcing  of  and  interac- 
tions between  the  scale  itself  and  larger  scales  or  larger  anomalies. 
Following  Cramer  (1955),  let  the  instantanious  magnitude  of  a  random 
variate  be    £    ;  and  let  the  change  of  this  variate  caused  by  a  for- 
cing t»   be  proportional  both  to  the  forcing  and  to  the  magnitude  3§£-\ 
which  had  in  turn  been  produced  by  the  previous  J>_,   impulses. 
Then     *.  -    t.    =  CS-  £:    c  =  coi\-*i&t±  9 


and      £.  J   =    7     2  3-^— S' 
i--i  ■  c     -'   5; 


fc 


which  in  turn,  in  the  limit  of  small  change,  "*"  4r   \   "^f    £"*•&) 


For  large  N,  the  central  limit  theorem  states  that   3-    is  normally 
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distributed.   Hence   <5;   has  the  lognormal  distribution, 

where  P  ^  (x)  is  the  probability  "§  =  x,  CS"2  is  the  standard  deviation 
of  3  >  E(x)  =  ^  §  >  ,  and  a  is  an  appropriate  normalizing  or  cut-off 
value  of   ^ 

The  lognormal  distribution  can  thus  be  expected  to  occur  when  a 
large  number  of  impulses,  or  a  large  range  of  scales,  interact  to  con- 
tribute to  the  increase  of  a  variable  at  a  rate  proportional  to  its  in- 
stantaneous value.   For  example,  Rao  and  Narayanan  (1971)  have  shown 
that  the  intervals  between  bursts  in  a  turbulent  boundary  layer  were  log- 
normally  distributed.   This  indicates  that  the  scales  of  the  bursts  and 
intervals  are  determined  by  coupling  with  the  scales  of  the  boundary  and 
of  the  outer  or  main  flow  as  suggested  by  Mollo-Christensen  (1970). 

Within  a  cat's  paw  or  gust,  one  variable  which  should  act  in  a 
lognormal  fashion  is  the  turbulent  energy  dissipation,  6  .   Dissipation 
is  primarily  dependent  upon  the  magnitude  of  the  velocity  fluctuation 
rate,  which  in  turn  depends  strongly  on  the  amount  of  'rubbing'  or  inter- 
action between  different  scale  eddies.   The  dissipation  is  due  to  the 
working  of  molecular  viscous  stress  against  this  coupled  rate  of  strain 
or  deformation  of  the  fluid  particles.   Provided  that  the  range  of  scales 
involved  in  the  deformations  is  sufficiently  large,  say  including  inter- 
action between  the  localized  high  frequency  instabilities,  fluctuations 
on  the  scale  of  bursts  and  intervals,  and  similarly  between  groups  of 
sub-bursts  and  bursts  up  to  the  scale  of  the  full  cat's  paw,  dissipation 
should  be  approximately  lognormally  distributed. 

Another  variable  which  should  be  lognormally  distributed,  and  for 
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much  the  same  reasons,  is  the  magnitude  of  the  squared  time  derivative 
of  horizontal  velocity.   Time  derivatives,  like  space  derivatives,  are  en- 
hanced by  the  interactions  between  a  variety  of  scales  of  eddies  or  in- 
stabilities.  The  interaction  between  high  frequency  velocity  fluctuations 
and  bursts  has  already  been  demonstrated  in  the  preceeding  section.   Define 
the  normalized  squared  magnitude  of  the  velocity  derivative  as 


X;  - 


<r  c  ^/*±Y  > 


Chen  (1971)  has  analyzed  the  distribution  of  3f*{'S\)      with  data  taken 
from  FLIP  during  BOMEX.   He  pointed  out  that  squaring  the  derivative 
emphasizes  the  large  fluctuations  and  minimizes  the  magnitude  of  small 
ones.   (Gurvich  and  Yaglom  [1967]  showed  that  a  lognormal  distribution 
with  large  variance  corresponds  to  a  high  probablity  of  values  close  to 
0  or  much  larger  than  the  rms,  but  small  probablities  of  intermediate 
values.)   The  natural  logarithm  of  such  a  variable  will  have   %  Cxj  — <=»<=° 
as   "5  — *  O   due  to  the   (^°/^0  *°        'tail'  imposed  by  the  rapid 
axis  crossing  of      /  9 1        .   Taking  the  average  of  a  small  number  of 
values  (i.e.,  ™-\Z   ;t«,  "^i  J  m  «=^a'1      )  wipes  out  this  tail.   In  physical 
terms  this  corresponds  to  averaging  the  fluctuations  over  a  very  small 
but  finite  volume,  or  numerically  to  limiting  the  analysis  to  scales 
slightly  larger  than  the  cut-off  imposed  by  the  sampling  rate.   This  aver- 
aging will  increase  the  tendency  toward  lognormality,  provided  that 
the  number  of  scales  interacting  is  sufficiently  large  and  that  the  averag- 
ing time  is  short  enough  to  include  the  important  structure.   Chen's 
data  were  from  a  roughly  4-minute  period  at  mean  heights  of  3  and  31  meters 
above  the  sea  surface.   The  resulting  distribution  was  almost  identically 
lognormal  for  averages  over  8  data  points. 
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Similar  results  from  runs  2  and  6B  of  this  study  are  shown  in  Figs. 
25  and  26.   The  burst/no  burst/both  character  of  the  T.P.  is  noted  in 
the  figures.   Fig.  25,  with  the  data  in  histogram  form  (run  2,  2-second 
T.P.),  clearly  shows  how  averaging  (over  5  points)  removes  the  'tail'  of 
the  distribution.  In  Fig.  26,  the  cumulative  probabilities  of  JL\   (  ^) 
from  both  runs,  2-and  1-second  T.P. ,  are  plotted  on  a  normal  probability 
scale.   Purely  lognormal  distributions  with  comparable  means  and  vari- 
ances are  indicated  by  straight  lines. 

Since  only  a  relatively  small  number  of  points  was  used,  the  central 
limit  theorem  can  not  be  strictly  applied  and  none  of  the  plots  are 
strictly  lognormal,  though  most  of  them  show  relatively  minor  departures. 
As  shown  best  by  the  1-second  periods  from  the  low  stability  run  2,  per- 
iods of  intense  activity  show  a  deficit  of  intermediate  values,  an  excess 
* 

of  large  and  small  ones.   The  quiet  periods  show  the  opposite  effect, 
particularly  in  T.P.  6,  taken  after  the  passage  of  the  gust  (past  the 
end  of  Fig.  2).   In  this  very  quiet  T.P.  the  fluctuations  are  more  nearly 
normally  than  lognormally  distributed.   The  2-second  T.P.  show  these 
effects  less  clearly  due  to  averaging  over  both  bursts  and  intervals, 
but  the  tendency  toward  normality  is  clear  in  the  quiescent  T.P.  5.   The 
run  6B  data  (high  bulk  stability)  shows  similar  results,  with  intermediate 
deficits  during  bursts,  excesses  during  intervals.   One  major  exception 
is  T.P.  2,  from  the  leading  edge  of  the  cat's  paw.   It  was  shown  in  the 
preceeding  section  that  in  this  time  period  the  major  effect  of  the 
stress  fluctuations  was  the  local  tumbling  and  destabilization  of  the 
boundary  layer.   The  excess  of  intermediate  departures  in  this  distribu- 
tion is  further  evidence  of  the  initial  effectiveness  of  the  buoyancy 
forces  in  retarding  momentum  transfer  and  limiting  fluctuation  strength. 
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The  distribution  of  the  squared  velocity  derivative  is  thus  an  ef- 
fective indicator  of  the  range  of  scales  involved  in  boundary  layer  tur- 
bulence.  Interactions  are  strongest  and  involve  the  largest  number  of 
scales  during  the  bursts.   During  quiescent  intervals  and  particularly 
between  gusts,  fewer  scales  interact;  coupling  is  weak  or  non-existant , 
and  Jw t%i)      tends  to  have  a  normal  distribution,  indicative  of  a  lack  of 
interaction. 

It  should  be  noted  that  for  isotropic,  homogeneous  turbulence, £ 
is  related  directly  to  dissipation  by  use  of  Taylor's  frozen  turbulence 
hypothesis,   -S-   =  -  "U  <L         .   Chen  in  fact  regarded  his  data  in 
these  terms.   Ruggles  (1969)  has  demonstrated  that  the  frozen  turbulence 
hypothesis  is  quantitatively  inapplicable  in  a  case  even  less  restrictive 
than  this  study;  and  indeed  the  above  results  indicate  that  it  is  not 
necessary  to  so  restrict  the  assumptions   on  the  nature  of  the  flow  in  or- 
der to  account  for  the  lognormal  behavior  of  fluctuation  magnitude.   The 
well-known  intermittency  of  the  small  scale  structure  of  turbulence,  the 
consequent  large  dispersion  of  £  ,  and  the  dependence  of  the  probability 
structure  of  €        upon  larger  scales  (Obukhov  [1962],  Landau  [1959]) 
have  been  theoretically  applied  to  the  statistics  of  turbulence  only  for 
the  locally  isotropic,  homogeneous  case  (Kolmogorov  [1962]).   The  present 
development  has  avoided  these  restrictions  and  emphasized  strong  non- 
linear coupling  between  scales  as  the  causative  phenomenon.   Further, 
it  has  been  demonstrated  that  local  departures  from  lognormality  are 
also  indicative  of  the  relative  strength  of  coupling  of  the  turbulent 
components,  and  that  many  features  of  the  boundary  layer  and  transfer 
processes,  such  as  the  spacing  of  bursts  and  intervals  and  the  frequency 
and  strength  of  local  velocity  fluctuations,  are  controlled  by  this 
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coupling. 

B.2.   An  Example  of  Weak  Resonant  Wave  Interactions. 
Nonlinear  wave-wave  interactions  have  been  theoretically  studied  in  con- 
siderable detail  during  the  last  10  years  (Kinsman  [1965]  provides  a 
good  summary  of  the  major  concepts).   Energy  exchange  among  gravity 
waves  of  different  wave  numbers  is  only  seen  when  the  perturbation  ana- 
lysis extends  to  third  order;  it  involves  triads  of  primary  waves  inter- 
acting to  feed  energy  into  a  fouth  component.   For  waves  in  the  cap- 
illary-gravity range,  the  analysis  shows  that  interaction  occurs  at 
second  order  (as  well  as  at  third  order);  two  waves  interact  to  feed 
energy  into  a  third  (McGoldrick  [1965],  Simmons  [1969]).   This  weak 
resonant  interaction  has  been  demonstrated  in  the  laboratory  (McGoldrick 
[1970]).        » 

McGoldrick  (1965)  has  discussed  the  implications  of  this  energy 
transfer  with  respect  to  wind-wave  generation.   Of  particular  importance 
is  the  fact  that  the  interaction  times  involved  in  the  growth  of  the  re- 
sonant wave  are  only  2-3  periods  of  the  initial  waves  creating  it.   Be- 
cause of  this  short  interaction  time,  the  amount  of  energy  that  the  wave 
field  can  receive  from  the  wind  at  a  given  frequency  is  determined  not 
just  by  breaking  but  also  by  nonlinear  transfer  to  other  components. 

Of  particular  interest  is  the  'self -resonant'  wave,  of  wave  number 
k,  =  J?/iT     (  T   =  surface  tension/density).   This  is  commonly 
known  as  Wilton's  wave;  he  noted  the  presence  of  a  singularity  at  this 
wave  number  in  his  (1915)  perturbation  analysis.   The  general  require- 
ments for  resonance  at  2nd  order  are  k-+  k~  +  k~  =  0,  w1  ±  w»  +  w~  =  0, 

where  *  =  *.vf  a*  J  k  =  **/»  *  ^ 


and  £•  >  3  kj  +rT  W*  \~-  \v,-z 
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k^  fulfills  this  relationship  for  k  =  k_  =  -1/2  k?,   or   k2  -    \  ■—= 

Simmons  noted  that  waves  shorter  than  k_  are  hydrodynamically  unstable 
at  2nd  order;  further,  k~  is  the  longest  wave  which  can  interact  with 
two  longer  waves.   Waves  at  k^  are  extremely  unstable  and  can  occur  only 
as  a  transient  state;  they  can  not  continuously  exist  as  free  surface 
waves  and  once  generated  will  rapidly  transfer  their  energy  to  the  k. 


'2 


• — i  -        3    2 
waves.   For  clean  water,  JL  =   75  cm  /sec  and  k^.  and  k„  roughly 


correspond  to  frequencies  of  10  and  20  hz,  respectively. 

Fig.  27  shows  the  wave  autospectra  from  the  2-second  segment  ana- 
lysis of  run  6B  (80%  confidence  limits  are  roughly  1  decade  wide).   All 
5  T.P.  display  the  peak  of  the  major  wave  field  at  around  1  hz.   T.P.  1, 
from  before  the  gust,  is  rather  smooth  in  the  band  above  3  hz.   T.P.  2, 
near  the  leading,  edge,  shows  a  slight  deficit  around  9.5  hz,  followed  by 
a  levelling  off  to  1A  hz,  a  rapid  drop  to  19  hz,  and  a  slight  increase 
around  20  hz.   During  the  'active'  T.P.  3  and  the  major  burst  in  T.P.  4, 
spectral  slope  is  less  steep  and  the  spectrum  shows  a  slight  bulge  near 
10  hz  (i.e.,  near  the  frequency  where   ®r  r?         leveled  off  and  reversed 
slope);  in  T.P.  4,  there  is  also  a  local  maximum  just  above  20  hz.   In 
T.P.  5,  immediately  after  the  burst,  there  is  a  significant  deficit 
near  10  hz,  and  again  a  local  minimum  near  19  and  maximum  just  above 
20  hz.   This  sequence  of  events,  particularly  the  change  between  T.P.  4 
and  5,  is  strongly  suggestive  of  the  picture  given  by  McGoldrick.   Dur- 
ing generation  the  k^  waves  are  not  only  present  but  near  saturation 
because  of  the  energy  input  from  the  wind;  as  soon  as  the  forcing  dis- 
appears, so  do  the  extremely  unstable  k..  waves.   Their  energy  is  'drained' 
to  higher  frequencies.   Fig.  28  shows  a  similar  example  from  run  2, 
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where  the  mean  wave  field  was  significantly  larger  and  the  results  con- 
sequently less  clear.   Here  the  major  burst  was  in  T.P.  3,  with  lesser 
activity  in  T.P.  4  and  quiescent  in  T.P.  5.   In  T.P.  5  the  spectral 
deficits  occur  at  10  and  18  hz,  with  the  local  maximum  around  19.5  hz. 
The  slight  frequency  shifts  of  the  postulated  k.  and  k0  waves  are 
readily  attributable  to  contamination  of  the  Charles  River.   The  very 
fact  that  the  predicted  effects  of  weak  resonant  interactions  stand  out 
at  all  in  other  than  a  'carefully  tuned'  laboratory  experiment  is 
strong  evidence  of  their  importance  in  wave  spectral  development.   Ev- 
idence that  perturbation  analysis  is  suitable  for  predicting  continuum 
nonlinear  effects  in  major  wave  fields  was  recently  provided  by  the 
bispectral  analyses  of  Garrett  (1970).   Nonlinear  transfer  of  energy 

is  thus  apparently  of  importance  at  all  stages  of  spectral  development. 
I 

C.   A  Further  Note  on  Time  Averaging.    The  results  of  the  preceed- 
ing  sections  have  emphasized  the  importance  of  high  frequency  fluctuations, 
the  intermittent  nature  of  generation,  and  the  coupling  between  scales. 
Averages  over  time  scales  comparable  to  the  entire  gust  or  to  a  combina- 
tion of  burst  and  interval  have  been  shown  to  obscure  the  selective 
high  growth  rates  associated  with  the  local  instabilities,  the  effec- 
tiveness of  wavelets  in  coupling  the  wind  and  the  waves,  and  the  nature 
of  interactions  between  scales.   One  further  example  will  serve  to  il- 
lustrate the  critical  importance  of  proper  averaging  in'  the  analysis  of 
wind-wave  generation  data. 

Consider  a  typical  cat's  paw  from  run  2.   Typical  mean  velocity 
values  for  different  scales  can  be  very  roughly  taken  as 
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Similarly,  temperature  means  may  be  estimated,  with   ^\;r  <ea  y  £  c   ,  as 


^T  > 


Next,  assuming  rms  turbulence  levels  of  10%  (a  very  low  estimate  for 

a  gust  or  burst),  3.75°  <:  T\    <  4.25°,   18  m/sec  <    u',    <      22  m/sec. 
6  '  burst  burst 

As  demonstrated  in  Section  A,  the  important  fluctuations  in  heat 
and  momentum  transfer  and  wave  and  turbulence  generation  are   /S  =  w'T' 
and  t   -  -p  u'w'  (assuming  symmetry  in  v'  and  V  =  0) .   Then 

I 
U  -  <U>  +  u*,     W  =  <^  W  >  +  w',    T=^T>+T' 

and  since 

«£  W  >  and<^w'>are  taken  =  0  over  whatever  averaging  period  is  chosen, 


u'w'  =  w'(U  -^U>  )     and   w'T'  =  w'(T  -  <:  T  >  )  . 

Because  of  the  impressed  zero  mean  of  W,  the  long  term  mean  fluxes 

will  be  independent  of  the  choice  of  < U >  and  <T>   ;  however,  both  the 

magnitude  and  the  sense  (sign)  of  the  instantaneous  fluxes  (or  of  very 

short  averages)  are  strongly  affected  by  this  choice.   For  the  example 

above,  choosing  U  =  <U>     ,  T  =  <T>     ,  and  plotting  a  joint  histo- 
b  gust  gust  t.    J 

gram  of  the  resultant  values  of  ZT  and  /3       for  the  burst,  the  result 
will  appear  as  in  Fig.  29B  (Fig.  29  is  from  run  2,  2-second  T.P.  2);  while 
if  the  means  are  chosen  as  the  burst  means,  results  are  as  shown  in  Fig.  29A. 
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When  viewed  against  the  general  gust  background  it  appears  as  if  the  f 
and  fi      fluctuations  are  negatively  correlated  (of  opposing  effect  in 
turbulence  production)  and  very  strong;  when  viewed  against  the  2-second 
background,  mostly  burst,  the  f luctutations  are  relatively  much  smaller 
and  show  a  weak  positive  correlation.   Similarly,  when  the  results  are 
applied  to  the  turbulent  kinetic  energy  balance,  the  effects  on  production, 
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Figure  29 
and  hence  on  the  stability  structure  of  the  boundary  layer,  are  opposite. 
This  is  an  extrememly  straightforward  but  seldom  applied  result.   The 
proper  choice  of  averaging  period  must  be  made  on  the  basis  of  the  im- 
portant scales  in  the  process  under  consideration.   For  wave  generation 
and  air-sea  momentum  transfer,  these  have  been  shown  to  be  the  scales  of 
and  within  the  bursts  and  the  main  wave  field.   Consideration  of  the  wave 
generation  process  in  terms  of  long  term  mean  effects  can  totally  obscure 
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the  physical  realities 


Summary  of  Results  and  Conclusions 

The  primary  facts  to  emerge  from  this  study  may  be  summarized  as  fol- 
lows, 

1.  Cat's  paws  are  highly  intermittent  phenomena  that  them- 
selves are  composed  of  localized  bursts  separated  by  quiescent  intervals. 
The  bursts  contain  high  frequency  fluctuations  of  stress  and  buoyancy 
3-4  orders  of  magnitude  larger  than  that  observed  in  the  mean.   The  var- 
iability in  time  was  shown  in  Fig.  2  -  Fig.  7;  representative  magnitudes 
are  shown  there  and  in  Fig.  8. 

2.  The  bursts  are  associated  with  the  generation  of  high 
frequency  waves  and  turbulence.   This  is  best  demonstrated  by  the  se- 
quence of  events  shown  in  Fig.  23.   Participation  of  the  fluctuating 
boundary  enhances  the  magnitude  of  the  fluctuations ,  as  shown  by  a  com- 
parison of  the  results  of  this  study  with  measurements  over  land. 

3.  The  high  frequency  instabilities  not  only  have  large  growth 
rates  on  their  internal  scales,  but  also  are  associated  with  the  growth  of 
lower  frequncy  components.   Both  high  and  low  frequency  fluctuations  are 
strongly  enhanced  in  the  bursts,  as  shown  by  the  changes  in  magnitude  and 
shape  of  autospectra  (Fig.  9)  and  on  the  detailed  SC-4020  plots.   This  is 
not  only  true  for  couplings  within  but  also  for  couplings  between  media  as 
demonstrated  by  the  diffeiEnces  in  coherence  between  stress  and  surface 
waves  in  the  2-  and  1-second  time  periods  (especially  T.P.  4,  Figs.  21  and 
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24).   The  bursts  are  also  associated  with  the  development  of  shear  layers 
in  the  wind  profile  and  with  large  values  of  turbulence  production  (see 
Fig.  10),  a  further  indication  of  the  importance  of  localized  instabilities. 

4.  The  boundary  layer  is  mixed  and  tunbled  by  the  instabil- 
ity fluctuations  in  the  bursts.   The  effect  of  strong  temperature  strati- 
fication in  the  boundary  layer  is  to  delay  the  onset  of  wave  generation  by 
requiring  intial  utilization  of  available  turbulent  kinetic  energy  for 
destabilization.   Buoyancy  forces  are  very  effective  in  degrading  the 
magnitude  and  effectiveness  of  momentum  flux,  and  quickly  cause  res tab il- 
ization  after  passage  of  the  gust. 

The  general  mixing  and  boundary  layer  inversion  patterns  are  shown 
in  the  joint  histograms  of  U  -  T'  (Fig.  11)  and  r - fB>       (Figs.  13,  15)  and 
in  the  phase  reversals  in  the  coherence  of  T'  -  7    and  "c  -p     (Figs.  12,  14). 
The  mixing  associated  with  the  bursts,  and  the  partial  restabilization 
during  intervals,  is  shown  by  the  changes  in  correlation  between  ~c     and 
/3   in  Fig.  20.   The  delay  of  generation  attributable  to  the  utilization 
of  kinetic  energy  for  mixing  is  demonstrated  in  Figs.  22  -  24,  where  the 
burst  at  the  leading  edge  (T.P.  2),  uncoupled  to  the  waves,  has  little 
effect  upon  the  water  surface  but  effectively  inverts  the  boundary  layer. 

The  effects  of  stability  have  also  been  demonstrated  by  the  ex- 
cess of  intermediate  sized  velocity  fluctuations  in  the  leading  edge 
(Fig.  26,  runs  T.P.  2),  by  the  rapid  rise  of  T  after  gust  passage  (Figs. 
11,  15),  and  by  the  temperature  increases  and  decreases  associated  with 
Z     and  ft     fluctuations  of  different  relative  strength  (e.g.,  Fig.  6). 

5.  The  coupling  between  scales  is  strongly  nonlinear,  as 
illustrated  by  the  lognormality  of  the  distribution  of  squared  velocity 
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fluctuations  during  bursts  (Figs.  25,  26),  the  strong  association  of  high 
frequency  instabilities  with  large  magnitude  velocity,  stress,  and  buoy- 
ancy perturbations  (Figs.  2-7  and  9),  and  the  association  of  burst  length 
with  the  period  of  the  major  underlying  wave  components  (Figs.  2,  4,  6). 

6.  Weakly  resonant  nonlinear  interactions  are  effective  in 

transferring  energy  among  capillary-gravity  waves.  This  was  demonstrated 

by  the  change  in  the  k.  and  k„  components  (Figs.  27,  28)  between* 'genera- 
te     2. 

tion'  (burst)  and  'detritus'  conditions. 

7.  The  important  momentum  transfer  processes  are  of  short 
duration,  highly  nonlinear,  frequency  selective,  and  strongly  associated 
with  intermittent  high  frequency  instabilities  and  coupling  between  scales. 
This  was  demonstrated  by  the  variations  in  ~c-'7    coherence  patterns  between 
different  averaging  times  and  burst/no  burst  conditions  (Figs.  21,  24),  by 
the  simultaneous  occurrance  of  strong  Reynolds  stress  and  capillary  wave 
activity, (Fig.  23),  and  by  the  selectivity  of  coherence  and  spectral  change 
(Figs.  24,  27,  28). 

8.  The  differences  between  'generation'  and  'detritus'  can  be 
described  by  the  changes  in  coherence  between  variables  and  in  the  proba- 
bility distribution  of  velocity  fluctuations.   Detritus  coherence  is  char- 
acterized by  low  magnitude  (^  )  and/or  by  continuous  phase  (  ©  )  slippage 
in  any  coherent  bands;  during  strong  coupling  and  generation,  T  is  high 
and  &  is  nearly  constant  in  the  frequency  bands  that  participate  most 
strongly  in  the  transfers  (Figs.  19,  21,  24).   Velocity  fluctuations 
(measured  as  ^   )  in  detritus  have  excesses  of  intermediate  values  and 
approach  normal  distributions;   those  from  generation  periods  are  more 
nearly  lognormally  distributed  with  excesses  of  very  large  values  and 
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values  near  zero.   These  differences  are  due  to  the  enhanced'  range  of  scales 
and  strength  of  the  couplings  involved  in  wave  generation  and  turbulent 

interactions. 

2 
The  intermittency  switch  function  (I  =  1,  ~C  >,   100  dynes/cm  ,1  =  0, 

2 
t    <      100  dynes/cm  ),  shown  only  in  Fig.  5  but  computed  in  all  runs  and 

used  for  selection  of  the  2  and  1  second  time  periods,  appears  to  be  a 

suitable  first  approximation  for  a  conditional  sampling  switch.  When 

smoothed  over  a  time  period  ^0.5  seconds  it  provided  a  reliable  indicator 

of  the  presence  of  a  burst  and  instability-associated  generation. 

9.   There  is  a  rough,  qualitative  similarity  between  scales 
of  motion  (e.g.,  gust  and  burst),  events  of  different  magnitude  (Fig.  7), 
and  between  turbulence  patterns  over  land  and  water  surfaces.   The  similar- 
ity lies  in  the  presence  of  a  large  range  of  scales.  The  significant  dif- 
ferences are  caused  by  the  strength  of  coupling  between  these  scales;  the 
larger  the  number  of  participating  scales  and  media,  the  stronger  the  coupl- 
ing and  the  more  efficient  the  associated  heat  and  momentum  fluxes.  These 
features  were  demonstrated  by  a  comparison  of  flux  measurements  over  land 
and  water,  and  illustrated  in  the  time  plots  of  the  variables.  The  simul- 
taneous enhancement  of  frequency  range  and  fluctuation  magnitude  in 
bursts  is  perhaps  best  displayed  in  Fig.  7,  where  the  'instantaneous' 
magnitudes  are  based  on  the  10  second  means. 

It  was  also  demonstrated  (Section  C,  Fig.  29)  that  the  effect  of   r 
and  A     fluctuations  can  be  different  when  their  magnitude  and  sense  are 
evaluated  over  different  scales;  e.g.,  a  burst  may  contribute  to  boundary 
layer  instability  as  seen  from  the  long  term  mean,  while  the  local  and 
immediate  effect  of  the  fluctuations  is  to  temporarily  reestablish  a 
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portion  of  the  initial  temperature  stratification.  Averaging'  over  time 
scales  other  than  those  of  importance  in  the  momentum  transfer  process 
can  lead  to  misinterpretation. 

10.   Evaluation  of  generation  in  terms  of  the  'mean'  can  be 
misleading  and  obscures  the  physics  of  the  transfers.   Since  much  wave  and 
turbulence  generation  does  not  occur  on  the  mean  but  in  localized,  fre- 
quency selective,  high  growth  rate  instabilities  (as  depicted  in  Figs.  23, 
24,  and  27),  and  since  the  mean  itself  is  a  seldom  if  ever  assumed  state 
(as  shown  in  the  time  series  diagrams  [Figs.  2-7]  and  the  histograms 
[Fig.  8]),  representation  of  the  generation  mechanisms  only  in  these  terms 
is  unsatisfactory. 

On  the  basis  of  these  results  it  is  concluded  that  microscales  play 
a  predominant  role  in  wave  generation.  High  frequency  waves  and  turbulence 
are  intermittently  generated  during  high  growth  rate,  local  instabilities, 
and  the  wavelets  subsequently  act  as  drag  elements  to  couple  the  longer  waves 
to  the  wind.   Further,  the  high  frequency  instabilities  aid  the  growth  of 
similar,  lower  frequency  fluctuations.   This  combination  of  events  leads 
to  spectral  growth  rates  much  larger  than  those  predicted  on  the  basis  of 
linear  and  weakly  nonlinear  theories . 
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III.   THE  ROLE  OF  MICROSTRUCTURES  IN  WIND-WAVE  SPECTRAL  DEVELOPMENT 


Abstract 

Measurements  of  the  wind,  wave  and  surface  current  fields  were  made 
from  a  38.5  meter  spar  buoy  at  the  mouth  of  a  tidal  channel.   Statistical 
analysis  of  velocity  profiles  showed  that  both  air  and  water  boundary  lay- 
ers were  predominantly  logarithmic  (neutral  stability) ,  and  confirmed 
Ruggles'  (1970)  finding  of  anomalous  behavior  of  atmospheric  profile  para- 
meters at  critical  wind  speeds.   A  lack  of  simple  relationships  among  pro- 
file measures  indicates  that  wind-wave  coupling  can  be  only  very  poorly 
characterized  on  the  basis  of  means,  that  the  time  history  of  spectral 
evolution  is  crucial  to  an  understanding  of  the  physics  governing  the  mo- 
mentum transfer,  and  that  more  than  one  wave  generation  regime  exists. 

A  two-day  history  of  spectral  development,  which  did  not  include  the 
major  profile  anomalies,  was  analyzed  to  determine  the  nature  and  scales 
of  momentum  transfer  processes.   Wave  generation  was  found  to  occur  pre- 
dominantly in  short,  intermittent,  burst-type  events,  and  to  involve 
strong  nonlinear  coupling  of  wind  and  wave  fields  with  associated  rapid 
growth  of  the  high  frequency  components.   These  important  events  involve 
coupling  between  wave  groups  in  wind  and  surface  wave  fields  and  strong 
interactions  between  the  media  and  between  frequency  bands  in  both  fields. 
Tidal  currents  modulate  the  wave  field  and  cause  imbalance  of  wind-wave 
equilibrium  conditions. 


Introduction 

Microscale  studies  in  cat's  paws  have  demonstrated  the  importance  of 
intermittent,  strong  interactions  between  scales  of  motion  and  between  at- 
mospheric turbulence  and  the  sea  surface  in  the  generation  of  high  frequen- 
cy waves.   Similar  types  of  interactions  are  also  significant  in  the  de- 
velopment of  the  major  portion  of  the  wind-wave  spectrum. 

o 

In  the  open  ocean,  typical  Reynolds  numbers  are  on  the  order  of  10  ; 
local  wave  generation  occurs  under  the  influence  of  surface  drift,  tidal 
currents,  propagating  swell,  and  stability,  precipitation,  wind  speed 
and  direction  shifts  associated  with  large  weather  systems.   The  large 
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number  of  participating  scales  and  phenomena  can  all  interact,  with  both 
strong  and  weak  coupling,  to  have  a  significant  effect  on  spectral  de- 
velopment.  It  will  be  shown  that  there  are  separate  regimes  of  wave  gen- 
eration, and  that  the  time  history  of  events  is  of  crucial  importance  in 
development  of  the  wave  field. 

Most  previous  open  ocean  studies  of  spectral  development  have  dealt 
with  mean  or  long  time  average  wave  energy  growth  and  have  examined  the 
influence  of  only  atmospheric  conditions  upon  the  wave  field.   Airborne 
measurements  (e.g.,  Barnett  and  Wilkerson,  1967;  Schule  et  al,  1971)  and 
other  laboratory  and  oceanic  studies  which  extended  over  a  range  of  fetch 
(e.g.,  Snyder  and  Cox,  1966;  Sutherland,  1967;  Mitsuyasu,  1969)  have 
dealt  extensively  with  the  major  features  of  spectral  development  with 
little  or  no  experimental  reference  to  the  causative  momentum  transfer 
mechanisms.   From  these  types  of  studies  spectral  growth  rates,  the  fre- 
quency shift  of  the  spectral  peak,  the  existence  of  an  equilibrium  range 
for  high  frequencies,  and  the  overshoot/undershoot  phenomena  have  become 
well  known  features. 

Studies  from  stationary  buoys  or  platforms  (e.g.,  Seesholtz,  1968; 
Ruggles,  1969)  have  evaluated  growth  rates  in  terms  of  wind  velocity  pro- 
file characteristics  with  minimal  consideration  of  the  effects  of  the 
underlying  current  and  subsurface  boundary  layer  conditions.   Such  studies 
have  shown  that  the  atmospheric  boundary  layer  is  essentially  logarithmic 
(at  least  for  near-neutral  stability)  and  that  spectral  growth  rates  are 
substantially  larger  than  those  predicted  by  the  inviscid  laminar  Miles 
theory  (Miles,  1957). 

Investigations  of  the  direct  interplay  between  wind  and  waves  have 
concentrated  on  the  transfer  mechanisms  suggested  by  this  theory.   One  of 
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its  critical  elements  is  that  turbulent  and  wave  induced  Reynolds  stresses 
do  not  interact,  and  that  momentum  transfer  occurs  by  means  of  an  insta- 
bility in  the  mean  shear  flow.   The  momentum  is  extracted  from  the  mean 
flow  in  the  critical  layer  and  transmitted  downwards  by  wave  induced  Rey- 
nolds stresses. 

The  recent  work  of  Kendall  (1970)  demonstrated  that  a  fixed  wave 
pattern  caused  significant  modulation  of  the  turbulent  structure,  but 
that  the  induced  fluctuations  did  not  have  a  constant  phase  relation  to 
the  waves.   Davis  (1970)  and  Stewart  (1970)  theoretically  and  experimen- 
tally investigated  the  vertical  structure  of  Reynolds  stress. and  found 
poor  quantitative  agreement  with  the  wave-turbulence  Reynolds  stress  pro- 
file predicted  by  instability  theories.   Wave  induced  Reynolds  stresses 
accounted  for  only  about  half  of  the  momentum  carried  to  the  wave  surface 
by  the  pressure  component  in  phase  with  wave  slope. 

Ruggles  (1969) ,  Davidson  (1970) ,  and  Shemdin  and  Lai  (1970)  have  all 
demonstrated  the  existence  of  wave  induced  perturbations  in  the  wind 
field  and  have  shown  that  these  wave  induced  Reynolds  stresses  do  interact 
with  turbulent  stresses.   There  is  thus  substantial  evidence  that  inter- 
action between  atmospheric  turbulence  and  surface  waves  is  one  of  the  cru- 
cial mechanisms  involved  in  the  transfer  of  momentum  to  the  waves.   This 
problem  is  "being  studied  by  Saeger  and  Reynolds  (1971);  their  initial  re- 
sults indicate  that  the  inviscid  model,  by  neglecting  the  important  struc- 
ture of  the  critical  layer,  gives  very  poor  predictions  of  the  Reynolds 
stress  and  hence  of  momentum  flux  to  the  waves.. 

The  cat's  paw  studies  have  demonstrated  the  importance  of  detailed, 
near-surface  measurements  in  determining  the  interactions  among  forces. 
Dobson  (1969,1971)  recently  made  the  first  successful  detailed  field 
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measurements  of  atmospheric  pressure  on  wind  generated  waves.   His  im- 
portant results  showed  that  the  observed  phase  shifts  were  an  order  of 
magnitude  larger  than  predicted  by  Miles'  theory,  and  that  energy  and 
momentum  flux  spectra  were  strongly  peaked  at  or  above  the  peak  of  the 
locally  generated  waves.   Further,  the  inviscid  laminar  Miles  model 
underpredicted  increases  in  wave  energy  by  factors  of  5-8.   These  facts 
suggested  a  strong  interaction  between  wave  induced  and  turbulent  Rey- 
nolds stresses.   Dobson  also  noted  that  energy  flux  and  wave  supported 
momentum  flux  were  negligible  above  relatively  low  frequencies,  contrary 
to  an  assumption  that  much  of  the  momentum  transferred  goes  into  short 
waves  (Wu,  1969;   and  the  cat's  paw  and  present  studies)  but  not  contra- 
dicting the  suggestion  that  short  waves  aid  the  growth  of  lower  frequency 
instabilities  by  nonlinear  coupling  between  scales  of  motion.   Near  sur- 
face pressure  measurements  over  waves  made  by  Elliott  (1970)  confirmed 
Dobson' s  pressure  shift  findings.   Elliott  noted  that  pressure  lagged 
the  wave  by  about  135°  during  generating,  and  by  180°  during  non-genera- 
ting conditions.   There  was  no  vertical  phase  shift  at  heights  above  the 
crests.   Elliott's  work  also  suggested  strong  interactions  between  wave 
induced  components  and  normal  boundary  layer  turbulence. 

The  basic  conclusion  from  detailed  investigations  of  near-naviface 
structure  and  air-sea  coupling  is  that  turbulence  and  waves  interact 
strongly  and  intermittently.   Only  a  relatively  small  portion  of  the 
wave  growth  is  due  to  processes  which  act  on  the  'mean',  and  these  pro- 
cesses are  themselves  strongly  dependent  upon  the  presence  and  location 
of  high  frequency  components  and  associated  strong  nonlinear  coupling. 

The  structure  of  the  underwater  boundary  layer  is  also  of  major 
importance  in  transfer  processes.   It  has  been  repeatedly  demonstrated 
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that  the  magnitudes  as  well  as  the  mechanisms  involved  in  transfer  of  mo- 
mentum, heat,  and  moisture  across  an  interface  are  radically  different  for 
flat  plates  or  solid  surfaces  and  fluctuating  or  moving  surfaces  (e.g., 
Chambers  et  al,  1970;  Mangaralla  et  al,  1970).   Over  water,  the  processes 
are  different  for  the  cases  of  wind  waves  alone  and  wind  waves  with  the 
added  presence  of  swell  and  surface  drift. 

Few  experimental  studies  have  considered  the  effect  of  the  water 
near-surface  layer  motions  upon  wind-wave  coupling,  although  the  weak  in- 
teractions between  waves  and  currents  or  very  long  swell  are  well  known 
theoretically  (e.g.,  Longuet-Higgins  and  Stewart,  1960,  1961).   The 
available  experimental  evidence  on  the  structure  of  and  interactions  in 
the  water  boundary  layer  is  somewhat  contradictory.   Shonting  (1967, 
1968)  noted  that  u-w  phase  relations  were  close  to  those  predicted  by 
classical  theory,  with  no  phase  shift  with  depth  and  only  minor  interac- 
tion between  wind  waves  and  swell.   Shonting,  Stewart  and  Grant  (1962) 
and  Yefimov  and  Khristoforov  (1969)  have  all  found  that  wave  stress  is 
strongly  concentrated  very  near  the  surface,  and  that  essentially  all 
dissipation  occurs  above  the  trough.   The  transfer  between  high  and 
low  frequencies  caused  by  dissipation  (breaking)  and  viscous  decay  of 
short  waves  (Longuet-Higgins,  1969b;  Hasselman,  1971)  is  also  a  surface 
phenomenon,  and  involves  non-linear  coupling  between  scales.   Shemdin 
and  Lai  (1970)  found  that  the  wind-induced  drift  in  a  wind/water  tunnel 
had  a  logarithmic  profile,  and  that  the  surface  stresses  computed  from 
the  air  and  water  profile  parameters  were  roughly  equivalent.   Intro- 
duction of  mechanical  waves  disrupted  the  wind  driven  current.   The 
couplings  between  the  adjacent  boundary  layers  and  the  interface  are  thus 
as  yet  imperfectly  described. 
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This  experiment  was  designed  to  investigate  the  nature  and  time 
scales  of  air-sea  coupling  during  wave  generation.   It  was  conducted 
during  a  variety  of  wind,  weather  and  current  conditions.   The  analysis 
was  divided  into  three  segments.   First,  a  study  was  made  of  the  air  and 
water  boundary  layers  in  terms  of  a  logarithmic  profile  model.   This 
section  extends  and  complements  the  investigations  of  Ruggles  (1970) 
and  Groscup  (1971) ,  and  describes  the  general  environmental  background 
of  the  more  detailed  analyses.   Second,  the  time  history  of  spectral 
growth  during  a  two  day  period  of  increasing  winds  was  investigated.   In 
this  section  the  general  nature  of  coupling  between  wind,  wave  and  current 
fields  is  described.   Finally,  segments  of  the  two  day  record  were  studied 
in  detail  to  emphasize  the  importance  of  strong,  intermittent  generation 
of  high  frequency  waves  and  resultant  wave-turbulence  couplings  that  were 
demonstrated  in  the  cat's  paw  study.   Generation  'events'  are  described 
and  their  influence  on  the  main  wave  field  is  displayed. 

The  major  objectives  of  the  study  were  to  demonstrate  the  wide 
range  of  interactions  involved  in  wave  generation  and  spectral  develop- 
ment, and  to  emphasize  the  importance  of  high  frequency  'microscales ' 
and  intermittency  in  the  momentum  transfer  process. 


The  Experiment 

Measurements  of  wind  and  water  speed  profiles  and  wave  height  were 
taken  from  a  stable,  38.5m  spar  buoy  (Mollo-Christensen  and  Dorman,  1971) 
(Fig.  1).   The  buoy  was  slack  moored  in  31.1m  of  water  (MLLW)  and  faced 
into  the  tidal  current  (NE  flood,  SW  ebb,  1.5m/see  maximum  current). 
The  experiment  ran  from  July  through  October,  1970.   The  site  was  at 
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Figure  1.   M.I.T.  Air-Sea  Interaction  Buoy 
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the  mouth  of  a  tidal  channel  (Vineyard  Sound,  Massachusetts);  minimum 
fetch  to  the  S  and  SW  (direction  of  prevailing  winds  and  sea  breeze)  was 
50km. 

Instruments  included  5  cup  anemometers  at  heights  of  2. Am  to  11.4m, 
3  ducted-propeller  (5cm  diameter)  current  meters  at  depths  of  -1.5  to 
-4.7m,  and  two  capacitance  type  wave  gauges  with  4.6m  sensor  length.   The 
primary  wave  gauge  was  mounted  directly  over  the  current  meters,  2.5m 
cross-stream  from  the  anemometers.   Thermistors  were  attached  to  arms  of 
the  velocity  transducers  at  +5  and  -2.5m.   Buoy  tilt  was  sensed  with  a 
pendulum  pot. 

Signals  from  the  transducers  were  mixed,  telemetered  to  a  shore  sup- 
port facility  on  Cuttyhunk  Island,  Massachusetts  and  recorded  in  analog 
mode.   Time,  voice  log,  and  25kc  clock  signals  were  transmitted  and  re- 
corded for  use  in  synchronization.   Velocity  information  was  in  pulse 
format,  wave  gauge  signals  a  variable  frequency,  and  thermistor  data  a 
variable  D.C.  voltage.   Data  processing  was  done  at  M.I.T.  and  used  a 
combination  of  analog  and  digital  methods.   Digitization  rate  for  the 
time  history  analysis  was  lOhz;  velocities  were  pre-filtered  at  2.5hz, 
wave  gauge  signals  at  5hz,  to  avoid  aliasing. 

A  total  of  146,  10  minute  'mean  condition'  runs  were  made.   There 
were  5  series  of  runs,  with  data  taken  once  each  hour  for  periods  ranging 
from  6  to  53  hours.   109  of  the  runs  contained  data  suitable  for  analysis. 
Winds  and  waves  were  from  the  S-SW  during  almost  all  runs,  with  10  min- 
ute mean  wind  speeds  (at  10m)  of  1.5  to  12.5m/sec.   Because  of  the  large 
fetch  the  atmospheric  boundary  layer  was  nearly  neutrally  stable  during 
all  runs  (AT  air-sea  £   3°C  and  -.02  <  Ri(bulk^  .05).   Similarly,  the 
water  surface  layers  were  generally  well  mixed  down  to  a  depth  of  at 
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feast  8-10m. 

The  analysis  was  conducted  in  three  parts.   First,  a  distribution 
analysis  was  made  of  the  parameters  obtained  from  a  comparison  of  velocity 
data  with  a  neutral  stability  logarithmic  profile  model.   Of  the  109, 
10  minute  anemometer  runs,  96  (88%)  had  correlations  ;>>  0.94  with  the 
model.   Of  the  51  minimeter  profiles  used  (from  the  longest  nearly  con- 
tinuous segment,  runs  109-159,  also  used  in  the  time  history  analysis), 
48  or  94%  were  logarithmic  to  within  the  required  contraint;  49  of  the 
anemometer  profiles  from  the  same  period  were  logarithmic.   Much  of  the 
analysis  was  based  on  30  second  averaging  times.   For  the  anemometers 
there  were  2467  samples  (81%  logarithmic;  1057  samples  with  91%  loga- 
rithmic for  runs  109-159) ,  for  the  minimeters  1060  samples  (85%  loga- 
rithmic) .   Ruggles  (  1969,  1970)  has  discussed  the  significance  of  the 
correlation  coefficient;  Groscup  (1971)  has  described  the  increasing 
tendency  to  a  logarithmic  profile  with  increasing  averaging  time  for  the 
runs  discussed  herein.   Only  the  runs  satisfying  the  /°   *  0.94  criterion 
were  used  in  the  histogram  analysis,  since  the  parameters  are  computed 
from  the  assumed  model.   This  procedure  is  really  only  satisfactory  for 
parameterization  and  description  of  the  general  background  conditions 
and  generation  regimes,  since  as  shown  in  the  cat's  paw  study  strong 
intermittent  events  are  associated  with  large  profile  deflections  and  the 
development  of  unstable  shear  layers.   The  approach  to  a  logarithmic  pro- 
file with  increasing  averaging  time  is  merely  another  indication  that  long 
term  averages  can  obscure  the  physics  of  the  transfer  processes. 

The  longest  set  of  runs  was  then  analyzed  using  standard  time  series 
techniques.   The  data  is  presented  in  terms  of  spectral  growth  and  the 
related  changes  with  time  of  atmospheric  and  water  variables,  and  the 
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coherence-phase  relationship  between  horizontal  wind  speed  and  waves. 
Spectral  estimates  have  6  to  20  degrees  of  freedom  (dof ) •  the  80%  and 
90%  confidence  limits  of  the  *>  distribution  (at  $  =  ion  n  integer") 

are  shown  on  the  spectral  growth  and  spectral  density  estimate  diagrams. 
Garrett  (1970)  has  recently  demonstrated  the  independence  and  normal  dis- 
tribution of  the  real  and  imaginary  parts  of  the  Fourier  coefficients  im- 
plicit in  the  use  of  this  measure  for  mean  wave  fields.   The  confidence 
limits  are  applicable  only  to  each  individual  realization  of  the  spectrum, 
and  do  not  provide  estimates  of  variance  of  the  wavef ield  .for  times  longer 
than  the  overall  averaging  period;  e.g.,  the  30  second  spectra  are  strong- 
ly affected  by  the  presence  or  absence  of  swell  and  major  wave  groups. 
Coherence  magnitude  (  5  )  was  tested  on  the  assumption  that  the  variates 
were  random  with  a  bivariate  gaussian  distribution  (Amos  and  Koopmans, 
1963).   (As  demonstrated  in  the  cat's  paw  study  and  further  displayed 

herein,  this  is^a  poor  assumption  for  strongly  coupled  variables;  un- 

P 
fortunately  it  is  the  only  readily  available  method).   The  results  may  be 

generally  stated  (for  5  minute  time  segments)  as,  at  the  .05  level  of 

significance,  the  probability  that  the  hypothesis  If  =  0  is  incorrect 

for  a  measured   °  ^  0.5  is  97%.   (For  the  30  second  segments,  this 

probability  has  dropped  to  70%). 

Finally,  several  10  minute  segments  from  the  start,  middle,  and  end 

of  the  set  of  51  runs  were  analyzed  using  30  second  averaging  periods. 

The  data  is  presented  in  a  similar  fashion,  with  the  added  utilization 

of  analog  strip  chart  displays  of  'events'. 
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Results  and  Discussion 

A.  Boundary  Layer  Profiles.  With  the  assumption  that  the  averaged 
stress  (  x  )  in  the  near  surface  region  is  constant  with  height,  a  charac- 
teristic velocity  of  the  turbulent  boundary  layer  may  be  defined, 


This  is  termed  the  friction  velocity.   Using  similarity  arguments,  for  a 
neutrally  stable  atmosphere  the  only  characteristic  length  is  height. 
The  resulting  profile  equation  is 

0  <3r)   =  t£  •*•»  l~)  where  r 

is  Von  Karman's  constant  -  0.40  and  z  ,  the  roughness  length,  is  general- 
ly not  added  to  the  height  in  the  numerator  (see,  e.g.,  Lumley  and  Panofsky, 
1964). 

There  is  general  agreement  that  this  simple  model  gives  a  good 
description  of  the  mean  velocity  profile  in  conditions  of  neutral  stability 
over  both  land  and  water,  but  there  is  still  considerable  debate  on  the 
meaning  of  and  relationships  between  the  profile  parameters  (U.  and  z  ) 
and  other  variables  computed  from  the  model.   From  the  development  of  the 
model,  U^  and  z   are  independent  similarity  parameters  characterizing 
the  boundary  layer;  however,  Sadeh,  et  al,  (1971)  have  demonstrated  that 
once  the  roughness  elements  are  sufficiently  high  compared  to  boundary 

layer  thickness,  U.  and  z  become  local  variables. 
3  '   *      o 

Measurements  of  the  profile  parameters  over  waves  have  been  re- 
viewed and  discussed  in  detail  by  Wu  (1969a,  1969b,  1970)  and  Ruggles 
(1969,  1970).   Wu  states  that  the  wind  stress  (t  = p  0fc  )  is  supported  by 
form  drag  provided  by  wavelets  superimposed  on  the  larger  waves;  air  flow 
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separates  from  the  wavelets  whose  phase  speed  is  less  than  U^.   z  becomes 

a  measure  of  the  height  of  these  roughness  elements  (roughness  height 

=  30  z  )  and  is  related  to  friction  velocity  by  Charnock's  relationship, 

z  /  (  U*  / <\       )  =  constant.   This  is  manifested  by  an  increase  in  the 
o        J 

drag  coefficient  (Cdz  =  T '/ P^i         =   (  u* /Ut  )    )  with  U  up  to  a  'sat- 
uration' value,  U  =  15  m/sec.   Above  15  m/sec  the  wind  speed  is  larger 
than  the  average  wave  phase  speed  and  Cd  and  z  become  nearly  constant. 
Effectively,  the  Charnock  relation  implies  that  a  single  parameter  can 
be  used  to  describe  the  flow  and  thus  the  momentum  balance  of  waves  and 
atmosphere.   The  development  of  the  profile  model  on  the  other  hand  as- 
sumed that  U^  and  z   are  statistically  independent. 

Ruggles  showed  that  over  waves ,  the  wind  speed  is  more  appropriate- 
ly represented  as  the  sum  of  a  mean  profile  component  (supporting  atmos- 
pheric dynamics),  and  a  term  representing  the  correlation  between  U.  and 
z   (from  wind-wave  interaction).   The  'correlation'  wind  component  repre- 
sents that  portion  of  atmospheric  momentum  required  to  support  air-sea 
interaction  dynamics.   It  is  roughly  10%  to  20%  of  U    1  above  4  m/sec, 
and  negative  below  2  m/sec  (indicating  a  minimum  wind  speed  for  generation). 
Ruggles  also  found  marked  changes  in  z  and  Cd,  with  excessively  high  peaks 
and  excessively  large  variances,  at  around  U.  _   =  2,  A,  and  8  m/sec.   These 
anomalies  correspond  to  the  division  points  of  the  Beaufort  scale,  i.e., 
to  marked  changes  in  the  physical  appearance  of  the  sea  surface.   Similar 
anomalies  are  shown  by  Wu,  who  attributes  the  changes  to  the  differences 
of  flow  in  aerodynamically  smooth,  transition,  and  rough  regimes.   Hicks 
and  Dyer  Q.970)  found  a  simpler  relationship,  with  Cd  =  .0011  and  aero- 
dynamically smooth  flow  below  6  m/sec,  and  Cd  =  .0015  with  rough  flow  for 
U-n  =  10  m/sec.   They  did  not  report  any  anomalies. 
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Additional  evidence  of  anomalous  wind  speeds  was  given  in  the  work 
of  Groscup  (1971)  and  Mollo-Christensen  and  Groscup  (1971)  (based  on  the 
data  of  this  study).  They  further  demonstrated  a  bi-  and  tri-modal  dis- 
tribution of  z   in  the  vicinity  of  the  U. ~  =  8  m/sec  discontinuity,  and 

a  bimodal  distribution  of  U.  at  or  near  the  z   and  Cd  discontinuities  at 

*  o 

4  and  8  m/sec.   The  anomalies  are  indicative  of  enhanced  air-sea  inter- 
action at  specific  wind  speeds  and  of  changes  in  the  nature  of  wind-wave 
coupling.   The  large  variances  and  bimodal  distributions  are  evidence  of 
enhanced  nonlinear  coupling  at  critical  wind  speed/surface  "condition  values. 

Fig.  2  shows  the  distribution  (smoothed  histograms)  of  30  second 
average  wind  speed  and  of  atmospheric  logarithmic  profile  parameters  from 
this  study,  both  for  all  runs  and  for  the  subset  109-159  (those  used  in  the 
time  history  analysis).  Fig.  3  through  Fig.  6  show  selected  joint  histo- 
grams of  the  variables. 

IL  „  was  found  from  the  profile  model,  using  the  computed  values  of 

u\  and  z  .   The  drag  coefficient,  Cd,  was  determined  for  both  the  computed 

3 
10  m  and  measured  11.4  m  speeds.   U^   is  a  measure  of  turbulent  kinetic 

energy  production, 

31  Ki 


and    Cu    is  the  variance  of  velocity  in  the  profile,  o\,  "  l^''"u' 
The  histograms  show  that  a  large  percentage  of  the  runs  had  very 

small  roughness  length  (  J^\  (O  <  -  ~l  S   or  £0  <  1 .8  v  io  cw^  values  of 

_3 

z   <  10   cm  were  lumped  into  the  In  [z  ]   =  -10  category  so  that  sta- 
tistics were  not  computed).   Stress  values  are  in  rough  agreement  with 
the  long  term  averages  found  by  other  investigators;  the  mean  is  largest 
for  the  continuous  strong  winds  of  runs  109  -  159,  and  decreases  with 
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increasing  averaging  time.   The  values  of  Cd  are  comparable  for  the 
measured  and  computed  heights;  Ruggles  (1969)  computed  Cd   for  5  anemometer 
heights  and  10  m,  and  found  similar  results >  with  as  expected  a  slight 
decrease  with  increasing  z.   The  Cd  values  are  concentrated  in  the  normal 
range  for  runs  109  -  159,  but  when  all  runs  are  considered  together,  the 
30-second  and  10-minute  means  are  excessively  large  and  there  are  numerous 
exceptionally  high  drag  values.   These  major,  average  changes  in  the  drag 
effect  of  waves  on  the  wind  suggest  the  existence  of  more  than  one  wave 
generation  regime.    CJo   ,  Cdlo  5  and  U^  (and  thus  U^  and  t    )  all  have 
decreasing  mean  values  with  increasing  averaging  times,  indicating  that 
profile  distortion  and  resulting  turbulence  and  wave  generation  occur  in- 
frequently; the  short,  strong,  important  anomalies  are  obscured  by  long 
term  averaging. 

Fig.  3  show's  JH  (U    UJ,  JH(U1Q,  ln(z  )  ),  and  JH(U   ,  Cd   )  for  the 

30-second  average  profiles.   In  this  figure,  anomalies  of  U.  ln(z  ),  and 

*     o 

Cd  „  all  occur  at  roughly  2.5,  3.5,  6.4,  7.7,  10,  11,  and  12.2  m/sec  (they 
are  particularly  striking  at  2.5,  3.5,  and  7.7  m/sec,  i.e.,  near  those 
found  by  Ruggles).   The  anomalies  are  in  the  form  of  exceptionally  large 
variances  and  wide  departures  of  the  interval  means  from  any  postulated 

simple  fit  of  the  data.   The  values  are  in  general  quite  scattered  and 

2 
poorly  represented  by  means.   Also  shown  is  JH(U^  /g,  K)  where  K  =  30  z 

is  height  of  the  roughness  elements  (wavelets)  as  postulated  by  Wu.   It 

is  obvious  that  a  simple  relation  such  as  suggested  by  Charnok  and  found 

by  Wu  does  not  hold  here.   There  are  bimodal  distributions  and  anomalies 

in  both  directions.   A  single  parameter  can  not  characterize  the  profile, 

because  of  the  interplay  between  atmospheric  dynamics  and  the  wave  field. 

None  of  the  data  can  be  fit  with  a  linear  relationship;  in  particular, 
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Cd  x  I03 
(Wu's   formulae) 
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Ruggles's  suggested  U^  =  0.04Uin  intersects  the  minima  of  the  mean  of 


U*  on  Uio( 


Fig.  4  shows  JH(U   ,  Cd   )  computed  using  the  formulae  of  Wu  (1969b) 
for  'breeze1  and  'light  wind'  regimes, 

U-  '     ,i.(<t.lu.  ..<>-)   .  u<  5  *«  >  Q~  =  M    ■■>-   )  ^«<<Kie  '  "'«• 

* V '  l     l  .OISCUaA}  ) 

The  anomalies  at  2.5  -  3.5,  7.5,  and  12  m/sec  again  appear,  but  much  less 
strongly  than  in  direct  computation  from  the  profile  model;  in  addition, 
the  mean  value  of  Cd  computed  using  this  technique  is  excessively  small. 
Large  scatter  was  found  in  joint  histograms  of  other  relationships  sugges- 
ted by  Wu.   The  failures  of  these  postulated  simple  models  on  the  relative- 
ly short  (30  sec)  averaging  times  indicates  that  the  dynamics  of  wind  flow 
over  open  ocean  wave  fields  at  large  Reynolds  numbers  can  be  only  poorly 
described  in  terms  of  long  time  averages,  even  when  the  significance  of 
wavelets  is  considered.   Wu's  contention  that  the  small  surface  components 
support  the  majority  of  wind  stress  is  not  disproven  by  these  data;  the 
results  merely  indicate  that  the  dynamics  of  the  wind-wave  coupling  is 
considerably  more  complex  than  that  caused  by  just  form  drag  and  air  flow 
separation. 

Fig.  4  also  displays  30  second  JH(  r  ,  ln(z  )  ),  JH(Cdin,  ln(z  )  )  and 
JH(  Cr0  ^04.  /k?  )•   The  T  -ln(z  )  relationship  shows  a  distinct  bimodality, 
with  generally  increasing  stress  with  roughness  length  and  stress  maxima 
for  z  ~z    1  cw\  (ln(z  )  =  0) .   High  stress  occurs  predominantly  over  large 

'roughness  elements',  but  only  under  certain  conditions  as  yet  undetermined. 

2 
The  Cd   =   (U./U   )   vs  ln(z  )  plot  is  very  closely  fit  by  the  profile 

relationship,    Jk  (tj  *  t."b  ~  T^fprtt* 

(using  K  =  0,40;  Businger,  et  al,  1971,  have  recently  noted  that  *c   =  0.35 
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is  a  better  value  of  Von  Karman's  constant  for  neutral  stability  over  a 
land  surface);  this  is  merely  a  check  on  the  consistency  of  the  logarith- 
mic profile  assumption  and  on  the  continuity  of  data  from  different  time 
periods,  and  serves  to  emphasize  the  significance  of  the  anomalies  of 
the  other  plots  (i.e.,  they  appear  not  to  be  due  to  instrumental  or  analysis 
peculiarities,  and  the  0.94  correlation  is  a  sufficient  criterion  for 
assuming  a  logarithmic  distribution).   The  Ct*u  -  U*5  /*\         relationship 
shows  that  turbulent  production  increases  with  increasing  fluctuations,  but 
with  a  very  strong  bimodality,  indicative  of  the  presence  of  more  that  one 
production-generation  regime. 

Fig.  5  shows  JH(Uin,  ln(z  )  )  and  JH(Uin,  Cd10)  for  averaging  periods 
of  2  and  10  minutes.   The  basic  patterns  of  anomalies  at  specific  values 
of  Uin  persist  at  the  longer  averaging  times,  but  with  less  variance.   (The 
10-minute  averages  have  been  compared  by  Groscup  [1971]  to  those  of  Ruggles 
[1970]  [on  the  basis  of  loglO  (z  )]).   The  basic  patterns  remain,  irrespec- 
tive of  averaging  time,  but  the  decreasing  intensity  of  the  anomalies  shows 
that,  as  in  the  cat's  paw  studies,  the  strong  and  intermittent  events  are 
blurred  or  obscured  by  long  time  averaging.   Even  with  these  long  averages 
and  highly  logarithmic  profiles,  it  remains  clear  that  the  'mean'  is  a 
very  poor  characterization  of  physical  reality.   The  time  history  of  events 
and  the  regimes  and  instantaneous  characteristics  of  wind  and  wave  fields, 
are  highly  significant  in  determining  the  nature  of  the  boundary  layer  dy- 
namics. 

This  point  is  further  emphasized  in  Fig.  6,  which  shows  JH  of  U  _  with 
U^,  ln(z  ),  and  Cd   for  runs  109  -  159,  and  JH(U   ,  U^)  for  all  runs  for 
the  two  subclasses   |  r-<O>>o'se'03o-5ec   ^?  5°  ^/sec  anA  <r  10  o*fcec.      The  JH  of 
runs  109  -  159,  characterized  by  strong  and  increasing  wind  and  a  large  wave 
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field,  occupy  concentrated  segments  of  the  overall  diagrams.   Signifi- 
cantly, most  of  the  anomalies  are  missing,  and  as  previously  noted  Cd1 _ 

-3 
for  these  runs  was  0.9  x  10   lower  than  the  overall  average;  the  value 

for  these  runs  is  in  rough  agreement  with  that  found  by  most  investiga- 
tors.  This  set  of  runs  apparently  satisfies  the  criteria  for  a  single 
or  homogeneous  wave  generation  regime  (the  bimodalities  of  the  other  rela- 
tionships also  disappear  or  decrease  when  these  runs  are  plotted  separate- 
ly).  The  division  by  mean  velocity  derivative  magnitude  shows  that  the 
parameter  relationships  do  vary  with  time  history,  and  that  (as  in  the 
cat's  paw  study)  higher  wind  speeds  are  in  general  associated  with  stron- 
ger fluctuations. 

Fig.  7  shows  examples  of  the  results  of  the  JH  of  the  parameters  of 
the  logarithmic  underwater  profile  (runs  109  -  159).   30-second  mean 
current  speed  (at  10m)  was  approximately  10  cm/sec,  with  a  range  of  -3 
to  27  cm/sec.   The  diagrams  show  quite  consistent  relationships  among 
the  profile  values,  with  no  indication  of  the  anomalies  or  bimodalities 
found  in  the  atmosphere.   U^  values  are  negative,  as  the  wind  drift  is 
strongest  near  the  surface.   Joint  histograms  between  air  and  water  pro- 
file parameters  from  this  set  of  runs  failed  to  show  any  overall  sta- 
tistical relationships;  this  is  not  surprising  in  view  of  the  fact  that 
most  of  the  momentum  from  the  wind  goes  directly  into  the  wave  field 
(Stewart,  1967;  Dobson,  1971),  so  that  the  effect  on  the  surface  drift 
current  is  not  immediate  or  else  (see  section  B.2,  Figs.  28  and  29)  causes 
the  short  average  profiles  to  become  highly  nonlogarithmic. 

In  summary,  the  mean  velocity  profiles  in  both  air  and  water  are 
predominantly  logarithmic  during  conditions  of  neutral  stability.   Ano- 
malies in  the  atmospheric  profile  parameters  demonstrate  the  importance  of 
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time  history  and  the  existence  of  more  than  one  wind-wave  generation 
regime,  with  enhanced  activity  at  certain  critical  wind  speeds  under  as 
yet  unknown  circumstances.   Anomalies  stand  out  most  strongly  and  the 
air  boundary  layer  becomes  much  less  logarithmic  at  short  averaging 
times,  indicating  the  importance  of  short,  intermittent  events  in  the 
development  of  unstable  shear  layers  and  associated  instabilities.   The 
profile  can  not  be  well  represented  by  a  single  similarity  parameter,  and 
there  is  no  simple  relationship  between  wind  velocity  and  the  profile 
measures.   Attempted  use  of  such  relationships  for  characterizing  the 
wind  regime  and  predicting  wave  growth  can  lead  to  significant  errors. 
Consideration  must  be  given  to  the  time  history  of  generation  and  to  the 
importance  of  intermittency ,  i.e.,  to  the  enhanced  air-sea  coupling  and 
momentum  transfer  at  certain  wind  speeds  or  in  short  events. 
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B.   Time  History 

The  51,  10  minute/hour  runs  (consecutive  except  for  1  hour  gaps 
between  runs  20  and  21  and  42  and  43)  selected  for  detailed  analysis  were 
taken  from  21  to  23  October,  1970.   The  wind  direction  had  switched  from 
N  to  S  30  hours  prior  to  the  first  run,  as  a  large  high  pressure  system 
(>  1040mb  at  the  center)  passed  the  site  moving  NE.   Winds  increased  and 
weather  deteriorated  throughout  the  three  day  period,  then  winds  slacked 
and  weather  cleared  on  October  24.   Analysis  of  current  records  and  of  the 
tide  record  at  Wood's  Hole,  Massachusetts,  showed  that  the  surface  drift 
continued  NE,  into  Vineyard  Sound,  after  the  NE  flood  tide  on  the  morning 
of  October  21.   The  buoy  thus  faced  into  the  wind  and  waves  throughout 

,  this  set  of  runs.   This  was  confirmed  by  visual  observation. 
B.l.   5  and  10  Minute  Averages 

Figs.  8  and  9  depict  the  time  changes  of  the  10  minute 
average  horizontal  wind  speed  at  10m  (computed)  and  2.4m  (measured)  and 
the  wind  extrema  at  2.4m.   The  figures  also  show  current  speed  at  Z  = 
-1.52m,  and  stress  (  t  zp  U^,  )  computed  from  air  and  water  profiles.   Short 
lines  connect  the  averages  of  2,  5  minute  periods  from  each  run  and  give 
an  indication  of  the  variability  on  shorter  averaging  times  (the  large 
offsets  in  Runs  111  and  142  are  due  to  noise;  those  in  120  and  157  are 
real;  2.4m  velocity  extrema  from  runs  129-133  were  unavailable). 

There  are  strong  periodicities  in  all  of  the  variables.   Mean  wind 

•  shows  an  overall  increase  from  8-12m/sec;  winds  increased  slowly  through 
run  132,  then  increased  rapidly  to  lOm/sec,  stayed  near  that  level  for  10 
hours,  decreased  and  remained  moderate  for  a  similar  period,  and  then 
again  increased  rapidly  starting  around  run  152.   Superimposed  upon  this 
long  modulation  is  a  shorter  but  strong  periodicity  of  roughly  2-4  hours. 
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Velocity  shear  and  extrema  also  show  an  overall  increase  with  time,  with 
superimposed  shorter  periodicities.   The  water  current  has  a  strong  semi- 
diurnal tidal  component  with  alternating  times  between  maximum  flood  of 
5  and  7  hours;  the  current  maxima  around  run  139  is  strongly  suppressed. 
Some  of  the  wind  fluctuations  occur  with  the  tidal  frequency  and  its  2nd 
harmonic;  local  wind  maxima  frequently  occur  during  maximum  flood. 

The  effect  of  the  tidal  modulation  is  vividly  displayed  in  the 
water  boundary  layer  profile  as  an  excess  during  ebb,  when  the  wind  and 
tidal  components  are  opposing  and  relative  wind  to  water  speeds  a  maximum. 
The  short  2-4  hour  fluctuations  of  the  wave  field  are  also  apparent  in  the 
water  stress.   Air  stress  has  a  diurnal  periodicity  and  a  small  2-4  hour 
component  as  well  as  local  maxima  during  the  two  major  wind  increases;  it 
also  temporarily  increases  in  runs   141  -  146  when  the  winds  decrease  ra- 
pidly.   Mean  stresses  are  thus  largest  when  there  is  a  significant  distur- 
bance of  the  equilibrium  balance  between  current,  wind  and  wave  fields. 
The  stress  values  from  air  and  water  are  in  rough  agreement.   Differences 
may  be  partly  accounted  for  by  the  momentum  carried  by  propagating  swell 
and  wind  waves,  and  by  the  momentum  and  equilibrium  inbalances  caused 
by  the  tides. 

Mollo-Christensen  (personal  communication)  has  noted  that  except  for 
the  tidal  variations,  10  minute  mean  air  and  water  velocities  are  compar- 
able when  normalized  by  their  dynamic  viscosities,  i.e.,  uA  U^  =-,Mw^vj  . 
This  is  shown' in  Fig.  10  (from  Mollo-Christensen  and  Groscup  [1971]). 
This  relationship  can  be  accounted  for  by  assuming  that  the  main  shear  is 
in  the  viscous  sublayer,  so  that 
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If   ow  =■  da.     >  e«g«>  if  the  sublayer  thickness  in  both  media  is  controlled 
by  the  same  elements  (say  the  wavelets),  then  since  the  stresses  are  in 
rough  balance  the  boundary  layers  are  shear  coupled  and  related  by  the 
ratio  of  their  dynamic  viscosities,  ^Ja  r  -rr~    ^^     .   The  long  term  ave- 

-A*  3k 

rage  values  thus  again  support  the  assumption  that  the  major  air  sea  coup- 


tir 


ling  can  be  simply  described,  in  spite  of  the  known  intermittent  discon- 

•  20 

uities.   Fig.  10  also  shows  the  weak  relationship  between  the  shorter 
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Figure  10 
term  variations  of  air  and  water  stress. 

Fig.  11,  Fig.  12  and  Fig.  13  show  the  variability  with  time  of  the 

wave  field  parameters.   Fig.  14  shows  the  corresponding  wave  spectra,  and 

Fig.  15  the  development  of  selected  frequency  components  (cross  sections 

of  Fig.  14).   80%  and  90%  confidence  limits  at  <£   =  10  are  shown  by  ticks 

inside  the  d3  axis  of  Fig.  15.   E  is  mean  energy  density,   E  =  /^l2 

and  momentum  density  M  =  E/c  (c  =  wave  celerity)  computed  by  spectral 

summation.   C  (9M&<  )  and  M  ($MA<  )    are  respectively  the  celerity  of  the 

wave  of  maximum  spectral  density  and  the  momentum  computed  by  assuming 

that  the  energy  density  is  concentrated  at  that  one  frequency.   The 
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information  in  these  figures  may  be  summarized  as  follows. 

During  the  two-day  period  the  mean  energy  density  of  the  wave  field 

increased  by  nearly  an  order  of  magnitude,  from  2  x  10  to  17  x  10  ergs/ 

2 
cm  .   The  period  of  the  wave  of  maximum  spectral  density  (T)  shifted  to 

lower  frequencies,  from  4  seconds  at  run  109  to  >  8  seconds  at  the  end. 
These  changes  can  be  clearly  seen  in  Fig.  11,  Fig.  12  and  Fig.  14.   The 
changes  in  T  and  E  are  slow  and  comparatively  steady  for  the  first  23 
hours.   E  then  increased  rapidly  in  conjunction  with  the  rapid  wind  in- 
crease that  started  around  run  132.   Fig.  14  shows  that  the  spectral  peak 
turned  sharply  toward  lower  frequencies  during  this  velocity  increase; 
after  run  132  the  separation  of  the  spectrum  into  two  peaks  (swell  and 
locally  generated  waves)  becomes  more  distinct  and  the  peaks  continue 
to  separate  as  the  swell  peak  shifts  to  still  lower  frequencies.   The 
equilibrium  range  above  the  local  peaks  has  a  very  rough  f   slope 
(Phillips  [1958])  with  occasional  local  maxima  in  the  range  of  0.5  to  1 
hz.   The  distinct  separation  between  local  and  swell  peaks  was  wiped 
out  during  the  rapid  wind  increase  of  runs  148  and  149. 

All  of  the  figures  show  that  the  growth  of  the  wave  field,  in  terms 
of  both  mean  energy  density  and  of  individual  spectral  components,  is  not 
steady.   There  are  large,  significant  variations  of  E,  T,  normalized  par- 
ameters, and  of   $y,o   (-f  ),  the  spectral  density  at  selected  frequencies 
(Fig.  15).   E    (1  Hour)  (Fig.  13)  becomes  as  large  as  +  .1%  of  the  hourly 
average  E  during  later  runs.   The  hourly  changes,  both  growth  and  decay, 
are  largest  when  the  wave  field  is  most  energetic  and  during  major  wind 
speed  changes. 

Fig.  15  demonstrates  that  even  the  higher  frequency  components  (.6310 
and  .3550  hz) ,  which  are  always  above  the  spectral  peak  and  in  the  'equili- 
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briura1  or  saturation  range,  have  significant  (at  the  90%  confidence  level) 
increases  and  decreases  of  energy  content  on  these  one  hour  averaging 
periods.   The  .2515  hz  component,  which  was  at  the  spectral  peak  at  the 
start  of  the  runs,  appears  to  'overshoot'  in  runs  110-112.   Its  spectral 
density  then  falls  off  to  a  lower,  'equilibrium'  value  during  the  slow 
growth  period  that  extends  through  run  132.   After  that  run,  the  4  second 
waves  again  show  major  growth,  'overshoot'  slightly  at  the  same  spectral 
density  they  had  earlier  reached,  and  then  settle  into  a  second,  higher 
'equilibrium'  level.   In  runs  154  -  156  their  energy  content  increases  to 
still  higher  levels.   The  wave  field  can  apparently  support  more  than  one 
overshoot  and  saturation  range  under  different  wind  conditions. 

The  .1583  and  .1254  hz  components  both  grow  until  around  run  139, 
when  the  higher  of  the  two  frequencies  reaches  a  quasi-steady  state  level. 
The  most  rapid  growth  of  the  low  frequency  waves  occurs  in  the  rapid  wind 
increase  starting  at  run  132.   The  major  wind  decrease  of  runs  142  -  143 
is  not  accompanied  by  significant  decay;  the     TV|(f)  decrease,  and  E(l  hour) 


for  this  period,  are  no  larger  than  those  found  in  the  minor  fluctuations 
during  growth  periods.   The  fluctuations  of  the  individual  frequency 
components  is  not  in  phase;  some  portions  of  the  spectrum  grow,  others 
decay,  in  what  appears  to  be  a  random  fashion  except  during  the  major 
wind  increase  at  run  132,  which  is  followed  around  run  135  by  growth  of 
all  but  the  highest  frequency  components.   The  locally  generated  waves, 
.2515  and  .3550  hz,  show  similar  periodicities  to  those  of  the  local  wind 
field.   The  individual  swell  and  high  frequency  components  are  essential- 
independent  of  these  short  term  variations,  while  the  total  wave  field  - 
E,  T  and  M  -  does  show  the  effect  of  local  generation. 
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A  comparison  of  the  plot  of  U   /C  (Q,*^*  )  with  those  of  E  and  M 
shows  that  spectral  development  is  most  rapid  when  the  wind  is  moving  fas- 
ter than  the  waves  of  the  spectral  peak.   Some  growth  still  occurs  even 
when  the  waves  move  faster  than  the  wind  (where  Miles  theory  is  clearly 
inapplicable  due  to  the  absence  of  a  'critical  layer').   The  general 
range  of  U   /C  (  $„Af  )   is  quite  low  for  the  entire  set  of  runs.   These 
features  indicate  that  swell  was  propagating  into  the  area,  but  do  not  ac- 
count for  the  lack  of  local  decay  of  low  frequency  components  between  runs 
137  and  158  when  U   /C  (  qS MAif  )  was  continuously  less  than  1.   It  may  be 
hypothesized  that  energy  was  being  transferred  to  the  low  frequency  waves 
by  nonlinear  resonant  interactions  within  the  wave  field,  (e.g.,  Garrett 
[1970])  or  that  wind  generation  occurs  at  these  frequencies  during  brief 
but  intense  periods  when  U  was  larger  than  C  ( <lwv  ). 

The  plot  of  C  ( <|    )/U,  shows  a  strong  diurnal  periodicity,  and 
lags  the  similar  variation  in  I   .by  roughly   "H/2.   There  is  no  ap- 

3.1  L 

parent  relationship  between  this  parameter  and  U._/C  (^     ),  or  between 

10       max 

it  and  energy  density  growth  rates. 

The  normalized  wave  momentum  plot  indicates  that  only  60-80%  of  the 

total  wave  momentum  is  carried  by  waves  at  the  spectral  peak;  most  of  the 

remainder  is  carried  by  waves  at  higher  frequencies  (lower  phase  speeds). 

Fig.  16  shows  the  momentum  spectrum  (5  minute  averages,  6  dof)  from  several 

runs  (ordinates  have  been  normalized  to  the  spectral  maximum).   These 

plots  show  the  overall  increase  in  M  and  M    ,  the  shift  of  M    to  lower 

max  max 

frequencies,  the  bi- and  trimodality  associated  with  the  separation  of 
swell  and  local  wind  wave  peaks,  and  the  occasional  very  strong  concentra- 
tion of  momentum  in  a  single  swell  frequency.   The  momentum  spectra  have 
roughly  the  same  shape  as  the  power  spectra,  but  tend  to  emphasize  the 
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concentration  of  the  wave  momentum  in  one  or  more  narrow  frequency  bands. 

The  10%  momentum  levels  —  the  range  of  f  for  M(f)  >  . 1M    —  were 

°  max 

approximately  .15  to  .55  hz  through  run  131,  and  then  decreased  slowly  to 
roughly  .09  -  .40  hz.   Wave  momentum  was  indistinguishable  from  noise 
level  above  1.40  hz  near  the  start  of  the  runs,  and  this  limit  decreased, 
starting  at  run  132,  to  approximately  1.00  -  1.20  hz.   Most  of  the  mo- 
mentum spectra  had  several  discrete  peaks,  and  there  was  usually  a  consid- 
erable difference  in  the  minor  features  of  the  2,5  minute  segments  in 
each  run.   Current  spectra  at  -1.5m  (not  shown)  had  maxima  at  both  swell 
and  locally  generated  wave  frequencies,  and  were  generally  coherent  with 
both  of  these  wave  field  components. 

Coherence  between  horizontal  wind  speed  at  2.4m  and  waves  at  the 
spectral  peak  (  Oun  ((S  ^)  measured  on  the  5  minute  subintervals)  was 
y>      0.65  for  50%,  and  >,  0.50  for  97%  of  the  102,  5  minute  samples. 
Phase  (  ©o«  )  at  that  peak  was  -t  IT    for  82%  of  the  coherent  runs,  and 
>:^   2  radians  for  all  of  them.   As  noted  by  Seeshctz  (1968),  wind 
speeds  are  maximum  over  the  wave  troughs. 

Examples  of  YUv?  and  G>y   are  displayed  in  Fig.  17.   The  shift  of 
the  peak  to  lower  frequencies  with  time,  and  its  bifurcation  into  local 
and   swell  segments,  are  clearly  evident;  there  was  a  secondary  maximum 
of  l£0   at  the  local  wave  frequency  (roughly  .32  hz)  for  about  80%  of 
the  samples  after  run  131.   In  the  majority  of  these  cases  (note  run 
137-2)  there  was  a  decrease  in  Ua  between  the  two  peaks,  with  a  shift 
in   ®uy,  from  Tl   to  Wi  and  then  back  to  "P"  .   Winds  were  generally  in- 
coherent with  the  waves  (  Yo«  statistically  indistinguishable  from  0) 
above  the  local  wave  peak  on  this  averaging  time;  one  notable  exception 
is  segment  157-1;  this  run  contained  a  very  strong  generation  event  and 
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is  discussed  in  detail  in  section  B.2.   In  roughly  30%  of  the  runs,uo«  was 

fairly  high  at  the  D.C.  level;  in  all  of  these  cases,  Qv„—0   .   With  the 

exception  of  the  occasional  D.C.  peak,  <S0«  was  low  for  f  <   f  J 

7        "  *  max 

Wind  spectra  did  not  consistently  display  distinct  peaks  at  the  wave 

-.      4 
spectral  maxima  (Fig.  18,  80%  and  90%  confidence  limits  at  <£no  =  10  and 

Q>0u  =  10  are  shown  on  the  ^  axis  of  run  142)  .   They  were  fairly  white 

(  <  1  decade  dropoff  in  the  . 2  to  2  hz  band).   Nevertheless  there  was 


nearly  always  a  spectral  increase  at  or  around  f  £    ,  and  in  cases  of 
a  double-peaked  coherence  the  wind  spectra  showed  a  local  minimum  between 
the  peaks.   The  examples  show  4?™  (f)  and  <P  o\j  (f)  for  the  same  run  seg- 
ments as  in  Fig.  17  (they  are  rough,  6  dof,  and  the  jaggedness  is  probab- 
ly not  significant) . 

From  these  facts  and  figures  the  following  inferences  may  be  drawn. 
Wave  generation  and  spectral  development  is  strongly  non-stationary.   There 
are  significant,  rapid  changes  in  spectral  growth  and  these  changes  are 
poorly  represented  by  long  term  means;  no  single  parameter  or  simple  com- 
bination of  parameters  can  describe  the  process.   Air  and  water  boundary 
layers  (both  velocities  and  stress)  are  viscously  shear  coupled.   Modu- 
lations at  the  tidal  frequency  and  its  harmonics  (  or  of  the  period 
ATide/  A  t)  appear  in  the  wind  and  wave  fields;  changes  in  the  tidal  flow 
disrupt  the  equilibrium  between  wind  and  wave  fields  and  cause  modifica- 
tion of  momentum  advection  patterns. 

Even  though  the  high  frequency  portions  of  the  wave  field  are  near 
saturation  and  account  for  insignificant  portions  of  the  spectral  density 
and  momentum,   they  show  significant , rapid  (hourly)  variations  and  can 
account  for  spectral  growth  when  the  U   /C($   ^relationship  would  normally 
predict  major  continuous  damping.   Measured  growth  rates  (long  term 
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averages)  of  these  high  frequency  components  are  inherently  incapable  of 
displaying  major  changes,  since  on  the  mean  these  waves  are  near  satura- 
tion, and  also  because  the  very  high  frequencies  axe  continuously  being 
generated  by  the  larger  waves  (Longuet-Higgins, [1963] ;Crapper  [1970]). 
The  importance  of  these  high  frequency  waves  lies  in  their  ability  to 
couple  the  wind  with  larger  waves  and  to  assist  the  growth  of  low  fre- 
quency instabilities  by  nonlinear  coupling  between  scales.   It  will  be 
shown  in  the  next  section  that  these  components  have  extremely  rapid 
growth  rates  in  intermittent  generation  'events',  but  again  these  impor- 
tant physical  momentum  transfer  features  are  obscured  by  the  averaging 
process. 

The  influence  of  waves  on  the  turbulent  atmospheric  boundary  layer 
is  strong  at  both  the  swell  and  the  locally  generated  spectral  peaks. 
Furthermore,  the  2-4  hour  periodicity  of  wind  and  both  the  locally  gen- 
erated and  overall  wave  fields   are  strongly  related.   Changes  in  measures 


of  the  wind  and  wave  fields  (e.g.,  E  (1  hour))  are  largest  when  the  mag- 
nitude of  the  measures  is  largest,  similar  in  nature  to  the  pattern  of 
higher  frequency  components  during  strong  gusts  that  was  noted  in  the 
cat's  paw  study. 

The  major  conclusion  from  this  section  is  that  study  of  the  time  his- 
tory of  generation  and  spectral  development  is  essential  for  an  understan- 
ding of  the  physics  involved.  Wind,  wave  and  current  fields  show  evidence 
of  strong  interrelationships  that  are  poorly  represented  by  parameteriza- 


tion of  the  mean.   The  rapid  variations  in  E,  E  (1  hour),  U   and  £   all 
indicate  that  generation  is  predominantly  an  intermittent  process,  char- 
acterized by  rapid  transfer  of  momentum  and  followed  by  weaker  couplings 
and  nonlinear  interchange  between  spectral  components. 
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B.2.   30  Second  Averages. 

Several  of  the  10  minute  runs  were  chosen  for  de- 
tailed analysis  to  help  determine  the  nature  and  scales  of  intermittency 
and  the  resultant  momentum  transfer.   Analysis  and  data  display  were  simi- 
lar to  those  used  for  the  10  minute  runs. 

|  Fig.  1Q  anH  Fig.  ?0  depict  the  changes  with  time  of  30  s econd 
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Figure  22 
averages  of  variables  in  the  air  and  water  for  runs  110,  127,  157  and  159. 


Fig.  21  shows  E  (30  sec)  and  C  (fa  ) /U .  for  the  same  runs.   The  manor 

K  max   * 

features  of  the  10  minute  runs  also  stand  out  in  these  displays  of  short 
term  change.   There  is  a  general  increase  in  U,  E,  M,  T,  wave  height  and 
spectral  density  and  momentum  extrema.   Stress  computed  from  the  air  and 
water  profiles  are  in  rough  agreement,  although  they  show  different  perio- 
dicities.  The  short  term  changes  are  largest  and  most  rapid  when  the  mag- 
nitudes are  largest;  changes  in  U  are  paralleled  by  increases  in  E,  though 
not  as  clearly  as  in  the  longer  averages  because  of  the  influence  of  wave 
groups  (these  had  a  roughly  40  second  periodicity  for  later  runs,  and  also 

significantly  affected  the  T..  measurements).   U    is  again  normally  over 

0         J  w  max     °  J 

the  troughs,  but  with  much  scatter;  this  is  illustrated  in  Fig.  22,  JH 

(U„  .,  V)    from  2  short  segments,  and  is  more  clearly  shown  in  the  coherence 


. 
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diagrams  (Figs.  29  and  30). 

In  these  runs,  the  waves  of  the  spectral  peak  carry  roughly  80%  of 
the  total  momentum  density  (again  implying  that  a  significant  portion  is 
carried  by  waves  of  higher  frequency,  or  slower  phase  speed) ,  while  the 
waves  of  the  momentum  peak  (when  different  from  the  spectral  peak)  over- 
account  for  the  total  momentum.   U,_/C(<$>     )  is  in  the  range  of  1  -  2, 

10      max 

and  is  in  fact  largest  during  the  earliest  run,  due  to  the  absence  of 
swell.   This  parameter  shows  no  significant  peaks  at  the  points  of  wind 
and  E  maxima;  this  may  be  interpreted  to  mean  that  the  wind  is  strongest 
over  the  largest  wave  components,  indication  of  selective  coupling. 
C  (  $     )/Uj.  shows  no  distinct  relation  to  spectral  growth  or  decay; 

HlctX 

E  (30  sec)  has  comparable  variations  in  runs  110  and  127,  and  the  para- 
meter is  similar  in  magnitude  in  runs  110  and  159  when  the  energy  growth 
patterns  are  radically  different. 

For  these  short  averaging  times  E(30  sec)  is  1%  to  10%  of  the  mean 
energy  density  for  the  corresponding  10  minute  periods.   The  rapid  changes 
in  spectral  shape  and  content  are  well  shown  in  Fig.  23,  from  run  157. 
Cross-sections  of  this  diagram  are  shown  in  Fig.  24  (n.b.  80%  and  90% 
confidence  limits  at  10  on  <f  „Jf )  axis) .   It  is  immediately  apparent  that 
the  lower  frequencies  (referring  here  both  to  swell  [.1254  hz]  and  to 
local  wind  generated  waves  [.2515  hz])  are  relatively  constant,  with 
variations  primarily  attributable  to  propagating  wave  groups  (this  was 
checked  by  examination  of  analog  strip  chart  displays  of  the  filtered 
signals).   The  change  in  the  high  frequency  components,  however,  is 
very  dramatic,  with  30-second  variations  of  up  to  2  orders  of  magnitude 
(3  orders  of  magnitude  for  the  4.46  hz  waves,  Fig.  23).   The  change  of 
shape  of  the  'equilibrium  range'  in  Fig.  23  also  demonstrates  this 
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variation;  there  is  a  repeated  pattern  of  large  excesses  in  specific  high 
frequencies  followed  by  a  filling  out  of  the  spectrum.   The  subjective  im- 
pression of  high  frequency  excesses  climbing  the  slope  of  the  spectrum 
and  transferring  their  energy  to  lower  frequency  components  is  hard  to  avoid, 
Fig.  24  also  shows  that  the  excesses  do  not  always  occur  simultaneously  in 
all  frequency  bands  and  that  on  the  30  second  averaging  period  the  high  fre- 
quency growth  does  not  always  correspond  to  the  presence  of  low  frequency 
peaks.   Further,  intermittent  high  growth  rate  events  are  of  short  duration, 
typically  lasting  between  30  and  60  seconds;  i.e.,  similar  to  the  period 
of  the  major  wave  groups. 

Changes  in  the  wave  field  and  the  wind-wave  relationships  through 
an  'event'  are  illustrated  for  the  burst  of  segment  8,  run  157,  in  Figs. 
25  and  26.   Fig.  25  shows  the  coherence  of  U  and  r\     from  before,  during, 
and  after  the  event  (segments  6-10);  Fig.  26  shows  the  momentum  spectra 
from  the  same  period.   Before  the  event  Ou^  is  high  and  ©o,,  = -TT  f or  the 
low  frequency  components;  Xo*     is  insignificant  for  most  higher  frequencies. 
During  the  event  Yo«j  increases  dramatically  in  the' 1. 25-3. 25  hz  range.  ©0*j 
shows  a  continuous  'slip'  through  this  range;  this  is  probably  accounted 
for  by  the  averaging  time  being  much  larger  than  the  high  frequency  insta- 
bilities so  that  both  generation  and  detritus  is  represented,  and  by  the 
deceleration  of  the  gust  parcel  during  generation.   Immediately  after  the 
event  Yo«  has  decreased  for  frequencies  >1  hz;it  is  still  lower  in  segment 
10,  except  for  the  local  peak  at  1.25  hz  where  the  waves  lead  the  wind  (an 
opposite  phase  relationship  from  that  of  segment  9).   Before  the  event  the 
wave  momentum  is  strongly  concentrated  in  the  local  waves  (0.32  hz)  and 
swell  (0,15  hz).   During  the  burst  the  high  frequency  components,  out  to 
over  3  hz,  carry  a  much  larger  portion  of  the  momentum  than  at  any  other 
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time.   Immediately  after  the  event  the  very  high  frequency  components  have 
disappeared,  and  the  0.45  hz  component  has  been  strongly  augmented.  This 
peak  shifts  to  a  still  lower  frequency  (0. 40  hz)  in  the  following  segment, 
and  eventually  nearly  all  the  momentum  is  again  concentrated  in  the  local 
waves  and  swell. 

Fig.  20  and  Fig.  21  show  that  the  event,  on  the  30  second  averaging 
basis,  was  notable  by  anomalously  high  values  of  U   >  wind  shear,,  wave  height, 
current  maximum,  E,  E(30  sec),  and  M.   T  increased  slightly  during  the 
event,  and  significantly  60  seconds  later.  ~C     .   is  not  excessively  large, 

cxXjT 

but  the  profile  for  segment  8  had  a  low  correlation  with  the  model 
( P  -   0.91)  so  that  the  computed  value  of  ~C    is  unrealistic.   The  large 
departure  of  profiles  from  the  logarithmic  mean,  and  the  development  of 
unstable  shear  layers,  are  characteristic  of  such  intermittent  events  (as 
noted  in  the  cat's  paw  study  and  further  shown  in  Figs.  28  and  29). 

Fig.  27  depicts  the  detailed  (analog  strip  chart  output)  structure 
of  such  a  burst  type  generation  event  from  run  116.   Isotachs  (from  profile 
measurements),  wind  velocity  fluctuations  at  two  heights  (high  pass  filter- 
ed at  0.05  hz) ,  wave  height,  and  a  measure  of  the  energy  content  in  low  and 
high  frequency  bands  of  wind  and  waves  are  depicted.   Scales  are  relative. 
The  major  burst  occurs  after  a  rapid  wind  increase  and  is  associated  with 
the  coupled  occurance  of  low  frequency  wave  groups  in  air  and  water.   The 
generation  of  atmospheric  turbulence  (very  high  frequency  wind  fluctuations) 
occurs  predominantly  over  the  wave  peaks  in  air  and  water,  and  these  peaks 
occur  simultaneously.   Both  low  frequency  and  high  frequency  wind  and  wave 
energy  are  maximum  during  the  burst,  and  some  intermittent  activity  con- 
tinues after  the  major  gust  had  passed.   The  low  frequency  wave  components 
are  slightly  damped  after  the  gust  leading  edge,  when  the  burst  is  strong- 
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est.   This  figure  clearly  demonstrates  a  very  strong  coupling  between 
media  and  between  widely  separated  frequency  components. 

Figs.  28  and  29  are  similar  diagrams  from  the  first  and  third  E  peaks 
(see  Fig.  20)  of  run  159  (the  *|  scale  is  comparable  in  Figs.  27  and  29; 
all  others  are  relative).   These  figures  again  show  that  intermittent,  high 
growth  rate  wave  and  turbulence  generation  occurs  during  the  coupled 
occurrance  of  large,  low  frequency  groups  in  wind  and  surface  wave  fields. 
The  event  is  preceeded  by  profile  distortions,  a  wind  decrease  then  rapid 
increase.   In  Fig.  28  the  low  frequency  components  are  again  seen  to  be 
slightly  damped  during  the  intense  high  frequency  instabilities.   Large, 
low  frequency  components  of  the  wind  and  wave  fields  that  occur  indepen- 
dently do  not  lead  to  this  intense  coupling  and  momentum  transfer  to  high 
frequency  components;   it  is  the  strong  nonlinear  coupling  between  the 
media  which  results  in  the  generation  of  high  frequency  instabilities  and 
enhanced  coupling  of  scales. 

The  strong  shear  layers  in  the  atmospheric  profiles  are  of  very 
short  duration.   They  are  almost  invariable  associated  with  large  wave 
groups  and  high  frequency  wave  energy  spikes,  and  then  quickly  revert  to 
logarithmic  form  as  the  mean  wind  energy  is  drained  and  the  atmospheric 
dynamics  recover  from  the  local  instabilities.   The  current  profiles  show 
the  effects  of  the  momentum  of  the  large  wave  groups;  they  have  little 
direct  relationship  to  the  instantaneous  wind  profiles. 

The  coherence  and  phase  diagrams  for  the  last  event  of  run  159  are 
shown  in  Figs.  30  and  31.   Most  of  the  large  energy  peak  for  this  time 
period  was  due  to  the  presence  of  an  exceedingly  large  wave  group,  and 
the  high  frequency  generation  was  less  intense  than  in  the  similar  set  from 
run  157.   Nevertheless,  the  sequence  of  events  associated  with  the  instabil- 
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ity  are  again  apparent.   There  is  the  same  pattern  of  high  Yon  at  low 
frequencies  and  momentum  carried  in  long  waves,  followed  by  enhanced  high 
frequency  coupling  and  temporary  momentum  excesses  in  higher  frequencies. 

These  figures  have  demonstrated  the  short  duration,  strong  coupling, 
rapid  growth  rate,  highly  nonlinear  and  intermittent  nature  of  wave  gen- 
eration and  spectral  development.   Even  the  30  second  averages  blur  the 
intensity  of  the  coupling  by  overlapping  between  generation  events  and 
detritus  or  propagation.   The  microstructure  is  of  critical  importance, 
both  because  the  majority  of  momentum  transfer  occurs  at  high  frequencies 
and  because  the  high  frequency  instabilities  and  waves  serve  as  catalysts 
in  the  coupling  and  low  frequency  interactions  between  scales  and  media. 
Long  sequences  of  measurements  and  long  averages  can  predict  the  background 
conditions  under  which  generation  occurs,  and  can  show  the  effects  of  very 
low  frequency  coupling,  but  are  incapable  of  depicting  the  essential 
physics  of  wind-wave  generation. 

Summary  of  Results  and  Conclusions 

The  cat's  paw  study  showed  the, importance  of  'bursts'  in  initial 
generation  of  high  frequency  waves  and  turbulence  and  indicated  the  signif- 
icance of  intermittency  and  strong  nonlinear  coupling.   The  data  from 
this  study,  where  a  much  larger  dynamic  range  of  scales  is  involved, 
support  those  conclusions  and  the  following  statements. 

1.   Wave  generation  is  intermittent  and  involves  strong  non- 
linear coupling  among  a  wide  range  of  frequency  components  in  the  wind  and 
wave  fields.   The  high  frequency  instabilities  are  critical  to  the  momen- 
tum transfer  process,  and  much  of  the  air  to  sea  momentum  flux  goes  directly 
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into  high  frequency  waves.   The  important  scales  are  those  of  the  major 
waves  and  wave  groups,  and  those  of  the  high  frequency  wavelets  coupled  to 
them.   These  facts  are  substantiated  by  the  large  and  rapid  fluctuations 
in  E,  M,  and  E   for  both  1  hour  and  30  second  averaging  periods  (Figs.  11-13, 
19,  20,  and  23),  by  the  simultaneous  occurrance  of  large  magnitude,  coher- 
ent fluctuations  of  wave  height,  velocity,  and  energy  content  in  widely 
separated  frequency  bands  in  air  and  water  (Figs.  25,  27-30),  and  by  the 
resultant  excesses  in  spectral  density  and  momentum  in  high  frequency 
waves  (Figs.  23,  24,  26,  31). 

2.  The  high  frequency  wavelets  act  as  drag  elements  to  couple 
the  mean  wind  to  the  major  wave  field;   they  also  aid  the  growth  of  lower 
frequency  instabilities  by  nonlinear  stress-strain  coupling.   This  was 
demonstrated  in  tiie  cat's  paw  data,  and  is  shown  herein  by  the  variations 

in  coherence  patterns  between  wind  and  waves  (Figs.  25,  30)  and  the  coupled 
occurrance  of  bursts  in  widely  separated  frequency  bands  (Figs.  27-29).  The 
results  of  Wu,  who  makes  a  similar  contention,  are  not  supported  by  the 
short  term  profile  averages,  primarily  because  of  the  complicated  dynamics 
of  air  sea  interaction. 

3.  The  wave  field  is  modulated  by  surface  currents;  this 
perturbation  in  turn  causes  an  imbalance  in  the  equilibrium  between  wind 
and  wave  fields.   The  wind  and  surface  current  "fields  are  viscously  shear 
coupled  and  their  profiles  yield  comparable  mean  stress  values  (Figs.  8-10). 

4.  There  appear  to  be  a  variety  of  wave  generation  regimes. 
These  regimes  can  be  only  poorly  represented  by  mean  parameters  or  by 
simple  fits  to  the  mean  profile  conditions.   The  different  regimes  are 
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characterized  by  bimodal  distributions  of  profile  parameters  (Figs.  2-5) 
that  are  not  present  when  a  'homogeneous'  set  of  runs  is  considered  separ- 
ately (Fig.  6).   Profile  parameters  for  a  non-homogeneous  set  of  runs  have 
larger  variances  and  a  much  more  complicated  pattern  of  values  than  do 
those  of  a  homogeneous  set,  and  there  are  no  simple  formulas  which  can 
adequately  describe  the  relationships  between  them. 

5.  Study  of  time  history  is  important  for  an  understanding 
of  spectral  development.   The  time  history  is  poorly  represented  by  mean 
parameters  or  long  averages.   The  approach  of  the  profile  to  a  logarithmic 
form  with  increasing  averaging  time  obscures  the  importance  of  the  insta- 
bilities associated  with  short,  intermittent  perturbations  and  development 
of  shear  layers  (Figs.  28,  29).   The  long  average  spectra  (Figs.  14-16) 
can  depict  the  general  growth  of  the  wave  field  but  fail  to  reveal  the 
significant  flux  to  high  frequencies  during  intermittent  generation  (Figs. 
23,  24).   The  rough  quantitative  agreement  between  stresses  and  viscously 
normalized  velocities  in  air  and  water  boundary  layers  is  strongly  modulated 
by  tidal  forces. 

6.  Boundary  layer  profiles  in  air  and  water  are  essentially 
logarithmic  for  averaging  times  >  30  seconds  and  for  neutral  stability. 
The  profile  parameters  U^  and  z0  can  not  be  expressed  in  terms  of  a  single 
parameter.   The  criterion  p  >*   0.94  for  a  fit  of  .measured  values  to  the  model 
is  satisfactory  (Fig.  4).  There  do  exist  anomalous  or  'critical'  wind 
speeds,  in  the  sense  that  around  these  speeds  the  variability  in  profile 
structure  is  increased  and  the  resulting  profile  parameters  have  anomalous 
means,  variances,  and  amplitude  distributions  (Figs.  3-5).   The  mechanics 
responsible  for  these  anomalies  have  not  yet  been  demonstrated. 
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7.  The  growth  of  the  major  portion  of  the  spectrum  is  hy- 
pothesized to  be  due  both  to  the  wind-wave  coupling  enhanced  by  the  high 
frequency  wavelets  and  to  nonlinear  transfer  of  energy  within  the  wave 
field.   The  shift  from  high  to  low  frequency  content  is  best  illustrated 
by  the  momentum  spectral  changes  in  Figs.  26  and  31.   Spectral  growth  can 

occur  even  at  U/C(<f>    )  <    1  (Figs.  11,  12)  and  this  fact  is  unexplained 

ni  9.x 

by  either  propagation  or  by  available  theories. 

8.  Both  the  overshoot  phenomenon  and  the  equilibrium  range 
are  meaningful  only  in  terms  of  long  averages.  The  wave  field  can  support 

a  range  of  equilibrium  levels  under  different  wind  conditions  (Fig  15,  par- 
ticularly the  .2515  hz  component),  and  the  wave  field  is  never  in  a  steady 
saturation  state  (Figs.  15,  23,  24). 
* 
On  the  basis  of  these  results  it  is  concluded  that  a  large  portion 
of  wave  generation  occurs  intermittently  and  is  due  to  strong  coupling 
between  a  large  range  of  scales  of  motion  and  between  wind  and  waves.   Local 
instabilities  and  direct  nonlinear  interactions  between  waves  and  turbu- 
lence are   crucial  factors  in  generation,  and  most  of  the  momentum  flux 
from  these  'bursts'  goes  directly  into  high  frequency  waves. 
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IV.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  questions  posed  in  the  introduction  may  be  answered  as  follows: 

1.  The  important  time  scales  in  the  transfer  of  momentum  from 
wind  to  waves  are  those  of  and  within  the  intermittent  turbulent  bursts. 
These  in  turn  are  determined  by  nonlinear  coupling  to  the  period  of  the 
major  waves  (localized  instabilities  or  subbursts)  and  the  wave  groups 
(the  gusts  or  bursts  themselves).   The  transfer  is  rapid,  produces  high 
growth  rates  of  high  frequency  waves,  and  quickly  drains  the  available 
energy  from  the  mean  wind  field.   Nonlinear  interactions,  and  transfer  of 
energy  within  the  wave  field  by  breaking  and  dissipation,  operate  on 
similar  time  scales.   Transfer  is  effectively  completed  within  two  or 
three  periods  of '.the  waves  or  wave  groups.   The  time  scale  of  equilibrium 
imbalance  and  restoration  between  the  mean  wind,  wave,  and  current  fields 
is  on  the  order  of  2  to  4  hours. 

2.  The  fluctuations  involved  in  wind  wave  generation  are 
intermittent  and  strongly  nonlinearly  coupled,  involving  interactions 
between  a  wide  range  of  scales.   The  short  instability  fluctuations,  and 
the  resultant  energy  density  growth  rates  of  high  frequency  waves,  are 
roughly  three  orders  of  magnitude  larger  than  those  observed  on  the  mean. 
The  generation  of  waves  and  the  generation  of  turbulence  are  inseparable. 

3.  There  are  several  types  of  coupling  between  media.  These 
include  the  strong  interactions  during  instability  fluctuations;  the  form 
drag  of  the  wavelets  on  the  mean  wind  flow  and  the  resultant  coupling 
between  low  frequency,  large  amplitude  wave  groups;  the  interaction  of 
wave  induced  Reynolds  stresses  at  swell  and  locally  generated  wave  fre- 
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quencies  with  turbulent  Reynolds  stresses;   an  undetermined  anomalous 
interplay  between  wind  and  waves  at  discrete  critical  wind  speeds;  and 
viscous  shear  coupling  between  mean  wind  arid  surface  drift.   Many  of  these 
interactions  appear  to  depend  strongly  on  the  presence  and  position  of 
high  frequency  fluctuations  (the  microscales) . 

Coupling  between  scales  involves  both  the  weakly  nonlinear  resonant 
wave  interactions,  and  the  strongly  nonlinear  stress-strain  relationships 
involved  in  the  interactions  between  Reynolds  stress  and  shear.  These 
latter  operate  both  within  and  between  media.  High  frequency  instabilities 
aid  the  growth  of  lower  frequencies;  fluctuation  magnitudes  are  largest 
when  the  range  of  scales  involved  in  the  fluxes  are  largest. 

4.  Generation  'events'  may  be  distinguished  by  the  short 
time  average  approach  to,  and  deviations  from,  lognormal  distribution  of 
velocity  fluctuations;   by  strong  coherence  between  high  frequency  waves 
and  turbulent  Reynolds  stress  or  horizontal  velocity;   by  temporary 
excesses  of  energy  and  momentum  in  high  frequency  (>  1  hz)  surface  waves; 
and  by  the  enhancement  in  magnitude  and  frequency  of  stress  and  buoyancy 
fluctuations  caused  by  the  strong  nonlinear  couplings.   These  features  can 
not  as  yet  be  parameterized  in  terms  of  the  mean  or  of  profile  parameters, 

but  can  be  used  for  conditional  sampling  in  a  real  time  data  collection/ 

2 
analysis  system.   The  criterion  T    <=  100  dynes/cm  ,  when  smoothed  over 

a  time  period.^  0.5  seconds,  is  a  good  first  approximation  for  a  switch 

to  determine  the  presence  of  a  burst  within  a  cat's  paw. 

5.  There  are  no  easily  measurable  similarity  parameters 
which  can  describe  the  dynamics  of  air-sea  coupling.  The  apparent  effect 
of  fluctuations  depends  strongly  on  the  averaging  scale  over  which  they 
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are  evaluated,  i.e.,  the  effect  is  different  for  different  physical 
scales.   There  are  more  than  one  distinct  regimes  of  generation  and  there 
are  large  anomalies  in  the  relationships  between  wind  field  parameters. 
There  is  a  gross,  qualitative  similarity  in  the  appearance  of  weak  and 
strong  fluctuations,  and  of  large  bursts  and  subbursts,  but  the  intense 
nonlinear  interactions  and  resulting  enhancement  of  the  magnitude  of 
fluctuations  and  fluxes  preclude  any  strict  similarity  representation  of 
the  physics  of  momentum  transfer. 

The  mechanics  of  the  development  of  a  wind-wave  field  can  be  heuris- 
tically  described  in  terms  of  these  features: 

Initial  generation  is  caused  by  the  forcing  of  random  turbulent 
pressure  fluctuations  (the  Phillips  resonance  mechanism) ;  the  critical 
time  scales  involved  in  this  coupling  are  those  of  the  waves  generated. 
In  the  presence  of  a  strong  thermal  stratification  the  generation  of  waves 
is  delayed  until  the  boundary  layer  is  destabilized.  The  work  done  on  the 
waves  creates  reaction  stresses  in  the  air  which  generate  small  scale 
turbulence.   The  small  scale  fluctuations  do  work  on  the  larger  scale 
shear  in  air  and  water,  which  aids  the  development  of  larger  scale 
instabilities  and  stronger  fluctuations.   The  energy  from  the  mean  flow 
is  rapidly  depleted  through  transfer  to  the  water  surface  and  dissipation 
to  heat.   The  subbursts  occur  predominantly  within  larger  scale  anomalies 
or  gusts  that  are  initially  triggered  by  atmospheric  dynamics. 

After  an  isolated  gust  has  passed,  buoyancy  forces  rapidly  reestablish 
boundary  layer  stability.   The  packet  of  free  surface  waves  generated  by  the 
gust  and  bursts  appears  as  roughness  elements  to  the  wind.  In  the  presence 
of  increasing  wind  or  recurrent  gusts  the  disruption  of  flow  caused  by  the 
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roughness  elements  acts  to  enhance  the  probability  of  occurrence  of 
instabilities   over  the  packet.   The  scales  of  interaction  between  media 
have  increased  to  include  not  only  the  bursts  and  subbursts  but  also  the 
group  or  patch  of  roughness  elements  and  the  major  waves  within  it.   In 
successive  interactions  the  scales  are  similarly  increased  and  the 
energy  in  the  initial  waves  and  groups  which  enhanced  the  coupling  is 
augmented.    The  scales  are  yet  further  enlarged  by  nonlinear  resonant  wave 
interactions  within  and  among  the  wave  groups.   The  very  high  frequency  wave- 
lets, which  received  much  of  the  momentum  during  the  localized  high  fre- 
quency instabilities,  either  break  or  otherwise  transfer  their  energy 
to  other  components,  or  are  'organized'  by  the  convergences  and  divergences 
of  orbital  velocity  in  the  underlying  longer  waves.   The  concentration  of 
high  frequency  components  on  the  crests  and  leading  faces  localizes  the 
drag  elements  in  a  favorable  position,  since  stresses  exerted  upon  them 
are  in  the  same  direction  as  the  orbital  velocity  and  celerity  of  the 
longer  waves.   Under  the  influence  of  the  drag  elements  the  mean  wind 
decelerates  over  the  crests,  accelerates  over  the  troughs. 

As  the  lower  frequency  waves  increase  in  size  their  effect  on  the 
atmospheric  boundary  layer  increases.   Fluctuations  in  the  wind  field 
become  coherent  with  these  components,  and  wave  induced  or  reaction 
Reynolds  stresses  develop  which  both  enhance  the  pressure  feedback  from 
the  mean  flow  to  themselves  (Miles 's  mechanism)  and  interact  with  atmos- 
pheric shear  and  turbulence.   The  surf ace- waves  are  generated  and  propagate 
in  groups,  and  enhanced  perturbations  of  the  boundary  layer  on  the  scale 
of  the  groups  causes  similar  forcing  and  feedback  on  that  scale. 

The  selectively  increased  forcing  of  the  wind  field  by  the  longer 
waves  and  wave  groups,  and  by  the  roughness  elements  associated  with 
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the  major  crests,  increases  the  probability  that  shear  layer  instabilities 
will  be  triggered  over  them.   The  gusts  are  thus  selectively  coupled  to 
these  areas.   During  these  instabilities,  as  in  a  cat's  paw  over  a  calm 
surface,  the  majority  of  momentum  transfer  occurs  during  the  high  frequency 
bursts  and  subbursts,  and  the  momentum  flux  goes  directly  into  the  high 
frequency  wave  components.   Energy  growth  rates  are  locally  three  orders  of 
magnitude  larger  than  the  mean,  and  the  available  atmospheric  energy  is 
rapidly  depleted.   The  high  frequency  waves  again  break  or  otherwise 
transfer  their  energy  to  the  longer  waves  and  swell.   The  process  occurs 
intermittently  and  can  cause  spectral  growth  even  when  the  largest 
waves  are  moving  faster  than  the  mean  wind. 

The  generation  and  spectral  development  process  is  self-enhancing 
and  involves  strong  coupling  between  scales  of  motion  in  both  wind  and 
wave  fields.   During  periods  of  relatively  steady  winds  a  quasi-equilibrium 
balance  is  maintained,  but  this  is  easily  broken  by  changes  in  surface 
current  or  by  wind  velocity  shifts.   A  large  portion  of  the  momentum 
flux  to  the  surface  occurs  during  the  intermittent,  large  magnitude, 
strongly  coupled  bursts;   and  most  of  this  flux  goes  directly  into  high 
frequency  waves.   The  microscales  play  a  dominant  role  both  by  being  the 
active  recipients  in  the  momentum  transfer  and  by  assisting  the  growth  of 
lower  frequency  components  by  stress-strain  coupling,  by  acting  as  roughness 
or  drag  elements,  and  by  weakly  resonant  interactions. 

Many  questions  concerning  the  physics  of  wind-wave  generation 
remain  unanswered.   This  study  has  demonstrated  the  importance  of  inter- 
raittency,  strong  nonlinear  coupling,  and  microscales  in  momentum  transfer, 
but  has  given  no  indication  of  the  percent  of  the  total  transfer  attribu- 
table to  these  features.   Further,  the  effects  but  not  the  physics  of 
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energy  transfer  within  the  wave  field  have  been  described. 

Meaningful  future  experimentation  must  provide  a  better  description 
of  the  horizontal  and  vertical  structure  of  Reynolds  stress  in  air  and 
water;   quantitative  measurement  of  wave  growth  during  intermittent 
generation  and  its  relation  to  the  development  of  the  wave  field;   a  sta- 
tistical description  of  the  structure  and  spacing  of  intermittencies  and 
their  relationship  to  waves  and  wave  groups;   simultaneous  measutements  of 
high  and  low  frequency  wave  energy  and  their  spatial  and  time  relation- 
ships;  a  detailed  examination  of  the  air-sea  coupling  at  the  anomalous 
critical  wind  speeds;   and  controlled  determination  of  the  effects  of 
high  frequency  degrading  slicks. 

Open  ocean  measurements  of  velocity  and  turbulence  in  the  troughs  of 
major  waves  are  still  unavailable;   the  role  of  breaking  and  the  effect 
of  the  resulting  vertical  vortices  on  spectral  transfer  is  poorly  under- 
stood;  and  the  effects  of  wind-surface  current  coupling  and  the  underwater 
boundary  layer  on  wave  generation  have  still  been  inadequately  described. 
Future  studies  must  consider  the  importance  of  time  history  and  inter- 
mittency  during  generation,  and  most  progress  will  be  made  through 
detailed,  near  naviface  measurements. 
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APPENDIX  I.  INSTRUMENTS  AND  CALIBRATION 

A.  Sites  and  Platforms. 

Buoy  Experiment,  Instruments  were  mounted  on  the  M.I.T.  Air-Sea 
Interaction  Buoy  (Mollo-Christensen  and  Dorman, (1971)).  The  site  was  at 
the  southwestern  end  of  Vineyard  Sound  (*H  20'  ^0"  N,  70°  55*  ^5"  W). 
The  buoy  is  shown  in  Fig.  AI.l.  The  site  is  8.9  km  south  of  Cuttyhunk 
Harbour,  7.8  km  west  of  Gay  Head  Light,  Martha's  Vineyard,  and  approxi- 
mately 50  km  ENE  of  Block  Island.  The  bottom  in  the  vicinity  is  very 
fine  sand  and  mud,  depth  31.1  m  at  MLLW,  and  mean  tidal  range  less  than 
1  m.  The  bottom  is  nearly  flat  to  the  south  and  southwest,  and  shoals 
slowly  to  the  north  and  northeast. 

The  experiment  was  conducted  from  July  to  November,  1970.  Winds  in 
the  area  are  dominated  by  a  land-sea  breeze  regime  during  the  summer,  and 
are  predominantly  southwesterly  in  the  fall.  Tidal  currents  are  primarily 
semidiurnal  with  a  clockwise  ellipse.  They  are  strongly  affected  by  the 
winds.  Recorded  current  maxima  were  around  1.5  m/sec  with  a  typical  slack 
duration  of  20  minutes  and  nearly  equal  duration  of  flood  (NE)  and  ebb 

(sw). 

A  shore  support  and  data  collection  facility  was  housed  in  the  Coast 
Guard  Station  (deactivated)  on  Cuttyhunk.  The  facility  had  line-of- sight 
communication  with,  and  visual  observation  of  the  buoy.  Ship  support  was 
provided  by  the  R/v  R.  R.  Shrock  (M.I.T.'s  oceanographic  research  vessel) 
and  by  project  owned  and  local  small craft.  The  Shrock  was  equipped  with  an 
electronic  bathythermograph  used  to  monitor  thermal  structure  during  ex- 
periments. A  small  weather  station  was  maintained  on  Cuttyhunk. 

The  buoy  was  slack  moored  (scope  of  2)  and  consequently  faced  into  the 
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tide  (one  notable  exception  being  runs  109-159  when  continuous  high  wind 
reversed  the  surface  flow  and  kept  the  buoy  pointing  southwest).  Torque 
balance  to  counter  the  drag  of  the  underwater  instrument  tower  was  accom- 
plished by  angular  preadjustment  of  the  mooring  bridle.  The  bridle  was 
placed  near  the  center  of  lateral  resistance  so  the  buoy  and  attached 
instruments  had  <  1°  mean  tilt  during  all  current  conditions. 

The  buoy  and  its  small  motions  had  little  effect  on  the  data.  This  may 
be  shown  by  the  following  considerations.  Static  stability  was  high  due  to 
the  k,27  m  metacentric  height  and  small  surface  peretrating  cross-section. 
A  vertical  load  of  100  lbs,  10*  from  the  centerline  at  platform  height, 
caused  only  0.23°  tilt.  The  large  virtual  mass  and  damper  plate  resulted 
in  a  critically  damped  natural  heave  period  of  21.3  seconds,  with  an  am- 
plitude response  for  6  second  waves  of  A  ==  (Buoy  Amplitude)/ (Sea  Surface 
Amplitude)  =  0.07^-.  Yaw  was  determined  by  field  experience  to  be  of  sig- 
nificance only  during  the  hour  or  so  at  and  around  slack.  Fluctuating 
tilt  was  measured  with  a  pendulum  pot  and  showed  a  maximum  response  at 
around  10  seconds  with  maximum  tilts  of  1.5  •  Instruments  were  positioned 
30°  f orward  of  the  normal  to  the  mean  current  to  minimize  interference 
calculated  on  the  basis  of  potential  flow.  Crude  estimates  on  this  basis 
indicated  that  velocity  distortion  (AU/u)  was  0.028  for  the  minimeters, 
0.009  for  the  anemometers.  The  estimates  are  of  course  worst  near  the  navi- 
face,  in  regions  of  high  shear.  Wave  gauges  were  placed  on  both  outriggers 
and  their  outputs  compared  to  eliminate  buoy  interference  during  adverse 
wave  directions. 

The  buoy  was  equipped  with  a  thermoelectric  generator  to  allow  unat- 
tended operation  for  up  to  2  weeks  under  mean  condition  run  schedules 
(12  minutes  eash  hour),  or  four  hours  of  continuous  operation  with  over- 
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CATAMARAN    and 
COUPLING  ARRANGEMENT 

Figure  AI.2. 
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night  recharge.  Primary  power  source  was  a  set  of  eight,  12  volt  120  amp- 
hour  nicad  batteries.  The  buoy  has  been  in  use  for  two  years  and  has 
operated  successfully  in  winds  over  20  m/sec. 

Catamaran  Experiment.  This  experiment  was  conducted  in  March,  1971  > 
from  the  dock  of  the  H.I.T.  Yacht  Club  boathouse.  Bulk  atmospheric  condi- 
tions were  generally  stable  (see  Table  1,  Section  II).  Data  were  useful 
only  when  winds  were  in  the  southern  quadrants,  such  that  the  fetch  was 
over  300  m. 

A  diagram  of  the  catamaran  is  shown  in  Fig.  AI.2.  It  was  designed  to 
support  up  to  1.5  kg  of  instrumentation  in  a  nearly  constant  position  rela- 
tive to  the  instantaneous  surface  components  of  f  <1.2  hz.  Heave  decoupling 
and  minimization  of  mean  instrumental  tilt  were  accomplished  by  the  use  of 
two  vertical  supports,  each  sliding  through  two  delrin  bearings  which  were 
separated  by  0.3  m.   and  affixed  to  a  carefully  leveled  support  structure. 
Pitch  and  roll  bearings  were  located  respectively  at  the  attachment  point 
of  the  fore  and  aft  crossbars  to  the  pontoons,  and  at  the  vertical  support- 
crossbar  attachment.  This  arrangement  gave  the  catamaran  a  60  cm  vertical 
range  and  allowed  pitch  and  roll  angles  of  up  to  20°  with  less  than  1.5° 
tilt  of  the  sensors.  Construction  materials  were  polystyrene  foam  and  .020"^ 
6061  aluminum  sheeting  and  tubing.  Total  weight  was  3*2  kg,  draft  5  cm, 
and  design  natural  frequency  1.5  hz. 

Instruments  were  attached  to  the  forward  support  tube  from  15  cm  below 
to  30  cm  above  the  surface.  They  were  positioned  within  5  cm  of  the  bow 
line  to  minimize  flow  obstruction.  Additional  sensors  were  mounted  on  the 
support  structure.  Cables  were  floated  aft  of  the  cat  to  lessen  weight  and 
permit  directional  control  for  a  full  3^0°.  Accelerometry  tests  were  con- 
ducted with  full  instrumental  complement  during  most  runs.  They  showed 
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that  the  cat  closely  followed  waves  of  f <  1.3  hz,  and  had  negligible  re- 
sponse to  those  over  1.85  hz.  Vlaves  excited  by  catamaran  motion  caused 
extraneous  spectral  peaks  at  50  and  70  hz  but  had  little  effect  on  frequen- 
cies of  interest. 

B.  Instruments,  l^st  of  the  instruments  were  designed  and  constructed 
at  M.I.T.  and  have  seen  several  years  of  use  and  development.  In  addition 
to  precision  and  accuracy,  a  major  design  criterion  has  been  systematiza- 
tion  of  the  entire  array  of  transducers  to  insure  that  data  formats  will 
be  compatible  and  easy  to  process.  Hardware  and  electronics  will  be  des- 
cribed in  Project  reports  now  in  preparation,  and  also  have  been  discussed 
at  length  in  other  reports  from  the  Project  (Seesholtz  (1968);  Ruggles  (19&9) 
Von  Zweck  (1969)).  The  essential  characteristics  and  calibration  results 
are  as  follows. 

Velocity  Sensors.  Profile  velocities  were  measured  with  a  set  of 
5  matched  cup  anemometers  in  the  air  and  3  ducted  propeller  'minimeters* 
in  the  water.  Eoth  instrument  types  use  chopper  wheels  to  interrupt  a  light 
beam  that  is  sensed  by  photoresistors.  Static  and  dynamic  anemometer  cali- 
brations were  conducted  in  the  M.I.T.  low  turbulence  wind  tunnel.  The 
results  indicated  that  the  instruments  had  a  distance  constant  of  101  cm 
and  a  static  linear  calibration  of  U(cm/sec)=2.325x(pulses/min)  +  62  in  the 
speed  range  0.8  to  3°  m/sec.  A  repeat  calibration  at  the  end  of  the  summer 
showed  deviations  from  initial  values  to  be  less  than  2$.  Directional  tests 
showed  that  the  units  had  little  sensitivity  to  vertical  components.  The 
anemometers  are  dynamically  nonlinear  (Wyngaard,  et  al  (1971))  due  to  a 
faster  response  to  wind  increases  than  to  decreases.  Overestimates  in  gusty 
winds  were  roughly  6$.  Ilinimeters  were  calibrated  in  a  tow  tank  at  V/oods 
Hole  Oc^anographic  Institution,  at  -needs  of  8  to  105  cm/sec.  The  static 
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response  was  linear,  U(cm/sec)  =  0.9111x(pulses/sec)  +  2.5.  Distance 
constant  was  8  en.  Directional  response  (Fig.  AI.3)  was  best  fit  by  a 
cos©  +  K  sinO   curve  with  K  differing  between  instruments. 

Velocity  measurements  from  the  catamaran  were  made  with  DISA  type  55A81 
industrially  coated  hot  film  probes  operated  as  constant  temperature  ane- 
mometers with  overheat  ratios  of  1.55A  in  air  and  1.08/1  in  water.  Static 
calibrations  were  simultaneously  conducted  with  all  5  probes  in  the  II. I. T. 
low  turbulence  wind  tunnel.  The  response  was  expressed  in  terms  of  a  best 
fit  (in  the  least  mean  squares  sense)  quartic  polynomial.  Directional  re- 
sponse (Fig.  AI.3)  was  similarly  fit  by  the  linear  relation 
U(cm/sec)  =  Q(1.0  -  0.002l(e  -25)),  where  Q  =  actual  velocity  and  <9  - 
attack  angle,  in  the  range  25° <  ©  ^  80°  (correlation  =  0.92).  Dynamic 
response  (Rasmussen  and  Bellhouse  (I968))  was  satisfactory  for  the  frequen- 
cies encountered  during  cat's  paws  and  in  fact  compensated  for  the  slight 
velocity  overestimate  at  high  attack  angles  inherent  in  the  directional  fit. 
Wave  Sensors.  Uave  gauges  were  capacitance  type;  they  were  design- 
ed and  built  at  M.I.T.  and  have  provided  satisfactory  results  in  the  field 
and  laboratory  for  4  years,  at  wave  heights  up  to  5  m  and  frequencies  up 
to  100  hz.  The  output  is  a  variable  frequency  that  is  created  by  beating 
the  frequency  from  an  oscillator  controlled  by  the  sensor  capacitance 
against  a  standard.  The  wave  gauges  used  on  the  buoy  had  a  4.6  m  long,  .95 
cm  diameter  stainless  steel  wire  core  steering  cable  sensor,  and  a  30.1  hz 
standard  oscillator.  Full  range  bandwidth  was  400  hz  with  resultant  re- 
sponse of  3.2  hz/cm.  The  capillary  wave  gauge  for  the  catamaran  had  a  15 
cm  long,  0.076  cm  diameter  teflon  coated  sensor  and  a  200  hz  standard  oscil- 
lator. Full  range  variation  was  1.4  to  2.2  khz,  or  amplitude  response  of 
3.88  hz/mrn.  Sensors  xvere  frequently  cleansed  with  alcohol  to  reduce  film 
effects. 
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Temperature  Sensors.  Two  types  of  thermistors,  both- with  15  kil 
resistance  at  25°C  were  used  to  permit  interchangability  in  the  bridge 
circuit  and  provide  roughly  equal  time  constants  in  air  and  water.  Water 
temperature  was  measured  with  matched  Fenwall  KkS>6   oceanographic  isocurve 
submini  thermistors.  Each  unit  had  two,  O.O76  cm  diameter  bead-in-glass 
tips.  Air  to  still  water  time  constant  was  70  ms.  The  Fenwalls  were  epoxy 
mounted  in  1.3  cm  diameter  stainless  steel  tubing  fitted  with  retractable, 
ventilated  shields.  The  shields  were  kept  forward  during  the  mean  condition 
runs  (removed  during  catamaran  runs),  and  decreased  the  time  constant  in 
moving  water  to  approximately. 2  seconds.  Air  thermistors  on  the  catamaran 
were  VECO  model  42A70,  .  0508  cm  diameter  bead-in-glass  probes  mounted  in 
.63  cm  diameter  tapered  wooden  dowling.  Their  air  to  still  water  time  con- 
stant was  35  ms.  Using  the  result  of  Aagaard  (1964)  for  similar  units, 
amplitude  response  at  nyquist  frequency  (125  hz)  for  the  VECOs  was  0.44. 

The  differential  bridge  circuit  was  designed  at  M.I.T.  It  was  used  in 
the  absolute  mode  with  a  6°  -  26°  C  window  for  mean  condition  runs;  accur- 
acy and  precision  with  both  thermistor  types  was  *  0.1°C.  Measurements 
from  the  catamaran  used  the  relative  mode  (thermistors  balanced  against 
precision  resistors  selected  on  the  basis  of  mean  temperature)  with  a  5°C 
window,  and  had  a  precision  of  i  0.02°C.  A  bomb  calorimeter  was  employed 
as  field  and  calibration  temperature  standard. 

Additional  instruments  included  a  pendulum  pot  for  measurement  of  buoy 
tilt,  a  tri-axial  strain  gauge  accelerometer  for  measurement  of  catamaran 
motion,  and  an  electronic  bathythermograph  connected  to  an  x-y  plotter 
on  the  Shrock. 

C.  Mounting. 

Buoy.  Cup  anemometers  were  attached  to  aluminum  channel  which  in 
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turn  was  bolted  to  the  meteorological  tower.  Cups  were  slightly  more  than 
1  m  from  the  tower.  Heights  (for  runs  67   -  159)  were  2.4,  3.7,  5.0,  7.5, 
and  11.4  m.  Minimeters  were  attached  to  carts  ;jhich  ran  on  one  set  of 
tracks  on  the  underwater  tower.  They  could  be  depth-positioned  by  hand  from 
the  buoy  outrigger.  Sensors  were  over  50  cm  from  the  tower  which  in  turn 
was  2.5  m  from  the  main  spar.  Mean  depths  (runs  109  -  159)  were  -1.52,  -2.74, 
and  -4.7  m.  Thermistors  were  mounted  on  the  arms  of  the  minimeters  and  an- 
emometers. The  primary  wave  gauge  was  located  directly  above  the  minimeters; 
the  secondary  was  on  the  opposing  outrigger.  Hot  film  and  thermistor  elec- 
tronics were  housed  in  a  wooden  box  located  between  the  braces  of  the 
tower  support,  and  the  main  electronics  package  was  in  a  compartment  in 
the  main  spar,  directly  under  the  personnel  platform. 

Catamaran .  The  accelerometer  and  wave  gauge  electronics  were  set 

6 

in  the  forward  cross-brace.  Orthogonal  hot  films  and  thermistors  were  at- 
tached to  the  forward  vertical  support  with  delrin  mounts.  The  wave  gauge 
sensor  was  stretched  between  two  delrin  braces  on  the  forward  support. 
Cross-stream  separation  between  the  wave  gauge  and  hot  film  tips  was  0.95 
cm;  between  the  wave  gauge  (air)  or  hot  films  (water)  and  thermistors  1.3  cm 
(see  Table  1,  Section  II).  Downstream  sepatation  was  <   1  mm.  Orthogonal 
positioning  of  hot  films  was  insured  by  pre-fixed  mounting  and  was  checked 
by  measurements  on  projected  photographs.  Separation  between  hot  film  tips 
was  3  roro* 
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APPENDIX  II.  DATA  COLLECTION  AND  ANALYSIS 

Both  hard-  and  software  packages  will  be  described  in  detail  in  forth- 
coning  reports  of  the  M.I.T.  Buoy  Project.  This  section  will  describe  the 
functional  nature  of  the  system  and  outline  the  major  subroutines  and  alo- 
grithms  used  in  the  analysis.  The  major  components  of  the  system  are 
depicted  in  Figs.  AII.l  and  All. 2. 

A.  Collection  and  Format. 

Buoy  System.  Information  from  the  transducers  passes  through  con- 
ditioning amplifiers  designed  to  produce  0-5  volts  variable  D.C.  output. 
This  output  drives  voltage  controlled  oscillators  set  on  standard  I RIG 
channels  3-15  and  17.  The  channels  were  matched  against  anticipated  data 
bandwidth  to  insure  at  least  a  5:1  bandwidth  ratio.  Accelerometer  pulses 
(and  minimeter  pylses  during  some  runs)  were  sampled  at  100  hz  and  combined 
on  a  single  IRIG  channel.  A  standard  25  khz  signal  was  also  transmitted 
and  later  used  for  synchronization  and  pulse  separation  during  breakout. 
Signals  from  up  to  22  separate  transducers  could  be  accomodated  simultan- 
eously. Signals  from  all  IRIG  channels  were  nixed  and  then  used  to  fre- 
quency modulate  the  transmission  signal  (243.8  mhz).  The  major  components 
of  the  buoy  data  handling  package  were  modified  from  an  Atlas  Ilissle  Telem- 
etry System  (USAF  excess  equipment). 

Data  received  at  the  base  station  on  Cuttyhunk  was  recorded  in  analog, 
direct  node  on  a  single  channel  of  ij-H  nagnctic  tape  on  a  Hewlett  Packard 
Model  39^03  (HP3960B)  4  track  recorder.  A  second. channel  was  reserved 
for  voice  narration,  weather  log,  and  time  signal.  Both  sides  of  2400* 
tapes  were  used  for  mean  condition  runs;  one  tape  held  four,  12  minute 
hourly  runs. 
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The  only  analysis  attempted  at  the  site  was  demodulation • and  monitoring 
of  selected  channels  to  insure  signal  quality  and  continuity.  The  buoy  was 
turned  on  and  off  remotely  by  a  modified  Heathkit  radio  control  set  which 
activated  a  micro switch  in  the  buoy  electronics  compartment. 

Catamaran  System.  Signals  from  catamaran  transducers  were  recorded 
on  a  Princeton  Instruments  Model  6100  (PI6100)  8  track  recorder.  Hot  film, 
thermistor,  and  accelerometer  data  were  passed  through  conditioning  ampli- 
fiers to  optimize  voltage  bias  and  gain,  and  then  recorded  in  FII  mode. 
Capillary  wave  gauge  signals,  in  frequency  format,  were  recorded  in  direct 
mode.  One  channel  was  reserved  for  voice  narration;  this  was  subsequently 
used  to  help  select  segments  for  digitization.  The  incoming  data  was  con- 
tinuously monitored  to  insure  signal  quality  and  continuity.  2*4-00*  tapes 
allowed  over  2  hours  of  continuous  recording  per  run. 

B.  Analysis. 

General.  A  combination  of  digital  and  analog  analysis  was  employed. 
Analog  processing  with  strip  chart  and  x-y  recorder  display  was  used  for 
initial  evaluation.  Energy  content  in  selected  frequency  bands  was  estimated 
by  filtering  and  squaring  the  signals  prior  to  plotting.  Analog  spectral 
analysis  of  potentially  interesting  events  was  performed  with  the  assistance 
of  a  PI  correlator  and  Fourier  analyzer  and  the  x-y  plotter.  The  initial 
processing  familiarized  the  user  with  the  data  and  permitted  early  estimates 
of  fluctuation  magnitude  for  later  use  in  digital  processing. 

Three  digital  computers  were  used  for  secondary  and  final  analysis. 
Digitization,  preprocessing,  and  simple  analyses  such  as  pulse  counts  and 
calibration  checks  were  performed  on  the  Project's  Honeywell  H-316  compu- 
ter. This  unit  has  an  8  k  memory  and  is  equipped  with  a  16  channel  Raytheon 
' miniver ter'  analog  to  digital  (a/d)  converter  Model  MADC  12A.  The  miniver- 


162. 


ter  outputs  12  bit  a/d  words  for  a  full  range  of  i  10  volts;. least  bit  pre- 
cision is  5  rov.  Programming  was  in  Honeywell's  DAP-16  Assembler  level 
language.  Output  from  the  Honeywell  was  in  two  forms:  the  digitized  data 
in  bit-string,  a/d  word  format  on  7  track  magnetic  tape,  and  pulse  count 
and  calibration  information  on  either  the  ASR-33  teletypewriter  or  paper 
tape.  The  latter  information  was  transferred  to  punched  cards. 

A  PDP-7  computer  (at  Linclon  Laboratory)  equipped  with  miltiple  tape 
drives  and  cathode  ray  tube  (CRT)  display  was  used  to  examine  in  detail 
the  tapes  produced  by  the  H-316.  The  bit-string  data  was  first  converted  to 
geophysically  meaningful  units;  statistics  were  calculated;  values  from  or- 
thogonally mounted  hot  films  were  converted  to  U  and  W;  and  the  resulting 
information  was  then  rerecorded  in  a  format  standard  to  the  facility.  The 
resultant  data  were  displayed  on  the  CRT,  limited  filtering  and  spectral 
analysis  performed,  and  hard  copy  graphical  output  obtained  for  further 
study.  Programming  was  done  in  Fortran  IV  by  Mr.  B.  Laird. 

Final  analysis  was  done  on  the  IBM  360  at  the  H.I.T.  Information  Pro- 
cessing Facilit}r  (IPC).  Programs  were  written  in  PL/l  to  reduce  l/O  prob- 
lems, and  extensive  use  was  made  of  a  Stromberg-Carlson  4020  (SC-4020)  plot- 
ter for  display  of  time  series,  probability  distributions,  and  spectra. 

The  data  handling  system  is  capable  of  statistical,  time  history,  and 
frequency  domain  analysis  of  up  to  16  channels  of  input  data  simultaneously. 
Analog  preprocessing  capabilities  include  rerecording  at  a  choice  of  speeds 
and  formats,  filtering,  integration,  squaring,  orthogonal  decomposition, 
amplification,  and  pulse  separation,  as  well  as  totally  analog  computation 
and  display  of  time  history,  correlations,  and  auto-and  cross-spectra.  The 
H-316  allows  digitization  at  up  to  250  hz. 

Buoy  Data.  The  mixed  signals  from  the  HP396OB  were  separated  in  a 
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set  of  8  demodulators  (from  the  Atlas  Tlissle  System).  Sampled  anemometer 
signals  were  separated  using  ring  counters  and  the  25  khz  clock  signal  for 
timing.  Sections  of  the  data  were  filtered,  squared,  amplified  and  display- 
ed on  a  ^-track  Sanborn  recorder  to  provide  an  initial  overview.  Tiro  sep- 
arate digitizing  passes  were  then  made.  The  first  yielded  velocity  profile 
information  from  a  pulse  counting  program  on  the  H-316.  Twenty,  30  secon  d 
counts  from  10  minutes  of  each  mean  condition  run  were  made,  transferred  to 
punched  cards,  and  run  through  a  two  part  program  on  the  IBM  3&0.  The  first 
digital  program  computed  profile  parameters  on  the  basis  of  the  logar- 
ithmic profile  model.  The  second  program  used  punched  card  output  from  the 
first  to  compute  single  and  joint  histograms  and  statistical  information, 
displayed  on  SC-^020  plots.  Velocity  values  from  selected  runs  (with  aver- 
aging times  as  short  as  5  seconds)  were  hand  plotted  in  profile,  isotach, 
and  vertical  shear  format  for  comparison  with  the  strip  chart  analog  data. 

In  the  second  digitizing  pass,  the  runs  selected  for  time  history  anal- 
ysis were  rerecorded  on  the  PI6100  in  demodulated  form.  The  data  was  fil- 
tered to  avoid  aliasing  (velocity  values  at  2.5  hz  and  wave  height  at  5  hz 
using  matched  Rockwell  I.odel  1010F  operational  filters),  amplified  to  the 
full  range  of  the  miniverter,  and  digitized  at  10  hz.  Calibration  during 
this  preprocessing  was  maintained  by  passing  sinusoidal  and  D.C.  signals 
through  the  entire  system.  The  calibration  data  were  displayed  on  the  ASR-33 
and  transferred  to  cards  for  use  in  the  main  analysis.  Additional  calibra- 
tion checks  during  digital  processing  were  made  by  comparison  of  pulse 
count  and  digitized  velocity  values.  The  time  signal  initially  recorded  on 
the  voice  channel  of  the  HP3960B  was  used  to  synchronize  the  two  passes. 
Final  analysis  used  the  software  package  described  below. 

Catamaran  Data.  All  of  the  cat's  paw  data  was  pre-viewed  using 
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analog  techniques.  Selected  portions  were  then  digitized  at  250  hz;  the 
data  had  been  low-pass  filtered  at  100  hz  using  the  internal  PI6100  filters. 
The  voice  narration  was  used  to  aid  in  distinguishing  the  presence  of  gusts 
until  their  imprint  on  the  analog  signals  became  readily  recognizable. 
Nine  tapes  containing  some  2.5  x  10'  data  points  from  each  of  8  channels 
were  obtained.  These  tapes  contain  information  from  over  75   well  identified 
catfs  paws  as  well  as  a  large  amount  of  background  data.  The  tapes  were 
examined  on  the  L.L.  PDP-7  and  segments  selected  for  final  analysis. 

IE!  360  Software.  The  analysis  package  was  designed  to  permit 
rapid,  standardized  analysis  of  time  series  information  sampled  at  equal 
intervals.  Flexibility  and  user  control  were  major  design  criteria.  The 
package  can  handle  up  to  16  channels  of  information  simultaneously,  and 
allows  for  the  inclusion  01  two  additional  channels  of  computed  values  such 
as  Reynolds  stress  and  buoyancy  fluctuations.  In  addition,  the  initial 
information  can  be  replaced  during  analysis  by  other  computed  values.  The 
package  may  be  of  use  to  other  investigators,  and  is  consequently  briefly 
described  herein. 

Program  flow  is  diagrammed  in  Fig.  All. 2.  It  is  highly  modular  and 
designed  to  handle  a  wide  variety  of  time  series  information,  with  a  mini- 
mal number  of  routines  specifically  restricted  to  air-sea  interaction 
studies.  Flexibility  is  maintained  by  user  option  control  from  punched  card 
parameters.  There  is  provision  for  a  user  supplied  routine  (OTHER)  which 
can  be  reached  from  both  time  and  frequency  domains.  Time  domain  analysis 
is  performed  first;  in  the  frequency  domain  segment  the  original  data 
(stored  in  a  three-dimensional  matrix)  is  replaced  by  its  fourier  transform 
coefficients.  The  major  features  and  alogrithms  are  as  follows;  details  will 
be  given  in  a  forthcoming  Project  report. 
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COIiTRL  handles  tape  positioning,  dimensions  the  major  data  matrix, 
selects  plot/no  plot  mode,  and  inputs  data  from  the  tape.  This  is  the  only 
routine  that  is  specifically  linked  to  the  H-316"  7  track  magnetic  tape 
format. 

DOIT  is  called  by  COMTRL  upon  completion  of  input.  This  subroutine 
does  little  computing,  but  determines  which  other  routines  will  be  called 
and  the  sequence  of  their  use.  It  does  contain  the  hanning  and  stress  for- 
mation routines. 

REALVAL  performs,  with  operator  option,  the  following  functions: 
conversion  of  a/d  words  into  voltage;  conversion  from  voltage  to  real  values 
with  up  to  a  quartic  fit  (required  for  linearization  of  the  hot  films);  de- 
termination of  trend,  extrema,  and  the  1st  k   central  moments;  computation 
and  plotting  of  a  histogram  with  user-supplied  number  of  intervals  and 
limits  set  either  to  computed  extrema  or  user-supplied;  adjustment  of  values 
to  a  calibration  standard  (used  for  comparison  of  pulse  count  and  digital 
velocities)  and/or  removal  of  mean  (supplied  or  computed)  and/or  trend, 
■'ean  and  trend  may  also  ba  removed  in  DOIT  immediately  prior  to  the  fre- 
quency domain  segmant. 

FILTER  operates  by  fast  fourier  transforming  the  data,  modifying  or 
deleting  the  selected  frequency  components,  and  then  retransforming.  The 
user  supplies  the  modification  factors  so  that  effectively  any  desired  win- 
dow may  be  employed.  A  standard  cosine  taper  of  adjustable  length  is  inclu- 
ded. This  routine  can  be  called  from  up  to  3  different  locations  in  DOIT 
so  that  the  data  may  be  used  in  both  filtered  and  unfiltered  form. 

ORTHOG  computes  horizontal  (U)  and  vertical  (W)  velocities  from  ortho- 
gonally mounted  instruments.  The  decomposition  is  done  in  two  steps;  first 
an  estimate  of  U  and  W  is  obtained  from  the  known  directional  response  and 
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velocity  readings  of  the  two  instruments,  and  then  an  alignment  correction 
is  made  to  compensate  for  mean  tilt  or  oscillatory  fluctuations  of  the 
probe  holder.  The  background  and  details  of  this  procedure  were  discussed 
at  length  by  Ruggles  (1969).  There  is  a  choice  of  either  a  cos©  +  K  sin© 
or  a  linear  directional  fit  to  the  data.  The  alignment  correction  assumes 
that  W  =  0  over  a  number  of  data  points  selected  by  the  user,  and  that  the 
instruments  are  truly  orthogonally  mounted.  A  known  mean  horizontal  vel- 
ocity, if  available,  can  be  used  for  additional  correction.  The  choice  of 
the  averaging  interval  for  VJ  =  0  is  critical,  and  is  made  by  a  considera- 
tion of  the  dynamics  of  the  system  and  by  post-facto  examination  of  plots 
of  the  data  to  insure  continuity  of  values  at  interval  end  points.  For  the 
catamaran  data,  60  points  -were  used  for  air  hot  films,  250  for  water 
hot  films. 

PROFIL  computes  parameters  of  the  assumed  logarithmic  profile  model; 
effectively,  it  compares  U(z)  to  (U*/k)  ln(z/z0)  by  a  best  fit  in  the  least 
mean  squares  sense,  taking  the  instrument  heights  as  known  (In  previous 
work  from  this  project,  e.g.  Ruggles  (1969),  velocity  values  were  taken  as 
known,  z  as  the  variable.  In  the  initial  stages  of  this  analysis  both  methods 
were  used,  and  for  profiles  showing  correlations  with  the  model  better  than 
0.9^  the  resulting  differences  in  U*  and  z0  were  always  less  than  5$»)«  If 
Reynolds  stress  values  are  available,  the  production  term  of  the  horizontal 


component  of  turbulent  kinetic  energy  equation  (-u*w'  dU/dz)  can  be  calcu- 
lated. The  separate  program  to  handle  the  pulse  count  data  use's  similar 
alogrithms. 

JOINTP  computes  and  plots  joint  histograms,  and  also  provides  one-way 
histograms  and  means-on-interval  of  the  input  variables. 

TIIIEPLT  plots  variables  against  time,  with  user  controlled  or  computed 
scaling  and  number  of  strips  per  SC-4-020  frame.  It  can  be  accessed  both 
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before  and  after  filtering.  This  subroutine  also  plots  histograms. 

F0URT2  is  a  USASI  Basic  Fortran  version  of  the  Cooley-Tukey  fast  fourier 
transform,  written  by  II.  Brenner.  It  can  handle  real  or  complex  data  in 
multidimensional,  arbitrarily  sized  arrays.  End  values  are  dropped  and 
tapered  in  DOIT  to  provide  a  window  prior  to  transforming;  in  this  analysis 
approximately  IO^j  of  the  data  was  so  treated  to  limit  Gibbs  effect.  The 
returned  coefficients  are  normalized  and  hanned  in  DOIT,  then  a  density 
factor  is  applied  to  result  in  the  power  spectral  density  estimate.  The 
hanned  estimate  has  a  "^  distribution  with  6  dof ;  restriction  of  averaging 
to  3  values  was  necessary  in  the  catamaran  study  to  preserve  detailed 
resolution  in  view  of  the  high  sampling  rate  and  relatively  few  points 
available  for  analysis  (core  and  averaging  limitations). 

XSPECTR  utilizes  the  normalized  fourier  transform  coefficients  to 
determine  cross-spectral  densities  (Co  and  Quad  spectra).  These  are  then 
used  to  compute  coherence  magnitude  ( V  )  and  phase  (©  ).  The  averaging 
interval  used  to  compute  coherence  is  selectable  by  the  user;  in  this  study 
9  points  were  used  for  samples  under  500  points,  15  for  the  5-roinute  mean 
condition  run  seg merits,  and  25  for  the  10-second  catamaran  time  periods. 
The  results  were  then  hanned  for  the  longer  samples. 

SPECOUT  prints  and  plots  spectral  data.  Auto spectra  were  plotted  on  a 
log (magnitude)  linear  (frequency)  grid  (log-log  is  also  available).  Magni- 
tude scaling  is  based  on  spectral  extrema  (up  to  10  decades  can  be  accomo- 
dated), but  the  user  controls  the  frequency  range.  Coherence  and  phase 
were  plotted  from  the  cross-spectral  analyses;  the  user  again  controls  the 
frequency  range  and  number  of  plots. 

OTHER  is  the  user  supplied  routine.  It  can  be  accessed  from  two  loca- 
tions in  the  time  domain  and  at  the  end  of  the  frequency  domain  computations 
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in  DOIT.  Two  major  versions  were  used  in  this  study,  one  each  for  the  buoy 
and  catamaran  experiments.  The  buoy  version  computed  mean  energy  density, 
both  directly  from  wave  height  data  and  by  spectral  accumulation,  and  wave 
momentum  at  individual  frequency  components  and  by  considering  the  energy 
to  be  concentrated  in  waves  at  the  frequency  of  energy  and  momentum  spec- 
tral maxima.  The  catamaran  version  had  four  segments;  one  computed  buoyancy 
fluctuations  and  determined  joint  histograms  between  buoyancy  and  other 
variables;  a  second  formed  the  natural  logarithm  of  squared  and  normalized 
velocity  derivatives  for  the  lognormality  study;  a  third  created  the  switch 
function  that  was  used  to  evaluate  intermittency;  and  the  fourth  contained 
a  band  pass  filter  for  rapid  examination  of  filtered  signals. 

Computational  Methods.  Let  x^,  i  =  1,  2,  ...  N  represent  the 
sampled  data  values  where  NAt  =  T  is  the  total  sample  length  and  l/At  =Af 
the  sampling  frequency.  Then  for  any  a/d  word  yi,  voltage 

v^  =  -10.0  +  yjj_  x  204. 8  from  the  known  characteristics  of  the  miniverter, 

I  i 

and    xi  =  ^  A 4  v*  in  geophysical  units  where  the  Aa*s   are  the  conver- 


sion  values  from  calibration.  Then 


M 


mean 


.th 


=  *i  =   *  *  u  21  *;  and 


N 


j   moment  =  -—  %   (*;  -  *\J         so  that 

variance  =   x;1  -  yz  =  &-z 

skewness  =   x?  -  3 7?  x  +1^   H  5^  (     and 

kurtosis  =   x?  -4x?  7  +  &***"-*?  *  15,  . 

Trend  (T)  is  found  from  a  best  fit  in  the  least  mean  squares  sense  to  the 
data.  Define  *r  -  T  «•  +■  /S   j  and  minimize 

6  a  (*;.-*x)  -  (*;-T  i  -/3s)7"  ^  where  both  Xj_  and  i  are  considered 

to  be  random  variates.  Then  the  expectation  of  0   is 
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vrhere  <**;  is  the  covariance  of  x^  and  i.  The. positive  definite  part  must 
=  0  for  minimum  S  ,  thus  /3  -    x  -  *V  \     . 

With  the  correlation  coefficient    /?  s -r.,   , 

adding  and  subtracting  ^ or/  ,    E  (SzN)  =  (Tc\  -  po^f  *■  (\  -p*) cr* 

so  that  for  minimum  o    f  =  f  ^  =  cr^j.     ancj  E(V>  =  f'-P)^7"  - 

*  °t     a': 

Thus,     -*- i  =  \ 

N   2-  L     _   lZ 


but 

i    = 

±LSJ             ,        E    il    - 

1- 

so  that 

T  - 

M 

o 

and  if  the  mean  and  trend  are  removed, 

For  the  calculation  of  profile  parameters,  the  assumed  model  is 

where  U*  is  friction  velocity,  K.  is  Von  Karmanfs  constant  (=  0.40)  and 
zQ  is  roughness  length.  Using  the  least  squares  method  developed  above, 
substituting  U^  for  x±9     S  e  U#/*i   for  T  ,  ln(z)  for  i  and 
(-U^/k  )  ln(z0)  for  /S      , 

r i   S  =  go-jucQ 


s 
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Further,  stress  ~c     =  DU*  where  D=  density  of  the  medium, 

U(iom\  5    Ute    =      ^    (3,  2.302?S>5     -  J*(*oV) 

for  velocity  in  cm/ sec  and  heights  in  cm,  and  drag  coefficient 

Figure  All. 3  shows  the  assumed  mounting  and  notation  for  orthogonal 

decomposition.  Assuming  that  the  hot 
films  have  a  linear  directional  response 
in  some  appropriate  range  of  <P    ,  let  the 
velocity  recorded  be  V^,  so  that 

where  Q^_  is  the  actual  velocity/  and  P 
and  V  are  determined  from  calibrations. 
For  the  DISA  hot  films,  V   =   25° ,  P  =  .0021 
and  the  relationship  holds  ( /°  -   0.92) 
in  the  attack  angle  range  25  €   &  *   80°. 

But   <£„  s  45°  +  0  and    4>A   =  AS°  -  Q 


a\\\ 

Hot        \$> 

Films  or      y 
M.M.      y^X 

B  /7^ 

/ 
/ 
/ 
/ 

•         \ 

'            1 

</>d=  45°  "    9 

Figure  All. 3 


so  that 


Qi 


V  ^  V 


9-  ■  ('-«*>('-&) 


2  ( i  -  lo?} 

Then  horizontal  (U)  and  vertical  (\J)   velocity  estimates  are  obtained  from 


W;  =  Q;  siv4  8; 


vJ;  =  Q;  ccs  6: 


For  minimeters,  the  directional  response  is  best  given  by 

in  the  range  0°  <  i>  <      80°,  with  the  value  of  P  varying  between  instru- 
ments. Substituting  for  <fc    in  terms  of  0'   and  then  in  terms  of  W  and  U, 
.£z/Va^V8A       y   _  jt  /  Va-  -v6a 
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To  correct  for  tilt,  assume  \J^   and  Ui  found  as  above  are  actually  at 
an  angle  «*;to  the  correct  axis.  Then  essentially,  letting  *U  and  °V  repre- 
sent the  corrected  values, 

Require  VJ  =  0;  then  if  a  known  U  is  available,  the  correction  to  lowest 
order  is  found  from 


5IN»<;     = 

<  M  > 

,         CO$=^;        = 

0 

;   nnt   kr 

lOT-jn,              Tii 

■Jo/.       —         

is  used.  The  averag- 

ing  times  are  selected  by  the  user  to  minimize  mean  or  fluctuating  tilt. 
Total  rms  error  in  computed  Reynolds  stress  using  this  lowest  order  tech- 
nique can  be  as  large  as  20$,  For  very  small  values  of  <fc>    ,  say  d>  <    10  , 
the  linear  hot  film  directional  response  fit  overestimates  V^  by  roughly 
1.5^.  Assuming  10$  rms  turbulence  levels  for  normal  atmospheric  conditions, 
the  probable  rms  deviation  of  <£  around  45  is  roughly  10  ;  thus  this  type 
of  fit  overestimates  only  extremely  high  turbulence  values.  This  is  somewhat 
compensated  for  by  the  frequency  response  of  the-  units,  but  still  tends 
to  emphasize  anomalously  high  stress  values  which  is  useful  for  detection. 

In  spectral  computations,  F0URT2  returns  the  real  and  imaginary  coef- 
ficients of  the  discretized  fourier  transform, 

V/here  t  =  j  At' and  s  =  -kAf,   F(«)  ■  *l*) +'<■/&  L*)  is  returned. 

The  coefficients  are  normalized  by  the  number  of  data  points,  e.g. 
a(^  =  ZotL4*)  /v4   ,  and  the  power  or  variance  spectral  estimate  is  formed: 
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Manning  is  applied: 

and  the  power  spectral  density  estimate  is  formed  by  multiplying  by  the 
density  factor,        *e  i4  6t  -T 

i.e.,     $,.«•  $0(»)--J>  *r   5„W-  Crffe^v^Crf)/^ 

where  f   is  the  Nyquist  or  folding  frequency  =  l/2At. 

Computation  of  the  cross- spectrum  (following  Wunsch,  (1970).)  is  based 
on  the  hypothesis  that  the  variables  are  related  by  the  convolution 

relationship,       y .  -  ][_    £  u 

Ji-  -«• 

for  all  realizations  of  x  and  y,  taking  x  and  y  as  random  variates  so  that 
S  is  a  known  transient.  Then 

where  cl    and  G  are  the  fourier  transforms  of  %   and  y.  Then 


For  the  hypothesized  convolution  relationship  to  have  meaning,  0(s)  must 
be  a  smoothly  varying  function  of  s,and  coherence  is  a  measure  of  this 
smoothness,  i.e.,  ut 

For  peaked  spectra  this  is  best  computed  in  a  form  which  includes  a  weigting 
by  the  power  at  the  various  frequencies,  i.e., 

where   denotes  the  complex  conjugate.  But 
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F„  l>)  ftJW  .  G^a  +  i  Q^  U>  ,  and   |f.  «|1G>(*|  -  (  $^  §   ft))* 


Coherence  magnitude  is  thus  defined  as 


i„   («>)  =     ; 


Using  p  and  x  to  denote  real  and  imaginary  parts  of  the  fourier  trans- 
form, the  alogrithm  is 

Similarly,  phase 

©  l^  --  -*  /  TftN-'  f 2 — 

M£-./,    l^W  J 

M  is  chosen  by  the  user  and  end  values  are  assumed  symmetrical  for  the 
averaging.  The  alogrithm  is  designed  so  that  positive  ©^IV)  means 
y  leads  x. 

Mean  wave  energy  density  is  computed  from  the  relationship 

t  -  rvi1 

where  y\      is  wave  height,  /»   density,  and  5  gravity.  But  ^  =  %7^ 
where   ^   is  the  autocorrelation  function  (not  computed);  in  turn, 

<*»  u>  00  I 

Both  direct  and  spectral  summation  techniques  were  used.  The  results  were 
always  within  1.5$  of  each  other  for  spectra  with  6  dof  or  more.  Departures 
of  20\j  with  mean  departure  roughly  5,j  occur ed  with  2  dof  (raw  or  unhanned 
short  spectral  estimates). 
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"omentum  was  computed  from  the  relationship 
v/here  c  is  linear  wave  celerity, 

Jk  =  wave  number  and  h  =  31  ™»  Thus, 

Mean  momentum  density  was  computed  both  by  summation  over  the  spectral 
components  and  by  assumption  that  E  was  primarily  attributable  to  the  wave 
of  maximum  variance  or  momentum,  i.e., 

m-  -£—       -    JL. 
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